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ABSTRACT 



MEASUREMENT OF CONTINUUM DIMUON PRODUCTION IN 
800-GEV/C PROTON-NUCLEON COLLISIONS 

BY 

JASON C. WEBB, B.S. 

Doctor of Philosophy 
New Mexico State University 
Las Cruces, New Mexico, 2002 
Dr. Vassili Papavassiliou, Chair 

Fermilab Experiment 866 has performed an absolute measurement of con- 
tinuum dimuon (Drell-Yan) cross sections in 800-GeV/c pp and pd interactions. 
Results differential in the mass, Feynman-x (xp) and transverse momenta (pr) 
of the dimuon pairs are reported. These results represent the most extensive 
study of the Drell-Yan process to date, based on a data sample of 175,000 dimuon 
events covering the widest range in kinematics yet achieved (4.2 < M < 16.85 
GeV and —0.05 < x F < 0.8) with this level of precision. The cross sections are 
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primarily sensitive to the magnitude and shape of the light antiquark distribu- 
tions (d(x) and u(x)) in the nucleon, but also provide important information on 
the valence quarks as well as the gluons. They are in good agreement with other 
existing proton-induced Drell-Yan experiments. There is also general agreement 
between the data and next-to-leading-order calculations based on various sets of 
parton distribution functions. Differences between data and theory are examined, 
and the potential impact of these data on future parameterizations of the parton 
distributions are discussed. 
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1 INTRODUCTION 



Investigation into the field of lepton-pair production in hadronic interactions be- 
gan with an experiment carried out at Brookhaven National Laboratory in 1970 
[Chr 70]. This experiment studied the reaction pXJ — > fi + fi~X for proton ener- 
gies between 22 and 29 GeV, resulting in pair masses between 1 and 6.7 GeV 1 . 
The data above ~ 3 GeV (where the J / tp resonance family would later be found) 
exhibited a rapidly decreasing continuum of muon pairs. The steeply falling na- 
ture of the cross section was explained later that year by Drell and Yan [Dre 70], 
who were interested in dilepton production as a possible application of the quark- 
parton model of hadron structure outside of deeply-inelastic scattering experi- 
ments. Their model of quark- ant iquark annihilation through the electromagnetic 
interaction, which has become known as the Drell- Yan process, was generally suc- 
cessful in describing the data and would only improve as our understanding of the 
strong interaction improved. 

The Drell- Yan process still remains an active area of experimental and the- 
oretical research some thirty years later. It has played a key role in developing 
the mathematical technology of perturbative Quantum Chromodynamics (QCD), 
being one of the first processes to be calculated to next-to-leading order 0(a s ), 

and remains one of the few processes to be calculated to next-to-next-to-leading 
1 Throughout this thesis, we will utilize "God-given" units, setting K = c = 1. 
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0{a 2 g). Experimentally it has provided a wealth of information about nucleon 
structure; its confirmation of the quark-parton model and its verification of the 
quark charge assignments being two notable early applications. In this thesis we 
describe and report the results of an absolute measurement of the dimuon cross 
sections from the interaction of an 800-GeV proton beam with hydrogen and deu- 
terium targets. As will be shown, these data provide important information about 
the distributions of antiquarks in the nucleon sea. 

1.1 Background 

The disciplines of high-energy and particle physics consist of the search for the 
ultimate constituents of matter and the study of the fundamental interactions 
between them. With the introduction of Quantum Chromodynamics in the mid- 
1970s, a consistent treatment of the weak, electromagnetic and strong interactions 
was finally realized. The gauge field theories of Quantum Electrodynamics (QED) 
and the Weinberg-Glashow-Salam model had long treated the electromagnetic 
and weak interactions felt by the leptons - the electron (e), muon (//), tau (t), 
and associated neutrinos (z/ e , u r ) - as the exchange of gauge field bosons 
(the photon (7) and the W ± and Z° bosons). Attempts to formulate similar 
theories for the strongly interacting hadrons, however, initially met with only 
limited success. The realization by Gell-Mann and (independently) Zweig in 1964 
that the hadron spectrum could be explained in terms of three fractionally charged 
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Table 1.1: Properties of the light quarks. Shown are the isospin (7) and its third 
component (Is), strangeness (S), baryon number (B) and charge (Q/e). 



flavor 


I 


h 


s 


B 


Q/e 


u 


1/2 


1/2 





1/3 


2/3 


d 


1/2 


-1/2 





1/3 


-1/3 


s 








-1 


1/3 


-1/3 



fermionic constituents, called quarks, set the stage for the development of QCD 
and its inclusion in what has become known as the Standard Model. 

Although Gell-Mann and Zweig's static quark model needed only three differ- 
ent quark flavors, experiment confirms the existence of six, with larger numbers 
being ruled out in the framework of the Standard Model by recent experiments. 
The six quarks are conventionally known as up (u), down (d), strange (s), charm 
(c), bottom (b) and top (t) and are found only in bound states called hadrons: 
mesons (7r ± , 7r°, J/ip, etc..) are quark- ant iquark (qq) bound states, while baryons 
(p, n, A ++ , etc..) are the bound states of three quarks (qqq). Table 1.1 shows 
the quantum numbers of the three quarks of the first three quarks. In order to 
avoid violating the Pauli principle in baryons such as the A ++ (uuu), the static 
quark model had to add an additional quantum number called color. Each quark 
carries one of three colors (red, green or blue). An object in which all three colors 
(or one color and its anti-color) are present is a color-neutral (color-singlet) state. 
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Color was an ad hoc addition to the static quark model. The role which it plays 
in QCD is more fundamental. QCD is a non-Abelian color gauge field theory, the 
three colors representing the fundamental basis of the SU(3) symmetry group. 
Color plays essentially the same role in the strong interaction as electric charge in 
electromagnetism, the principal difference being the number of "charge" states - 
two electric charge states (+ and — ) versus six color states (r, g, b, f, g, b). 2 The 
quarks interact through the exchange of QCD's gauge field bosons, the gluons 
(g). Unlike the photons in QED, which are electrically neutral, the gluons carry 
color charge. This leads to some important differences. A well-known property 
of QED is the running of the electromagnetic coupling constant a ~ j^f. A test 
charge polarizes the virtual e + e~ pairs surrounding it, leading to screening of the 
test charge at larger distances, resulting in an effective reduction in a. Probing 
the test charge at smaller distance scales, i.e. larger momentum transfers (q 2 ), 
results in an observed increase in the coupling. The same effect occurs in QCD, 
but the color-charged gluons are self interacting. This leads to an antiscreening of 
the color test charge, resulting in a decrease in the strong coupling constant (a s ) 
with increasing q 2 . This property of the strong interaction is known as asymptotic 
freedom. 

Owing to the mathematical complexity of the SU(3) ® SU(2) ®U(1) Standard 

Model, exact solutions are often not possible for physical interactions. We must 
2 Red, green, blue and their anti-colors. 



4 



therefore rely upon perturbation theory to obtain meaningful predictions from the 
model. This has never presented a problem in either the weak or electromagnetic 
interactions, where one can simply carry out the perturbative expansion in the 
coupling constant to whatever level of precision is desired. However, the strong 
interaction derives its name from the fact that it is a strong interaction. For 
momentum scales on the order of the proton mass, the size of the strong coupling 
constant {a s {M 2 ) ~ 1) renders fixed-order perturbation theory meaningless - all 
orders in the perturbative expansion are significant. Fortunately, the asymptotic 
freedom exhibited by QCD results in a strong coupling constant which decreases 
with increasing q 2 . At momentum-transfers typical of high-energy fixed-target 
experiments, a s becomes sufficiently small that it may be reliably used as an 
expansion parameter. Nevertheless, it remains sufficiently large that terms next- 
to- leading order in a s may remain important. 

An additional complication in hadronic interactions is that we cannot be cer- 
tain which of the interacting partons (the quark, antiquark and gluon constituents) 
are participating in the interaction. We can at best speak of the probability of 
finding a given parton carrying a fraction x (called Bjorken-x) of the interact- 
ing hadron's momentum, which we in general denote by f q /A( x ) f° r a parton of 
type q in a hadron of type A. Calculations of physical processes in the Stan- 
dard Model are thus predicated on prior knowledge of these parton distribution 
functions (PDF's). 
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Of particular importance to many experiments are the parton distributions 
in the nucleon (i.e. u(x) = f u /p(x), u(x) = fu/ P (x), etc.). Early parame- 
terizations of the PDF's relied on fits to the structure functions measured in 
deeply-inelastic lepton-hadron scattering experiments (DIS), which are primarily 
sensitive to the light quark distributions (u(x) and d(x)). More modern param- 
eterizations, such as those performed by the CTEQ, MRST and GRV collabo- 
rations [MRST, CTEQ, GRV] use several different physical processes to extract 
complementary information about the parton distributions. The lepton-charge 
asymmetry observed in W ± production provides additional information about the 
light quark distributions, while jet production and prompt photon measurements 
are used to constrain the gluon distributions. The Drell-Yan process for dilepton 
production, which involves quark- ant iquark annihilation, provides constraints on 
the light antiquark distributions (u(x) and d(x)) in the nucleon sea. The CTEQ 
5 set of parton distributions is shown in figure 1.1, along with the light antiquark 
distributions from the MRST 98 and GRV 98 distributions. 

1.2 The Drell-Yan Process 

Prior to the introduction of QCD, hadronic interactions were calculated with 
some success using the quark-parton model, whereby the interacting partons 
(quarks and antiquarks) were treated as being free of the strong interactions on 
the timescale of a hard interaction (the impulse approximation). Cross sections 
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X 



Figure 1.1: CTEQ5 parton distributions. Shown are the valence quark and light 
anti-quark distributions, and the gluons. The inset compares d(x) and u(x) from 
the CTEQ5 [CTEQ], MRST 98 [MRST] and GRV 98 [GRV] fits. 



7 



involving hadrons were then reduced to the comparatively simple problem of cal- 
culating the parton-level cross sections and folding in the probability to find the 
particular parton configuration. Figure 1.2 shows the Feynman diagram for the 
Drell-Yan process for lepton pair production in proton-nucleon collisions in the 
dimuon channel studied by E866. In the initial state, a quark (antiquark) carrying 
a fraction X\ of the beam's momentum annihilates an antiquark (quark) carry- 
ing a fraction x 2 of the target's momentum. This results in a time-like (q 2 > 0) 
intermediate photon, which decays into the final state fj, + fj,~ pair. 

1.2.1 Kinematics 

In the lab, we measure the invariant mass of the muon pair, and their total 
longitudinal momentum. The mass of the dimuon is related to the momentum 
fractions of the interacting partons by 



where s is the total four-momentum squared of the interacting hadrons. The 
longitudinal momentum of the dimuon, as a fraction of its maximum possible 
value (p™ ax y/s/2), is referred to as Feynman-x (x F ). It is related to the initial 
state momentum fractions by 



q 2 = M 2 = sx±x 2 




PL 



(1.2) 



Xp = 



.max 



— Xi- X 2 . 
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Thus, measurement of the invariant mass and of the muon pair allows us to 
determine the momentum fractions carried by the interacting partons. Solving 
for x\ and x 2 in equations 1.1 and 1.2 we find 

x 1>2 = ±^(x F ± ^x 2 f + At) (1.3) 

where we have taken r = X\X 2 = M 2 /s. 




Figure 1.2: The Drell-Yan process for dimuon production. 

In order to completely specify the final state of the interaction, we need four 
more variables. One of these is the momentum of the muon pair transverse to 
the beam, denoted by px- This is expected to be small in the Drell-Yan process: 
< Pt > ~ 0.3 GeV, non-zero only due to Fermi smearing of the momenta of 
the interacting partons. The larger values of pt which are observed in dilepton 
production require additional partons in the final state, which are provided by the 
higher-order QCD processes which we discuss later. 
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H+H~ CM. 





-y 




/ V 







pN plane 



Figure 1.3: Definition of the Collin-Soper frame. 

The three other variables which characterize the process are the polar decay 
angle 9d, and the azimuthal production and decay angles <fi p and <fid respectively. 
These variables are most naturally defined with respect to the qq annihilation 
axis. However, when is non-zero, this axis becomes difficult to measure. The 
angular variables are therefore measured with respect to the Collins-Soper frame 
[Col 77], which is shown schematically in figure 1.3. In this frame, the z-axis is 
taken to be parallel to the bisector of the angle between the interacting hadrons 
in the rest frame of the virtual photon. 

1.2.2 Cross Section 



The cross section for the Drell-Yan process is easily constructed in the quark- 
parton model. The assumption that the interacting quark and antiquark are free 
on the timescale of the hard interaction allows us to factor the cross section into 
a short distance qq annihilation term (<J q q), and a term expressing the probability 
of finding the particular partonic configuration. The qq annihilation cross section 
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is simply given by [Dre 70] 

14 7ra 2 2 A . 
a -"=3 3M^ (L4) 

which is just the cross section for the annihilation of two fermions in QED, mul- 
tiplied by a factor of 1/3 required by color conservation and a factor of e 2 q to 
account for the fractional charge of the quarks. Multiplying by the probability of 
a quark (antiquark) in the beam annihilating an antiquark (quark) in the target 
and summing over all parton species we have an expression for the Drell-Yan cross 
section differential in x\ and x 2 

J^T = \jM^ S e\[f q/A (x{)f q / B (x 2 ) + f q , A (x{)f q i B (x 2 )]. (1.5) 

q£u,d,s,... 

Expressed in terms of quantities measured in the lab we obtain 

M2 JM^t F = sW^fk £ e * lf*/*Mh*M + f*/AMM*2)] (i.6) 

where equation 1.3 is used to evaluate X\ and x 2 . Alternatively, we can express 
the cross section in terms of the dimensionless variables r and y = | In ^ (the 
rapidity) 

S ~c^y = \~^T ^ ^ (/«m( Xi )/9/b(s2) + fq/A(x 1 )f q/B (x 2 )) . (1.7) 
This is known as the scaling form of the cross section. 

1.2.3 Predictions 

The Drell-Yan model makes some important predictions about the dilepton cross 
sections. The left-hand side of equation 1.7 is a function of the variables r, y 
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and s, while the right-hand side is only a function of r and y. Thus, the model 
predicts that the cross sections multiplied by s, and measured at different beam 
energies, but the same kinematics, should be the same - the cross sections should 
scale with ^, and this is generally observed [Yoh 78, Con 89, Ant 82]. 

Another important prediction is the shape of the angular distributions. The 
angular distribution of the fx + in the pair rest frame can be expressed as [Cle 78] 

— — ; — oc 1 + A cos 2 9d + fx sin 29 d cos0d + v 12 sin 2 9d cos 26 d (1.8) 

dcpd d cos 9 d 

where A, fx and v are functions of the other kinematic variables. For the Drell- 
Yan process, the assumption of massless quarks implies that the virtual photon is 
transversely polarized (A = 1, fx — v — 0), resulting in the prediction 

-ocl + cos 2 ^ 1.9 

d cos 9d 

when 9d is measured with respect to the qq annihilation axis. Several experiments 
[Con 89, McG 99, Cha 99, Bro 01] have confirmed that the dileptons are produced 
transversely polarized. 

Despite these successes in describing experimental data, the Drell-Yan model 
also had some initial failures. Because there are no particles in the final state off 
which to recoil, any transverse momentum in the dilepton must originate with the 
initial state. The mechanism for this is most likely Fermi smearing of the parton 
distributions, and the expected size of this intrinsic kx is small (~ 300 MeV). The 
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observed (s-dependent) values of (jpr) ~ 1.2 GeV cannot be reproduced by this 
mechanism in the Drell-Yan process. 

A second problem came about when absolute measurements of the dilepton 
cross sections were compared to Drell-Yan calculations based on the early PDF 
fits to DIS structure functions. Table 1.2 shows the iCfactors, defined as the ratio 
of measured to calculated cross sections (K = o" m cas./o"DY), from several different 
experiments using different beams and targets. Typically, experiments measured 
a cross section which was a (nearly constant) factor of two larger than the Drell- 
Yan model and DIS structure function measurements implied. Along with the 
discrepancy in the pr distribution, this was taken as a sign that QCD would play 
an important role in dilepton production. 

1.3 QCD Modifications 

With the introduction of QCD, the basic physics behind the Drell-Yan process was 
confirmed theoretically. That quark-antiquark annihilation gives rise to dilepton 
pairs is a consequence of the electromagnetic interaction. QCD provides substan- 
tial, albeit well-known, corrections to the leading-order (LO) electromagnetic con- 
tribution described by Drell and Yan. The diagrams which contribute at next-to- 
leading order (NLO or 0(a s )) are shown in figure 1.4. The gluon-bremsstrahlung 
and vertex-correction diagrams represent corrections to the basic qq annihilation 
picture of the Drell-Yan process, increasing the magnitude of the theoretical cross 
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Table 1.2: Experimental K-factors. 



Experiment 


Interaction 


Beam Momentum 


K — <7 meas ./<7DY 


E288 


[Kap 78] 


p Pt 


300/400 GeV 


~ 1.7 


WA39 


[Cor 80] 


7T ± W 


39.5 GeV 


~ 2.5 


E439 


[Smi 81] 


p w 


400 GeV 


1.6 ±0.3 






(p - p)pt 


150 GeV 


2.3 ±0.4 






p Pt 


400 GeV 


3.1 ±0.5 ±0.3 


NA3 


[Bad 83] 


7T ± Pt 


200 GeV 


2.3 ±0.5 






71- Pt 


150 GeV 


2.49 ±0.37 






71- Pt 


280 GeV 


2.22 ±0.33 


NA10 


[Bet 85] 


71- W 


194 GeV 


~ 2.77 ±0.12 


E326 


[Gre 85] 


71- W 


225 GeV 


2.70 ±0.08 ±0.40 


E537 


[Ana 88] 


p w 


125 GeV 


2.45 ±0.12 ±0.20 


E615 


[Con 89] 


71- W 


252 GeV 


1.78 ±0.06 



section and providing a mechanism to reach larger values of p?- Resummation 
techniques have been developed, which enable the calculation of such corrections 
to all orders in a s . Taking K equal to the resummed cross section divided by 
the LO contribution gives K 1.8 for a s evaluated at q 2 typical of fixed-target 
experiments. This already accounts for much of the experimental K-factor. 
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The contributions described above are simply modifications to the basic qq 
picture of the Drell-Yan process - the main interaction is still the electromagnetic 
annihilation of the interacting partons, with the strong interaction only modifying 
the initial states of the qq pair. The gluon Compton-scattering diagram, however, 
is fundamentally different. Here, the incoming quark (antiquark) scatters off a 
gluon through the strong interaction, emitting the virtual photon which decays 
in the final state of the interaction. Theoretical studies of this process by Berger 
et al. [Ber 98] have shown that for transverse momenta large compared to the 
mass of the pair (p T > M/2), the gluon Compton scattering diagrams dominate 
the cross section (accounting for pa 80% of the NLO calculation). Therefore, with 
increasing pt, the dilepton cross section becomes less sensitive to the antiquark 
distributions, and more sensitive to the gluons. 

1.4 Motivation 

The primary goal of Fermilab Experiment 866 (FNAL E866/NuSea) was to mea- 
sure the ratio of continuum dimuon cross sections in proton-deuteron (a p d) in- 
teractions to those in proton-proton (cr pp ) interactions. These measurements are 
sensitive to the x-dependence of the ratio d/u in the nucleon, and were moti- 
vated by several experiments [Ama 91, Bal 94, Ada 95] which demonstrated a 
large and previously unanticipated flavor asymmetry d/u > 1. The E866 re- 
sults [Haw 98, Tow 01], which were based on a large sample of dimuon events 

15 



Gluon Bremsstrahlung (qqg) 




Gluon Compton Scattering (qg and qg) 




Vertex correction (qq) 




Figure 1.4: Next-to- leading order contributions to the dilepton cross section. 
The final state muon pair is omitted. 
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(« 360, 000) induced by an 800-GeV proton beam (y/s = 38.8 GeV) interact- 
ing with hydrogen and deuterium targets, provided the first determination of the 
strong x-dependence of this flavor asymmetry. 

Absolute measurements of the hydrogen and deuterium cross sections, the 
topic of this thesis, provide additional information on the magnitude and absolute 
shape of the light antiquarks. This can be easily shown by writing the leading- 
order cross section expression in equation 1.6 for various targets. Evaluating the 
sum for the hydrogen cross section, and taking the limit of large Xp appropriate 
for much of the E866 data 3 we have 

{ M2 i£k) n>> „ {r {M + l d(x ^) ■ (hl0) 

Owing to the larger charge-squared factor in front of the uu term and the size of the 
up distribution compared to that of the down distribution in the proton, about 
80% of the proton cross section at leading order is due to the uu annihilation. 
Thus, to the extent that the valence distributions are known (DIS measurements 
constrain these well in our kinematic range), the proton cross section provides a 
direct measurement of the u(x) distribution. 

The deuterium cross section can be similarly written. Nuclear effects in deu- 
terium are small, and confined to small x 2 [Aid 90] . This allows us to express the 

deuterium cross section as the sum a p d = cr pp + <j pn . If we assume that charge 

3 In the large xf limit, x\ is sufficiently large compared to x 2 that we can neglect terms in 
the cross section with the antiquark in the beam and the quark in the target. 
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symmetry holds, then the parton distributions in the proton and neutron are re- 
lated through u(x) = d n (x), u(x) = d n (x), etc.... This yields an expression for the 
deuterium cross section (per nucleon) 

iuM + dM }[3( Xl ) + (1.U) 

where we have again taken the large-x^ limit. Thus, the deuterium cross section 
is directly proportional to the sum of the light antiquarks d(x) + u(x). 

In this thesis we describe and report results on the absolute measurement of 
continuum dimuon production in 800-GeV pp and pd interactions. The doubly- 
differential (M 3 d 2 a/dMdx F ) and triply-differential (Ed 3 a/dp 3 ) cross sections have 
been measured over a large range in the mass (4 < M < 16.85 GeV), x F (—0.05 < 
xf < 0.8) and pr (0 < pr < 7 GeV) of the muon pair. These data represent 
the first such measurement of the pp cross sections over an extended range of 
kinematics, while the pd cross sections have been measured over a wider kinematic 
range with greater statistical precision than has previously been achieved. 

Traditionally, the Drell-Yan process is discussed in the context of probing the 
light antiquark distributions in the proton. And indeed, that was the primary 
motivation of this experiment. But it is clear from equations 1.10 and 1.11 that 
dimuon production also provides sensitivity to the valence distributions. While 
the kinematic range covered by previous dimuon measurements [Mor 91] has been 
limited to regions where the valence quarks are well constrained by DIS measure- 
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merits, the E866 data cover a much wider kinematic range. Our data probe the 
valence out to x ~ 0.8, which provides an important measurement of the large- 
x behavior of the valence quarks complementary to existing DIS measurements 
[Zha 01]. 
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2 APPARATUS 



Fermilab Experiment 866 used a modified version of the dimuon spectrometer 
[Cri 86] located in the Meson East experimental area at the Fermi National Ac- 
celerator Laboratory. Previous experiments which used this spectrometer were 
E605, E772 and E789. Figure 2.1 shows a drawing of the E866 spectrometer, 
which was designed to detect oppositely-charged muon pairs while minimizing 
sensitivity to everything else. The coordinate system of the spectrometer was 
defined in the following way: the z-axis of the spectrometer was parallel to the 
nominal direction of the beam; the x-axis laid parallel to the floor of the experi- 
mental hall, pointing to the left facing station-1 from the targets; the y-axis then 
formed the right-handed Cartesian coordinate system by pointing up, perpendic- 
ular to the floor of the experimental hall. The origin of the coordinate system was 
chosen as the center of the upstream face of the SM12 magnet. 




Figure 2.1: FNAL E866/NuSea Spectrometer. 



20 



2.1 Beam line 

An 800-GeV proton beam was extracted from the Fermilab Tevatron, split at the 
switchyard and sent to several fixed target experiments which ran concurrently 
during the 1996 fixed target run. The beam was transported down the east beam 
line of the Meson experimental area and arrived at the target in 20-second-long 
spills, with new spills arriving once every minute. The accelerator's frequency 
was 53 MHz, which meant that during the spill the protons arrived grouped into 
"buckets" separated by 19 ns. The typical intensities which were used for E866 
were between 1 x 10 11 and 2 x 10 12 protons per spill, depending on trigger rates 
and radiation-safety requirements for a particular spectrometer configuration. 

The position, size and intensity of the beam was constantly monitored at 
various points along the beam line by several detectors. The beam's position and 
size were monitored by segmented wire ion chambers (SWICs). The last SWIC in 
the beam line was located approximately 188 cm upstream of the target. It had 
a 2 mm horizontal wire spacing and a 0.5 mm vertical wire spacing. During the 
experiment, the beam spot on the target was approximately 6 mm wide by 1 mm 
high. 

Three different types of detectors monitored the intensity of the beam. The 
primary detector was a secondary-emission monitor (SEM), located at about 100 
m upstream of the targets. In addition to the SEM counter, an ion chamber (IC3) 
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and a quarter-wave RF cavity were used. These additional detectors allowed us 
to study the linearity and offsets of the SEM. 

The luminosities of the targets were also monitored during the experiment 
with a pair of four element scintillator telescopes located about 85° from the beam 
direction in the horizontal plane. The data from these detectors (referred to as 
AMON and WMON) were used to determine beam duty factor, data acquisition 
live-time and which target was in the beam during the spill. 

2.2 Targets 

Three identical target flasks, one of which is shown in figure 2.2, were mounted 
on a movable table located at z rs —400 cm. They consisted of a 17-inch long 
cylindrical shell with a 3-inch inside diameter. Hemispherical end caps with a 
3-inch diameter sealed either end, creating a 20-inch (50.8-cm) long target along 
the beam axis. 1 Facing the targets from upstream, the leftmost target was filled 
with liquid hydrogen and the rightmost target was filled with liquid deuterium. 
The center target was empty and held at vacuum to measure the background rates 
due to the beam interacting with the SWIC, the material in the target flasks, and 
the air and other materials downstream of the target. Immediately downstream 

of the target, filling the aperture of the SMO magnet, was a bag filled with helium 

1 It was confirmed by visual inspection after the experiment that the z-axis of each target 
was aligned properly with the beam. 
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Figure 2.2: One of the three identical target flasks used in E866. 

to reduce background interactions in the vicinity of the targets. 

The targets were cycled periodically to minimize the systematic uncertainties 
between the hydrogen and deuterium measurements. The programed sequence 
was 5 spills on the deuterium target followed by 1 spill on the empty target, 5 
spills on the hydrogen target, and one more spill on the empty target. The target 
control computer counted only "good" spills - spills which the DAQ determined 
had a sufficiently large number of protons. Upon receipt of the appropriate number 
of good spill signals, which were sent after the end of the spill, the target-control 

23 



computer initiated the move to the next target in the sequence. The move was 
completed before the beginning of the next spill. 

A set of four switches on the target table were engaged by a stationary roller 
mounted in front of the table. The positions of these switches corresponded to the 
beam striking the center of each of the targets, plus a fourth position where the 
targets were clear of the beam. When engaged, each switch completed a circuit 
which was read out by the DAQ, recording which target was in the beam. The 
switches were also tied into the beam interlock system, which required that one 
of the four switches be engaged for beam to be allowed to enter the target area. 
This ensured that the beam would not be able to hit the sides of the target flasks 
and cause a radiation hazard. 

2.3 Spectrometer Magnets 

The E866 spectrometer used three dipole magnets whose magnetic fields were 
oriented in the x direction (bending charged particles in the y direction). The two 
magnets closest to the targets were SMO and SM12. These magnets were used to 
optimize the acceptance of the spectrometer for different ranges of dimuon masses. 
The SM12 magnet was used to focus dimuon pairs on the downstream detector 
elements. Access to higher-mass events, which tend to have larger opening angles, 
was possible by using larger magnetic fields in SM12. This would ensure that 
such events were kept inside the aperture of the SM12 magnet and the detector 
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stations. 

The SMO magnet was used to optimize the spectrometer for lower mass dimuon 
events in the low-mass data set. In this data set, the SMO magnet was used 
to defocus (increase the opening angle of) the pairs before they entered 

SM12. This increased the probability that the muons would miss the dump. 
More importantly, this increased the probability that the events would miss the 
dump cut used to eliminate beam-like muons, which were indistinguishable from 
the high flux of muons produced in the dump. 

The SM3 magnet was located between the station-1 and station-2 detectors. It 
was used in the experiment to analyze the momentum of the muons by measuring 
the deflection it produced in the tracks. It was operated using a single current of 
4230 A throughout the experiment. 

2.4 Beam Dump and Hadronic Absorbing Wall 

In order to prevent damage to the downstream detectors from the beam the target 
was followed by a large, water-cooled beam dump whose upstream face was located 
at z ~ 173 cm. The beam dump consisted of 255 cm of copper, which spanned 
-15.25 cm < y < 15.25 cm at its thickest point. At its widest, the dump blocked 
the entire x-aperture of the magnet. The beam dump was nominally centered on 
the z-axis, but had settled over the years and developed a small slope relative to 
SM12 such that it was approximately 0.2" low at z — 68" and 0.4" low at z — 530". 
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A large hadronic absorbing wall was located immediately downstream of the 
beam dump. The purpose of this wall was to fully attenuate the intense flux of 
hadrons produced in the target and the dump. The absorbing wall consisted of 
61 cm of copper 2 in the most upstream section, followed by a 205.8 cm section of 
carbon, a 68.6 cm long section containing both carbon and polyethylene, and a 183 
cm section of polyethylene. The configuration of the beam dump and absorbing 
wall is shown in figure 2.3. 

2.5 Detector Stations 

Four detector stations were responsible for triggering on dimuon events and track- 
ing the muons which traversed the spectrometer. The first three detector sta- 
tions consisted of hodoscope- and drift chamber-planes, while station-4 consisted 
of hodoscope- and proportional-tube planes. The fast response time of the ho- 
doscopes was utilized to trigger on events with the expected dimuon signature, 
at which point the slower responding drift chambers could be read out. The 
proportional-tubes and hodoscopes at station-4 were used to provide some posi- 
tion information to the trigger, and discriminate against any hadrons which made 

it through the absorbing wall. 

2 There was a 7.6 cm gap between the copper and carbon sections of the absorbing wall in 
the region behind the dump. 
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Figure 2.3: Side and elevation views of the absorbing wall and beam dump. 
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2.5.1 Hodoscope Planes 

Hodoscope planes were located at each of the four detector stations and were used 
to trigger on dimuon events. Each hodoscope plane was split into two half-planes 
of parallel scintillator paddles attached to photomultiplier tubes by plexiglass light 
guides. Stations 1, 3 and 4 each had two hodoscope planes, with their paddles 
oriented in the X and Y directions respectively. Station-2 had only one hodoscope 
plane, with its paddles oriented parallel to the floor. A vertical gap between the Y 
half-planes and the central X hodoscopes was inserted in order to avoid triggering 
on the high rate of beam-like muons from the decay of pions which were produced 
copiously in the dump. The design specifications of the hodoscope planes can be 
found in table 2.1. A more precise alignment of the hodoscopes was achieved by 
examining the distributions of positions of muon tracks at each hodoscope plane 
when a given hodoscope in that plane fired. 

2.5.2 Drift Chambers 

Each of the first three detector stations contained six drift-chamber planes, ar- 
ranged in three pairs of planes with parallel wire orientations (referred to as 
"views"). Wires were oriented horizontally in the "Y" view, and at an angle 
of +14° and —14° with respect to the X axis in the "V" and "U" views, respec- 
tively. The wires in the second plane in each pair were offset by half the cell size 
of the drift chamber. The plane in each pair closest to the target was denoted as 
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Table 2.1: Specifications of the hodoscope planes. Distances measured in inches. 



detector 


z position 


# of counters 


counter width 


gap 


Ax 


Ay 


Yl 


769.78 


16 


2.5 


0.38 


47.50 


40.75 


XI 


770.72 


12 


4.0 


0.47 


47.53 


40.78 


Y2 


1114.94 


16 


3.0 


0.66 


64.625 


48.625 


X3 


1822.00 


12 


8.68 


1.0 


105.18 


92.00 


Y3 


1832.00 


13 


7.5 


0.0 


104.00 


92.00 


Y4 


2035.50 


14 


8.0 


0.0 


116.00 


100.00 


X4 


2131.12 


16 


7.125 


0.0 


126.00 


114.00 



the "unprimed" plane, while the plane in the pair furthest from the target was 
denoted as the "primed" plane. These planes provide information on the X and Y 
intercept of the muon tracks at the detector station, with redundant information 
about the Y position. 

The drift chambers were operated with a gas mixture of 49.7% argon, 49.6% 
ethane and 0.7% ethanol, mixed by volume at a constant temperature of 25°F. 
The anode wires at Station- 1 were gold-plated tungsten wire. Stations 2 and 3 
used silver-coated beryllium-copper wires as anodes. All of the anodes were 25 /xm 
in diameter. The cathode wires were all 62.5-/zm silver-coated beryllium-copper 
wire. The drift chambers were operated at voltages between 1700 and 2200 volts. 
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Typical drift velocities were w 50 /zm/ns. Drift chamber specifications can be 
found in table 2.2. 

The signals from the drift chambers were read out by a fast amplifier and 
discriminator system. Single-hit time-to-digital converters (TDCs), which only 
record the first hit on the wire during an event, were used to measure the drift 
time. The combination of good hits together with their associated drift times in 
all three views result in a "triplet" hit in each station. The bank of the triplets 
was saved to provide information used later in reconstructing the muon tracks in 
the analysis. 

2.5.3 Proportional Tubes 

The detectors at station-4 provided both trigger-information and muon discrim- 
ination capabilities. Located downstream of the electromagnetic and hadronic 
calorimeters, 3 station-4 consisted of two hodoscope planes (X4 and Y4) and three 
proportional-tube planes (PTY1, PTX and PTY2). Each of the proportional- 
tube planes had two layers of 1 x 1-inch cells. Adjacent layers were offset by half a 
cell to cover the dead region between adjacent cells. The proportional-tubes were 
operated using the same gas mixture used in the drift chambers. To reduce the 

probability of any hadrons making it through the calorimeters, a 3 foot wall of zinc 

3 The calorimeters were left over from a previous experiment. For E866, they were un- 
necessary and were only used to eliminate any hadrons which made it through the absorbing 
wall. 
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Table 2.2: Specifications of the drift chamber planes. Distances measured in 
inches. 



1 J_ J_ 

detector 


/-position 


# of wires 


1 1 

cell size 


J_ / "\7" A T \ 

aperture(XxY) 


operating voltage 


VI 


724.69 


200 


0.25 


48x40 


+ 1700 


VI' 


724.94 


200 


0.25 


48x40 


+ 1700 


Yl 


740.81 


160 


0.25 


48x40 


+ 1700 


Yl' 


741.06 


160 


0.25 


48x40 


+ 1700 


Ul 


755.48 


200 


0.25 


48x40 


+ 1700 


Ul' 


755.73 


200 


0.25 


48x40 


+ 1700 


V2 


1083.40 


160 


0.388 


66x51.2 


-2000 


V2' 


1085.52 


160 


0.388 


66x51.2 


-2000 


Y2 


1093.21 


128 


0.40 


66x51.2 


-2000 


Y2' 


1095.33 


128 


0.40 


66x51.2 


-2000 


U2 


1103.25 


160 


0.388 


66x51.2 


-1950 


U2' 


1105.37 


160 


0.388 


66x51.2 


-1975 


V3 


1790.09 


144 


0.796 


106x95.5 


-2200 


V3' 


1792.84 


144 


0.796 


106x95.5 


-2150 


Y3 


1800.20 


112 


0.82 


106x91.8 


-2200 


Y3' 


1802.95 


112 


0.82 


106x91.8 


-2200 


U3 


1810.24 


144 


0.796 


106x95.5 


-2200 


U3' 


1812.99 


144 


0.796 


106x95.5 


-2200 
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Table 2.3: Specifications of the proportinal tube planes. Distances measured in 
inches. 



detector 


Z-position 


No. of wires 


cell size 


aperture XxY 


PTY1 


2041.75 


120 


1.0 


117x120 


PTX 


2135.875 


135 


1.0 


135.4x121.5 


PTY2 


2200.75 


143 


1.0 


141.5x143 



and a 4 inch wall of lead were placed between the calorimeters and the station- 
4 detectors. Furthermore, 3 feet of concrete was placed between PTY1 and X4 
and between PTX and PTY2. This provided a total of 16.6 hadronic-interaction 
lengths upstream of the Y4 plane. The only detectable particles which could reach 
the station-4 detectors were muons. Specifications of the proportional-tubes may 
be found in table 2.3. 

2.6 Trigger 

Despite the fact that the absorbing wall all but eliminated the rate of hadrons 
produced in the dump, the rates of single and dimuon events from both the target 
and dump were too large for the drift chambers to be read out on every muon 
which was detected. Therefore the information from the hodoscopes was used to 
trigger on events which had the expected signature of a dimuon event originating 
in the target. The trigger was also configured to take prescaled samples of events 

32 



which would be used later to reconstruct the combinatoric background. 

The trigger system used in E866 was upgraded from that used for previous 
experiments [Gag 98]. The purpose of the upgrade was to increase the flexibility 
of the trigger, and to improve the acceptance of high p T dimuon events. The main 
elements of the trigger system were the Trigger Matrix Modules and the Track 
Correlators. Figure 2.4 shows a diagram of the trigger electronics for the left hand 
side of the spectrometer. 
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Figure 2.4: Block diagram of the E866 trigger. Only the inputs and electronics for the left side of the 
are shown. 



The Trigger Matrix Modules were responsible for determining whether a track 
was likely to be from a dimuon event which originated at the target. The outputs 
of the Y hodoscopes at stations 1, 2 and 4 were input into the Trigger Matrix 
Modules. The hit patterns registered by the hodoscopes were tested against a 
preselected set of hit patterns which corresponded to a muon from a dimuon event 
in the target traversing the spectrometer. These hit patterns where determined 
from Monte Carlo studies prior to the experiment. Some matrix elements had 
considerable rates of muons passing through them which actually originated in 
the dump. In order to take larger rates of target events, these matrix elements 
were turned off. 

Four sets of Trigger Matrix Modules were used in the trigger system. These 
were referred to as matrix-up-left (MUL), matrix-down- left (MDL), matrix-up- 
right (MUR) and matrix-down-right (MDR). They corresponded to single muons 
on the left- or right-hand side of the spectrometer which went above or below the 
dump. These, along with the outputs of the Terminator/OR's and the Station-4 
XY (S4XY) Track Correlators, were the inputs to the Track Correlators. The 
job of the Track Correlators was to filter the inputs and determine if a trigger 
should be fired and the event written to tape. There were three Track Correlators 
capable of triggering an event: Physics A, Physics B and Diagnostic. 

Physics A (PhysA) was the main physics Track Correlator. Like all of the 
Track Correlators, it could handle four sets of trigger conditions. Table 2.4 shows 
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the triggers that were used in each of the data sets in this analysis. The trig- 
gers in each of the data sets performed the same functions. The PhysAl trigger 
selected oppositely-charged muons which traversed opposite sides of the spectrom- 
eter, while the PhysA3 and PhysA4 triggers selected oppositely-charged muons 
which traversed the same side of the spectrometer. The PhysA2 trigger was the 
like-sign trigger. It selected same-sign muon pairs on opposite sides of the spec- 
trometer. 

In addition to the requirements on the charge and geometry of the muons 
applied by the matrix, the Track Correlators used in the low- and intermediate- 
mass data sets added constraints on the X and Y positions of the muons at 
station-4 (the S4xxx entries in the table.) These constraints were designed to 
veto events which populated areas where the dump dominated the event rate. 

The Physics B (PhysB) track correlator was used to sample events which were 
used to monitor the efficiencies and systematics of the detector. An example of one 
of the many B triggers used in the experiment appears in table 2.5. The PhysBl 
trigger was designed to read out any event with muons traveling along paths on 
opposite sides of the spectrometer, regardless of their relative sign. The PhysB2 
trigger was often configured to trigger on events with a single muon traveling 
along one of the predetermined matrix roads that a muon in a dimuon event 
would follow. These data were used to extract the combinatoric background, as 
will be discussed later. 
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Table 2.4: Trigger configuration (PhysA). Configurations shown for the low-, 
intermediate- and high-mass data. The symbols "*", "+" and "!" correspond to 
logical AND, OR and NOT, respectively. 



Mass Setting 


Number 


Prescale 


Requirements 


low 


PhysAf 


f 


(MUL * MDR) + (MUR * MDL) 




PhysA2 


4 


(MUL * MUR) + (MDL * MDR) 




PhysA3 


f 


(MUL * MDL) * (S4UL2 * S4DL2) 




PhysA4 


f 


(MUR * MDR) * (S4UR2 * S4DR2) 


int 


PhysAf 


f 


( (MUL * MDR) * ( !S4DLf + !S4URf ) ) + 








( (MUR * MDL) * ( !S4DRf + !S4ULf ) ) 




PhysA2 


2,1 


( (MUL * MUR) * ( !S4DLf + !S4DRf ) ) + 








( (MDL * MDR) * ( !S4ULf + !S4URf ) ) 




PhysA3 


f 


( (MUL * MDL) * 








( !S4ULf + !S4DLf ) * ( S4UL2 * S4DL2 ) ) 




PhysA4 


f 


( (MUR * MDR) * 








( !S4URf + !S4DRf ) * ( S4UR2 * S4DR2 ) ) 


high 


PhysAf 


f 


(MUL * MDR) + (MUR * MDL) 




PhysA2 


4 


(MUL * MUR) + (MDL * MDR) 




PhysA3 


f 


(MUL * MDL) 




PhysA4 


f 


(MUR * MDR) 
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Table 2.5: Trigger configuration (PhysB). An example of a PhysB trigger used in 
the high-mass data. The symbols "*", "+" and "!" correspond to logical AND, 
OR and NOT, respectively. 



Mass Setting 


Number 


Prescale 


Requirements 


high 


PhysB 1 


800 


X134L * X134R 




PhysB2 


1000 


( MUL + MDL + MUR + MDR ) 




PhysB3 








PhysB4 







2.7 Data Acquisition System 

The E866 data acquisition system (DAQ) was essentially an upgraded version 
of the system used in E789. The DAQ had three areas of responsibility: event 
readout, data archiving and online analysis. The backbone of the event readout 
was the Nevis Transport System. Upon receipt of a signal from one of the Track 
Correlators, a busy signal was raised which inhibited further triggers from being 
accepted. Simultaneously, the first word of the event was inserted onto the trans- 
port bus. Signals were sent to the TDC readouts on the drift chambers and the 
coincidence registers (CR's) on the hodoscopes. When signaled, the TDC's and 
CR's began digitizing their signals for insertion onto the bus. Each hodoscope and 
proportional-tube hit resulted in the insertion of an identifying word on the bus. 
The TDC's would also start a timer which would be terminated by the amplified 
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signal from the drift chamber. This provided a measurement of the drift time in 
the chamber, which was inserted onto the transport bus along with the number 
of the struck wire. 

At this point, the data stream on the transport bus was fed into the Versa 
Module Eurocard (VME) based archiving system. During the 20s spill cycle, data 
were first transported to a pair of high-speed memory boards. Once one of the 
boards was full, it was drained into a large memory buffer while the other board 
was filled. Upon completion of the spill, the data were then formatted and sent 
down the VME pipeline to the tape archiving system. Here the data were written 
out to 8mm Exabyte tapes. Additionally, a fraction of the event data and all of 
the spill information (number of triggers fired, SEM counts, target position, etc..) 
were sent to the online monitoring systems. 

The online monitoring system consisted of a database system interfaced to 
several graphical tools, and the E866 online analysis code. The database sys- 
tem allowed us to monitor the status of various components of the beamline and 
spectrometer. Graphical displays of the luminosity, magnet voltages, livetimes, 
etc., gave an overall indication of the health of the spectrometer. The online 
analysis code analyzed the data that were sent to it via the VME pipeline. The 
histograms generated by the analysis code were accessible in real time by making 
use of PAW (Physics Analysis Workstation) [PAW] global sections. This allowed 
us to monitor the detector planes, watching for and correcting any inefficiencies 
which might develop. oq 



3 EVENT RECONSTRUCTION 

By the time the experiment ended, approximately 250 GB of data had been 
recorded to tape. The data on these tapes represented the hits on each cham- 
ber and hodoscope plane for each event, the general characteristics of the beam 
in any given spill and the overall state of the spectrometer. The reconstruction 
of the kinematics of the events from these raw data required several steps. After 
configuring the analysis for the specific data set, the tracks were reconstructed 
from the chamber hits. Once the tracks had been found, their paths through the 
SMO and SM12 magnets had to be reconstructed while accounting for energy loss 
and multiple scattering in the absorbing wall and beam dump. These procedures 
are outlined in greater detail below. 

3.1 Data Sets 

During the nine months starting in September 1996 in which FNAL E866/NuSea 
took data, over 360,000 continuum dimuon events were written to tape. These 
events were grouped into several sequentially numbered data sets which were dif- 
ferentiated by the polarities and currents used to energize the three spectrometer 
magnets. Each data set fell into one of three categories, defined by the mass 
range for which the spectrometer was optimized. These categories were referred 
to as the low-, intermediate- and high- mass data, corresponding to spectrometer 
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settings optimized for masses near the J/ip, between the J/ip and T masses, and 
near the T mass, respectively. Table 3.1 lists the data sets used in this analysis 
and the mass setting to which they belong. It also tabulates the currents used to 
power the spectrometer magnets, and whether the targets were filled with either 
the first or second sample of cryogenic liquids. 

Table 3.1: Definition of the data sets. The SM3 magnet was operated at a single 
current of 4230 A, with the same polarity as the SM12 current. 



mass setting 


data set 


SMO current 


SM12 current 


target fill 


low 


5 


-2100 A 


2800 A 


first 




10 


2100 A 


-2800 A 


second 


intermediate 


9 


A 


2800 A 


second 


high 


7 


A 


4000 A 


first 




8 


A 


4000 A 


second 




11 


A 


-4000 A 


second 



3.2 First Pass 

Since only about ~ 1% of the events recorded to tape were due to a continuum 
dimuon event originating from the target, multiple analyses of the entire data sam- 
ple were not feasible. It was therefore necessary to reduce this data sample to more 
manageable levels. A first pass analysis was performed which filtered out many 



41 



of the bad (non-dimuon) events in the data sample, while preserving almost all of 
the events we were interested in. The computing resources for this analysis were 
provided by the Fermilab Computing Division. A number of IBM workstations 
were linked together into four parallel computing clusters called "farms", which 
were able to process the complete data sample in approximately two months. 

3.2.1 Track Reconstruction 

The analysis began by reading in the configuration of the spectrometer from var- 
ious files, some specified by the user and some loaded automatically based on the 
data set being analyzed. Once configured, the analysis code examined the hits 
in the drift chambers at stations 2 and 3 to find a set of candidate tracks. If 
four of the six planes in a station registered a hit, the position of that hit was 
considered a possible track reconstruction point. Points which were inconsistent 
with the trigger which fired were dropped. The remaining reconstruction points 
were iteratively combined to form track segments, so long as those segments were 
consistent with a single particle forming both points. 

The track segments were then compared with hits in station 1 to reduce further 
the number of track segments. Each segment was first extended to the SM3 bend 
plane, located between stations 1 and 2. The charge and momentum of the particle 
was not yet known, so a window-of-interest at station 1 which pointed back to the 
target in 
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the X direction was searched. If none of the hits at station 1 met the search 
criteria, the track segment was discarded. 

Once the track segment between stations 2 and 3 was connected with the 
station 1 reconstruction point (or points), the charge and momentum of the track 
were determined from the direction and size of the deflection of the track at SM3. 
The track segment was then extended back to station 4 to verify that it was 
indeed a muon. The concrete and lead shielding present at station 4 absorbed any 
hadrons and electrons which managed to penetrate the absorbing wall, but also 
caused any muons to multiple scatter. Since at this point the analysis knew the 
momentum of the track, a window consistent with the expected multiple scattering 
could be searched for hits in the five station 4 detectors. If the track fired less 
than three of the station 4 detectors, it was rejected on the grounds that it was 
not likely to be a muon. 

3.2.2 Event Kinematics 

At this point in the analysis, the charge, momenta and positions at the various 
detectors of the tracks in the event were known. It had also been determined 
at this point that each of the tracks was due to a muon. The tracks were then 
propagated back through the SM12 and SMO magnetic fields, accounting for the 
energy lost in the hadronic absorbing wall and (if applicable) the dump. The 
paths of the tracks were reconstructed inside the magnetic field in 46-cm steps, 
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using field maps which had been measured prior to the experiment. A correction 1 
was also made to compensate for the energy lost by the muon where the track 
passed through either the absorbing wall or the beam dump. 

In addition to suffering energy loss inside the absorbing wall and dump, the 
muons also underwent multiple Coulomb scattering (MCS). The correction for 
this began by reconstructing the track back to an effective scattering bend plane 
(ZSCPLN) and recording its position and momentum there. The track was then 
completely analyzed, noting the deviations in X and Y from the nominal target 
center. A correction to the momentum vector at the ZSCPLN was calculated 
from these deviations, and the reconstruction of the track repeated from that 
point. This procedure was iterated until the track converged on the nominal 
center of the target. 

The production of pions and other hadrons in the target contributed to a 
large rate of single muons in the spectrometer as they decayed between the target 
and dump. A significant fraction of these muons followed a beam-like trajectory 
through the dump and absorbing wall. Because there was no way to differentiate 
between one of these single muons and a muon which was a part of a valid dimuon 
event, these muons could cause ambiguities in the analysis that could not be re- 
solved. To prevent this from happening, muons which passed through a 10.2 cm 

1 The most probable energy lost in a given section of absorbing wall or dump was parameter- 
ized versus the momentum of the track. This energy loss parameterization was used to estimate 
the energy lost by the muon. 
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horizontal band centered on the nominal beam position at the dump were elimi- 
nated from the event sample. 

After the kinematics of the individual tracks were reconstructed at the target, 
they were combined to form dimuon pairs. Less than 0.08% of the reconstructed 
events contained more than two tracks. From the charge of the tracks and their 
momenta, the kinematics of the event were determined. Events which recon- 
structed to a mass < 2.0 GeV were rejected on the grounds that the acceptance 
for such pairs arising from the target was zero. 2 Events whose uniterated vertex 3 
indicated they were unlikely to have originated in the target were also cut. This 
cut was realized in terms of the z-coordinate of the uniterated vertex (ZUNIN), 
which was restricted to values less than 350 inches from the center of the target. 

3.3 Second Pass 

The goal of the first pass analysis was to reduce the data sample to a more man- 
ageable level. The purpose of the second-pass analysis was to reconstruct the 
kinematics of the events as accurately as possible, minimizing any systematic er- 
rors which would affect the cross sections. In the remainder of this chapter, we 

2 There was a problem in the application of this cut which applied to events with more than 
two tracks. If one of the pairs which could be combined from the tracks in the event was below 
the mass cut, the entire event was cut rather than the particular combination of tracks. Studies 
indicate that the error introduced by this is on the order of 0.1%, which is negligible compared 
to other uncertainties in the experiment. 

3 The position at which the tracks made their closest approach before the MCS corrections 
were made. 
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will discuss the differences between the algorithms the first- and second-pass anal- 
yses used to correct for multiple scattering and energy loss. We will also outline 
the procedures used to determine the strengths of the fields in the spectrometer 
magnets, and the position and angles of the beam at the target. 

3.3.1 Multiple Scattering 

The ZSCPLN was used to correct for the multiple scattering that the muons un- 
derwent in traversing the absorbing wall and beam dump as described above. The 
first pass-analysis used a single position for the ZSCPLN, which was empirically 
determined by minimizing the correlation between the mass of the dimuon pair 
and the mean ZUNIN. This approach ignored the additional scattering that events 
which passed through the dump were subject to. The second-pass analysis based 
the position of the ZSCPLN on the length of the dump which the muon pene- 
trated. It was found that for every 1 cm the track traveled through the dump, 
the ZSCPLN should be moved upstream 0.7 cm. 

The second pass analysis also employed a different technique for finding the op- 
timal position of the ZSCPLN. Monte Carlo studies indicated that the correlation 
between the mass of the events and ZUNIN was due to energy loss fluctuations in 
the beam dump and hadronic absorber, rather than difficulties in reconstructing 
the track angles as had previously been believed. The ZSCPLN in the second- 
pass analysis was placed at the position which optimized the reconstruction of the 
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track angles at the target for Monte Carlo events. This resulted in a ZSCPLN 
located 12 m downstream of the target, near the end of the carbon section of the 
hadron absorber. The method used in the first pass analysis placed the ZSCPLN 
22 m downstream of the target, well outside of the absorbing wall. 

3.3.2 Energy Loss 

As will be discussed below, determination of the strength of the magnetic fields 
from the data relied on an accurate reconstruction of the masses of the J/ip and 
T resonances, and the ZUNIN distribution of the events. Deviations between the 
average energy lost by the muons and the parameterization used to correct for the 
energy loss in the analysis could cause systematic shifts in both mass and ^-vertex 
reconstruction, and thus the measurement of the magnetic fields. 

Every muon which passed through the spectrometer was subject to a large 
amount of energy loss in the absorbing wall and possibly the beam dump. A 
typical 200 GeV muon lost on average ~ 2 GeV in the absorber. Figures 3.1 and 
3.2 show the energy lost in traversing the beam dump and absorbing wall for 
muons with different incident momenta. The calculations were performed using 
the TRAMU muon transport program [Gin 86]. The large, non-Gaussian tail 
skews the mean to much larger energy losses than the most probable value. It 
was important to parameterize the energy loss in the reconstruction of the data 
with emphasis on accurate (on average) reconstruction of the kinematics of the 
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events rather than optimal resolution. Thus, the mean energy loss was used in 
this analysis. 

Given the incident momentum of the muons on a specified material, the pa- 
rameterization of the mean energy loss is a well defined problem. However, the 
incident momentum of the muons on the absorbing wall or dump was the quantity 
which we wished to estimate from our measurement of the momentum downstream 
of SMI 2 and some parameterization of the random distribution of energy loss in- 
side the absorbing wall and dump. To further complicate the problem, muons 
which lost a large amount of momentum in the dump/absorber could be swept 
out of the acceptance. Muons with lower incident momenta would be more sus- 
ceptible to this, introducing an acceptance effect into the parameterization of the 
energy loss. 

Figure 3.3 shows the mean energy lost by muons in the 24" section of the ab- 
sorbing wall according to the Monte Carlo (where we knew the energy lost by each 
muon thrown). We plot the mean energy loss versus the momentum incident upon 
the 24" of copper, and versus the momentum after the 24" of copper. The "after" 
parameterization is generally flatter than the "incident" parameterization. This is 
due to the fact that larger mean energy losses in the "incident" parameterization 
will be shifted to smaller values of Paftcr- It should be clear from the figure that 
this effect is important, especially at lower momenta where we have most of our 
statistics in the data. 
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Figure 3.1: Energy loss for muons passing through the dump (144" of copper). 
Various incident energies are shown. 
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Figure 3.2: Energy loss for muons passing through the absorber wall. Various 
incident energies are shown. 



50 




Figure 3.3: Mean energy lost by muons in copper section of the absorber. Shown 
is the energy loss versus the momentum of the muon entering the section (circles) 
and the energy loss versus the momentum of the muon exiting the section (tri- 
angles). The solid line is the fit to the second set of points which was used in 
reconstructing the events. 
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Figure 3.3 also shows the fit to the energy loss in the 24" copper section of 
the absorbing wall. The parameterizations which we used for all of the absorber 
sections in the high-mass data are given by 

dEcu = 1-2202 + 0.2228 x 10~ 2 p a ftcr - 0.2472 x 10~ 5 ^ fter 

dE c = 0.2461 + 0.3135 x 10" 3 p aftcr - 0.4362 x 10~ 6 p aftC r 

dE CH2 = 0.1857 + 0.2073 x 10" 3 p aft c r - 0.3501 x 10- 6 p 2 ftcr 
where p a fter is the momentum immediately downstream of the absorbing wall sec- 
tion being traced through. Deviations in the shape of this parameterization, but 
not intercept, were found for the low- and intermediate-mass data sets. When used 
to reconstruct Monte Carlo events, systematic shifts in the masses and ZUNIN of 
the events are held to acceptable levels. 

3.3.3 Magnetic Fields 

The magnetic fields in the three spectrometer magnets were measured prior to 
the installation of the absorbing wall. Measurements were made at each of the 
planned operating currents, except for the 4000 A setting of the SM12 magnet. 
This setting was not measured due to concerns over the structural integrity of 
certain parts of the magnet under the increased strain of higher magnetic forces. 
At the time, it was anticipated that the 2800 A setting of the SM12 field would 
be the largest current at which that magnet would be operated. Later during the 
experiment, the SM12 magnet was successfully operated at the 4000 A setting, 
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providing the largest and cleanest sample of continuum dimuon events in the E866 
data sample. 

The magnet map used in the high-mass data utilized a measurement of the 
SM12 field from a previous experiment. However, the upstream portion of the 
magnet was reconfigured for E866. This region of the magnet map relied on an 
OPERA [Ope 94] calculation of the magnetic field, matched to the field map in 
the known region of the magnet. A similar procedure had to be performed to 
determine the magnetic fields in the beam dump, where it was not possible to 
measure the fields. 

Although the shapes of the magnetic fields were well determined, the absolute 
strengths of the fields were uncertain at the ~ 1% level. Monte Carlo studies 
indicated that this level of uncertainty in the SMO and SM12 magnets could be a 
significant source of uncertainty in the acceptance of the spectrometer, especially 
near its low mass edge. Figure 3.4 shows the ratio of Monte Carlo events gener- 
ated with magnetic fields differing by 5% using a fast version of the E866 Monte 
Carlo. The relative acceptance falls by nearly 40% over the given mass range. A 
systematic uncertainty of ~ 1% in the SM12 field would correspond to a ~ 8% 
point-to-point systematic uncertainty in the acceptance, which would have been 
unacceptable given the precision of our data. Furthermore, systematic shifts in 
the SM3 field would cause systematic shifts in the reconstructed kinematics of the 
events. Such shifts were a problem, considering the steeply-falling nature of the 
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Figure 3.4: Sensitivity of the acceptance to the magnetic fields. The ratio of 
Monte Carlo events generated with the nominal SM12 field to those generated 
with the SM12 field augmented by 5% is shown as a function of mass in the range 
0.5 <x F < 0.55. 

cross section we were trying to measure. 

The general procedure for determining the absolute field strengths from the 
data, once the impact of the energy loss parameterization had been limited, was 
as follows. We began by fitting the masses of any resonances (J/ip, T) in the mass 
spectrum to a Gaussian plus a quadratic background, noting the reconstructed 
mass of the peaks. We also estimated the mean of the ZUNIN distribution. Devia- 
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tions of these quantities from the nominal position of the target and the published 
resonance masses indicated an incorrect value for the SM12 and SM3 fields respec- 
tively. The data set was then re-analyzed with a multiplicative constant - referred 
to as a tweek (sic.) - applied to each of the magnetic fields to try to compen- 
sate. This was performed iteratively until the resonance masses converged to their 
published values, and ZUNIN was centered on the nominal target position. 

The procedure outlined above was essentially the procedure used to determine 
the tweeks on the SM12 and SM3 magnets in the high-mass data. It is a slight 
simplification, since a shift in the SM3 magnet would also have caused a shift in 
the ZUNIN of the events due to the increased momentum and resulting decrease in 
the deflection of the muon by SM12. The one important change to note was that 
only events which missed the dump were studied. Muons which hit the dump went 
through a somewhat uncertain amount of material, contributing an uncertainty to 
the amount of energy loss which should have been corrected for. Additionally, the 
ZUNIN of the events with one track in the dump was largely determined by the 
other track in the event. This reduced the sensitivity of ZUNIN to the magnetic 
field in SM12. 

Additional problems arose in the analysis of the low and intermediate mass 
data sets. In the intermediate-mass data set, we had only a limited number of 
events to work with. The reconstructed T mass, the only resonance available in 
that data set, was too uncertain to provide an accurate assessment of the SM3 
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tweek. Therefore, the value found in the high mass analysis was used for SM3, 
and the usual procedure of finding the SM12 tweek using the ZUNIN distribution 
performed. 

A similar problem occurred in the low- mass data set, where we had an extra 
magnetic field (SMO) to worry about. In principle, since we could resolve both 
the J/ip and ip' resonances, we had the extra information needed to tune the SMO 
field. In practice, the ip' mass was relatively insensitive to the SM3 field, for a 
reasonable set of tweeks. Once again, the SM3 field tweek from the high-mass 
data was used in this data set. The SMO tweek was found from the fits to the J/ip 
mass, and the usual ZUNIN fitting procedure used to extract the SM12 tweek. 

The procedures described above provided an estimate of the average field 
strength in a given data set. It also provided a way to estimate the systematic 
uncertainties on the tweeks by studying the impact of a given change in tweek 
on the ZUNIN distribution and the resonance masses. Uncertainties in the mass 
and vertex fits could then be translated into uncertainties in the magnetic fields. 
Once the uncertainties in the magnetic fields were determined, the correspond- 
ing acceptance uncertainty could be calculated by running the Monte Carlo using 
different input fields. We will describe this in greater detail in the next chapter. 

Figure 3.5 shows the J/ip and T mass regions from the low- and high-mass 
data, when reconstructed using the optimal magnetic field tweeks. The recon- 
structed J/ip and T masses of 3.095 and 9.461 GeV compare quite well with the 
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Figure 3.5: Reconstructed J/ip and T resonances. Left panel shows the low- 
mass data, right panel the high-mass data. Data is the sum of the hydrogen and 
deuterium targets. 

published values of 3.097 and 9.46 GeV [Gro 00]. The intrinsic widths of the res- 
onances (which are best expressed in keV), are smeared out by energy loss and 
multiple scattering in the absorbing wall and dump - our mass resolution is 92 
MeV at the J/ip mass and 183 MeV near the T. Figure 3.6 shows a fit to the 
ZUNIN distribution of the continuum dimuon regions, where the data has been 
summed over all data sets. In practice, fits to the ZUNIN distributions and res- 
onance regions were performed separately in each data set. The deviations of 
the fits from their expected values amount to uncertainties of only ±0.10% and 
±0.15% in the SMO/12 and SM3 tweeks, respectively. 

One other source of systematic error in the tweeks was uncovered using the 
Monte Carlo. When Monte Carlo events were retracked through the reconstruc- 
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Figure 3.6: Uniterated z-vertex distribution. Data is the sum of both targets, 
corrected for empty-target background. The fit is to a Gaussian curve. 

tion algorithm, the momentum found at the SM3 bend plane was ~ 0.06% larger 
than that which was actually generated by the Monte Carlo. Thus, the correct 
field tweeks were subject to a small offset which had to be accounted for to obtain 
a more accurate simulation of the spectrometer. The origin of this effect was un- 
clear. It may have been due to slight differences between the actual and simulated 
geometry of the apparatus, resolution effects coupled to the single bend plane ap- 
proximation used for SM3, or some combination of these and other issues. This 
additional source of systematic uncertainty in the SM3 field does not contribute 
significantly to the overall size of the systematic uncertainties in the acceptance. 

In table 3.2 we list the tweeks found for each mass setting using the previously 
described procedures. Studies of the stability of the mean ZUNIN and T masses 
versus data set and run number in the high-mass data suggest that a single value 
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Table 3.2: Field strengths for each data set. 



Mass setting 


Data set 


SMO tweek 


SM12 tweek 


SM3 tweek 


low 


5 


1.0140 ±0.0010 


1.0068 ±0.0010 


1.0049 ±0.0015 




10 


1.0140 ±0.0010 


1.0068 ±0.0010 


1.0049 ±0.0015 


int 


9 


N/A 


1.0093 ±0.0010 


1.0049 ±0.0015 




7 


N/A 


0.9865 ±0.0010 


1.0049 ±0.0015 


high 


8 


N/A 


0.9865 ±0.0010 


1.0049 ±0.0015 




11 


N/A 


0.9865 ±0.0010 


1.0049 ±0.0015 



for the SM3 field tweek is reasonable. These studies also indicated that the SM12 
field was reliably stable, and a single value was used for each data set within a 
given mass setting. 

3.3.4 Beam Position and Angles 

The position of the beam was determined by fitting the distribution of the x 
and y vertices of the raw data at the target position. The uniterated x and y 
positions of the two tracks in the event were averaged together to yield the x and 
y vertices. Run-to-run variations in beam position were found to be negligible for 
our purposes. 

The beam angle was found through an iterative process involving the Monte 
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Carlo. In principle, the Y component (and similarly the X component) of the 
beam angle could be found from just the data by looking at the muon momentum 
vectors Q^ muon = (p+ + p-y(p+ + p~). However, asymmetries in the dump cut 
in the first pass analysis limited the accuracy of this technique. An independent 
test of the validity of the beam angle extracted from the data was to check the 
(j) p distribution in the data against that of the Monte Carlo. Since we know that 
the 4> p distribution must be isotropic, and having generated an isotropic distri- 
bution in the Monte Carlo, any difference between the distribution in the data 
versus the Monte Carlo could only result from a difference in the beam angle be- 
tween the experiment and the Monte Carlo. The beam angles were improved by 
minimizing these differences, and we estimate that the uncertainty in the beam 
angle contributes no more than 5% (~ 0.7 GeV) of the transverse momentum of 
an average event. 
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4 ANALYSIS 

The differential cross section is defined as the number of interactions per target 
particle leading to an event with kinematics between Q and Q + Sfl per number 
of incident particles per unit area. We express this as 

*Z N ^ (41) 
dtt N T {Ni/A)' v ' ; 

where N is the number of dimuon events from the target, N T the number of target 

particles in the path of the beam, Ni the mean number of protons incident on the 

target and A the cross sectional area of the beam. 

The simple expression in equation 4.1 is somewhat misleading, since we did 

not measure the number of pp and pd interactions which led to a dimuon event 

directly. Instead, we measured the yields of oppositely charged muon pairs from 

the interaction of a proton beam with target flasks containing liquid H 2 and 

D 2 . These yields were subject to contamination by background events originating 

from the interaction of the beam with the SWIC and the front and back faces 

of the target flasks, and from the random coincidence of two uncorrelated muons 

produced in the same beam bucket. The yields were also subject to losses due to 

the acceptance of the spectrometer and inefficiencies in the detection of the muons. 

Equation 4.1 also relies on quantities such as the number of target particles and 

incident protons, which were also not directly measured, and the cross sectional 

area of the beam, which was not measured with any great precision. 
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Fortunately, we can rewrite equation 4.1 entirely in terms of quantities which 
were directly measured, while eliminating terms which were not measured well. 
We define the number of dimuons, corrected for backgrounds, as iV M + M - . This can 
be related to the number of interactions in equation 4.1 by 

jytargct jybackground 

ae ae 

where a is the geometric acceptance of the spectrometer, defined as the probability 
that both tracks in a given event would traverse all of the active detector elements 
of the spectrometer and survive all data cuts, e is the efficiency with which the 
given event was detected, and the "target" and "background" superscripts refer to 
events which originated in the liquid target and background events respectively. 
The number of target particles can be expressed in the usual way: 

N T = N A pLA (4.3) 

where Na is Avogadro's number, p is the density of cryogenic liquid, L is the 
length of the target, and A is the cross sectional area of the beam. 

To get the mean number of incident protons, we need to average over the 
attenuated flux in the target. The number of protons at any point along the 
length z of the target can be expressed in terms of the incident number of protons 
and the hadronic absorption length A, as 

Ni(z) = N p e~ z/ \ (4.4) 
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Averaging over the total length of the target L yields 

Ni = y N t (z) dz= 1 -^{\- e- L ' x ). (4.5) 
Li Jo L, 

Substituting equations 4.2-4.5 into equation 4.1 results in an expression for 

the differential cross section: 

da = iV> M - 

dn CaeAtt' 1 ' ; 

where we have grouped several of the factors in the denominator into a single 
quantity, called the luminosity, given by 

C = N A p\(l-e- L/x )N p . (4.7) 

We will now explain in detail how each of the quantities in equation 4.6 was 
determined. 

4.1 Event Yields 

The total event sample produced by the second pass analysis contained several 
questionable events. These included backgrounds from uncorrelated muon pairs 
and correlated pairs originating outside of the target volume, events which were 
reconstructed by the analysis but were not responsible for the trigger which caused 
the DAQ to read out, and events which occupied regions of the spectrometer which 
complicated the reconstruction of their kinematics. These events were either fil- 
tered out of the second pass event sample, or independently binned and subtracted 
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from it, to obtain the event yields used to calculate the cross sections. The same 
set of cuts were used on the Monte Carlo events used to calculate the acceptance. 
Simulations showed that contributions from the interactions of secondary hadrons 
produced in the target were negligible. 

In order to avoid the J/ip and T resonances, we selected events within the 
mass range 4.2 < M M + M - < 8.7 GeV, and events with masses M > 10.85 GeV. As 
the acceptance above the T resonance was quite small in the low and intermediate 
mass data sets, we only used the high mass data in that region. Due to uncer- 
tainties in the magnetic fields, data from the intermediate mass setting below 6.2 
GeV were also discarded. Figure 4.1 shows the total E866 mass spectrum, with 
the contributions from the low-, intermediate- and high-mass data sets in their 
respective allowed mass ranges. The resonance regions have been retained, but 
all other cuts were applied. The number of events surviving all cuts is shown in 
table 4.1. 

4.1.1 Target Flask Background 

One of the larger sources of contamination was due to the beam interacting with 
the SWIC (rs 188 cm upstream of the target) and the walls of the target flasks 
which contained the cryogenic liquids. Figure 4.2 shows the yields from the cryo- 
genic liquid and empty target flasks plotted versus ZUNIN. The SWIC is clearly 
visible in the empty target distributions, and much of the contamination was re- 
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Figure 4.1: Total E866 mass spectrum surviving all cuts. The contribution of 
the low, intermediate and high mass data sets superimposed in their respective 
allowed mass ranges. Inset figures show the low, intermediate and high mass data 
sets on a linear scale over their entire mass range. 
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Table 4.1: Number of events in each data set which survived all data cuts. 



Set 


H 2 target 


D 2 target 


Empty target 


5 


17217 


36962 


699 


7 


9548 


21030 


425 


8 


20049 


44129 


1007 


9 


1768 


4966 


63 


10 


1658 


3509 


60 


11 


4850 


10545 


230 


total 


55090 


121241 


2484 



moved by simply applying cuts (ZUNIN > —50" in the low mass data, and ZUNIN 
> —60" in the intermediate and high mass data). Events with ZUNIN > 60" in 
the low mass data, and ZUNIN > 90" in the intermediate and high mass data 
were also rejected, as they were likely to have originated in material downstream 
of the target. 

The remaining background from the target flask was measured during the 
experiment, allowing us to subtract it from our event yields. A correction was 
made to the empty target yields based on the uniterated ^-vertex of the event to 
account for the ~ 7% and 14% attenuation of the beam between the front and 
back end caps of the flasks containing liquid hydrogen and liquid deuterium. This 
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Figure 4.2: Comparison of target yields to empty-target background. The circles 
represent the sum of the hydrogen and deuterium yields, the histogram represents 
the normalized empty-target contamination. 
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correction was given by 

f N empty x e- L '\ if ZUNIN > 

N ernpty = J ^ ^ 

1 7V empt2/ , if ZUNIN < 

where we corrected events which appeared to come from downstream of the nom- 
inal target position z > by the attenuation of the cryogenic liquid appropriate 
to the target. 

4.1.2 Combinatoric Background 

Another important background source was the random coincidence of oppositely 
charged, uncorrelated muons produced in the target in the same beam bucket. 
The numbers and kinematic distributions of these random events were estimated 
using event samples triggered on the like-sign (PhysA2) and single-muon (PhysB2) 
triggers. These random events were typically confined to masses below 5 to 6 GeV 
and were only important in the low- and intermediate-mass data sets, the latter 
being a non-issue due to the 6-GeV cut in the intermediate mass. 

The single-muon trigger was used to trigger on what appeared to be single 
muons from the target traversing one of the matrix elements. Due to the high rate 
of such events from the target, these triggers were prescaled when in use during 
normal run conditions. Single-muon event samples were taken either during the 
normal course of the data collection or in special low-intensity runs. The resulting 
event samples were subjected to the same analysis chain as the opposite-sign data. 
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The single muons were then randomly combined with other single muons to form 
an uncorrelated dimuon event. This produced a sample of random dimuon events 
which, when properly normalized, could be subtracted from the data sample. 
Figure 4.3 compares the total dimuon yield in the data (summed over both targets 
and the high and low mass data sets) to the combinatoric background. The low 
mass data set contributes the largest contamination - approximately 5% of its 
total yield was due to randoms. 

The normalization of the randoms was accomplished by comparing a subset 
of the randoms (those with two muons with the same charge), to events which 
satisfied the like-sign trigger in the data. These like-sign events were analyzed 
in the same way as any other event, except that one of the tracks was reflected 
about y = when the kinematics of the events were calculated. Figure 4.4 shows 
the real like-sign events compared to the combinatoric like-sign events, scaled to 
match the real like-signs. 

4.1.3 Trigger Cuts 

For any given event, the possibility existed that the muons which fired the trigger 
and resulted in the event being read out by the DAQ were not the same muons 
which were reconstructed by the analysis. The easiest way for this to occur was to 
have a single, beam-like muon fire a set of hodoscopes, one of which was outside of 
the matrix, while simultaneously a second, low-momentum muon fired the missing 
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Figure 4.3: Comparison of the target yields (dots) to the combinatoric back- 
ground (histogram). The top panel shows the unsubtracted mass distributions, 
with the inset showing the same distributions down to 2 GeV. The bottom panel 
shows the fractional contribution of randoms vs. mass. 
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Figure 4.4: Comparison of the data (circles) and combinatoric (histogram) like- 
sign events. The data and combinatoric like-signs are summed over the low and 
high mass data sets and both targets. The inset figures show the mass, xf and 
Pt like-sign distributions over the mass range used in this analysis. 
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matrix element. When combined with a second, uncorrelated, oppositely charged 
muon which fired a valid matrix element, this would appear to be a high pr 
continuum dimuon event. We therefore had to check that the muons in each 
dimuon event satisfied the trigger condition which had caused the reading out of 
that event. This was done by comparing the reconstructed positions of the muons 
at each hodoscope plane to the aperture of the hodoscopes in the trigger which 
fired. If the track was outside that aperture (which was widened slightly so that 
we would not lose events due to chamber resolution effects), it was discarded. 
Events which fired PhysA2 (the like-sign trigger) were binned separately from 
events which fired the other PhysA triggers. 

In addition to the trigger cuts outlined above, comparisons of the data and 
Monte Carlo suggested that there was a problem with the trigger hardware cor- 
responding to certain matrix elements. To first order, the acceptance of the 
spectrometer should be left-right symmetric. A separate analysis [Bro 01], which 
studied events from the dump in the high mass data, noticed that several ma- 
trix elements yielded much more data relative to Monte Carlo on one side of the 
spectrometer than the other. This discrepancy was much larger than could be 
explained by left-right asymmetries in the spectrometer, and thus these matrix 
elements were suspect. Similar studies were performed for the target data, con- 
firming that there was a problem with several matrix-elements. These matrix 
elements, which tended to be at the edge of the acceptance where the statistics 
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were generally poor. Events with muons occupying these elements were cut from 
the final event sample. 

4.1.4 Event Geometry 

Several events which passed through all of the second pass cuts occupied regions 
of the spectrometer which complicated the reconstruction of their kinematics, or 
were judged to raise questions about their origins. Events which had a muon 
which appeared beam-like (tracks with \0 t y ar '9 et \ < 0.006 or tracks which were on 
the wrong side of the dump given their charge) were discarded. Such tracks 
could have originated from the decay of a hadron produced in the dump, and also 
followed a trajectory which skimmed the edge of the dump, calling into question 
the amount of energy lost by the track. 

Another problematic class of events was found by studying the Monte Carlo. 
The energy lost by tracks in the Monte Carlo was plotted versus the x position 
of the tracks at SM3. Tracks near the edge of the acceptance, |xsm3| > 24", had 
significantly more energy loss in them than tracks with |xsm3| < 24". The source 
of this excess energy loss was the walls and teeth of the SM12 magnet, indicating 
that this class of tracks passed through a region of the magnetic field which was 
highly uncertain. 

As discussed in the previous chapter, energy-loss fluctuations and multiple 
coulomb scattering cause the events to appear to reconstruct away from the actual 
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Table 4.2: Spill quality cuts. Events were kept if they obeyed the inequality. 





mass setting 


Cut 


low intermediate high 


SEM 
IC3 

AMONSB / AMON 
DUTY 


> 150 > 100 > 100 

> 5000 > 5000 > 5000 

> 0.9 > 0.9 > 0.9 

> 25 > 50 > 25 



interaction point in the target. To limit the fluctuation of energy loss and the 
effects of multiple scattering in the data, and to provide additional discrimination 
against events which originated outside the target volume, we required the event 
to reconstruct to within 2.5" of the nominal transverse position of the beam. 

4.1.5 Spill Quality Cuts 

During each beam spill, a number of quantities were monitored which told us the 
overall quality of the data being collected. These included the intensity of the 
beam (SEM and IC3), the luminosity of the target (AMON), the live-time of the 
spectrometer (AMONSB/ AMON) and the beam duty factor (DUTY). Table 4.2 
lists the spill quality cuts which were required for each spectrometer setting. 
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4.2 Luminosity 

The intensity of the beam in each spill was measured with the SEM6 counter on 
the Meson East beam line. This recorded a SEM count which was proportional 
to the number of protons in the spill. The SEM6 counter was calibrated during 
special calibration runs where a thin copper foil was inserted into the beam and 
a SEM count accumulated during a number of spills. By studying the rate of 
the 1368 keV 7's emitted by the 24 Na produced in the interaction, the number of 
protons incident on the foil could be deduced, and the SEM response calibrated. 
These calibrations rely upon a measurement of the 24 Na production cross section 
[Bak 84] of 3.90 ±0.11 mb 1 . 

Table 4.3 shows the results of several calibration measurements dating back 
to the original E605 experiment. Based on these measurements (taking the more 
recent measurements with greater weight), we estimate a SEM response of (0.79 ± 
0.051) x 10 8 protons / SEM[Bro 97]. The uncertainty in the normalization was 
determined by adding the uncertainty due to the 24 Na cross section in quadrature 
with the standard deviation of the measurements in table 4.3. 

Studies of the linearity of the SEM response were also performed by comparing 
the SEM counts to the luminosity as measured by the AMON and WMON detec- 
tors, and comparing to a redundant measurement of the beam intensity using the 

1 Although reference [Bak 84] measures the 2A Na cross section using 400-GcV protons, it is 
generally believed that the cross section is independent of energy. 
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Table 4.3: SEM calibration measurements on the Meson East beam-line. 



Date 


Protons 


Protons/SEM xlO 8 


01/26/84 


2.8E13 


.83 


02/11/85 


6.9E12 


.78 


03/25/85 


2.0E12 


.85 


05/06/85 


2.0E12 


.90 


07/25/85 


9.9E13 


.80 


08/29/85 


9.6E13 


.80 


11/10/87 


3.4E13 


.79 


11/16/87 


2.4E13 


.77 


02/11/88 


3.6E13 


.78 


10/29/96 


1.7E12 


.885 


10/29/96 


8.5E11 


.915 


12/20/96 


5.1E13 


.76 


08/15/97 


5.0E11 


.86 


08/29/97 


2.0E11 


.77 
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ion chambers. These studies showed the SEM response was linear with increasing 
beam intensity, but subject to a constant offset which varied from run to run. 
These offsets, shown in table 4.4, were subtracted from the SEM counts before 
calculating how many protons were incident on the target in each spill. 

Both the density and attenuation lengths of the targets depended upon the 
exact composition of the targets. Therefore, assays of the cryogenic liquids were 
performed using samples taken when the target flasks were emptied for the Christ- 
mas shutdown (end of data set 7) and at the end of the planned run (end of data 
set 11). The cryogenic liquids used in the first fill had negligible amounts of con- 
taminants. The same was true for the liquid hydrogen in the second fill. However, 
two separate analyses of the second fill's liquid deuterium showed a significant HD 
contamination. Table 4.5 shows the measured contamination and the estimated 
nucleon content in the target. Once the final pp cross sections had been measured, 
the deuterium cross sections were corrected for the ~ 3% contamination. 

In order to determine the target densities, the vapor pressures of the cryogenic 
targets were measured and periodically recorded. The high circulation rate of the 
liquids, combined with the comparatively low beam current ensured that the cryo- 
genic liquids did not boil during the spill which would have reduced the densities 
of the targets. Table 4.6 shows the mean target pressure for selected data sets in 
the experiment. 

Cryogenic data tables [Rod 73] were used to calculate the target densities from 
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Table 4.4: SEM offsets. These were subtracted from each spill. I Tun is the run 
number. 



Starting run 

o 


Ending run 

o 


SEM offset 


X OO X 


1 ^88 


97 


1 ^89 


1 620 


1 479 x T — 2436 


1 631 

X \J*J X 


1 690 


28 


1 87Q 


1 884 


98 

Z/O 


1 88fi 
loot) 


1 887 
loo ( 


100 


1888 




275 


1889 


1944 


375 


1945 


2020 


130 


2021 


2046 





2049 


2053 


98 


2056 


2078 


100 


2079 


2115 


30 
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Table 4.5: Percent molecular and atomic abundance of second deuterium fill. 



Molecule 


percent abundance 


percent abundance 




(first analysis) 


(second analysis) 


D 2 


94.f3±0.58 


92.7 ±0.8 


HD 


5.82 ±0.58 


6.89 ±0.69 


H 2 


0.05±0.0f 


0.f47±0.0f5 


other 






Atom 


estimated abundance 


D 


97.00 ±0.6 


H 


3.00 ±0.6 
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Table 4.6: Target pressures (in psi) for selected data sets. 



Data set 


hydrogen 


deuterium 


5 


14.97 


14.92 


7 


15.04 


14.96 


8 


15.11 


15.17 


11 


15.15 


15.21 



the average vapor pressures. Taking the pressure P to be measured in psi, and 
density p in g/cm 3 , the expression for the density of hydrogen is given by 

— = 62.473 {0.2115 + 1.171 x 10~ 3 P - 1.109 x 10" 5 P 2 }, (4.9) 
Ph 

and that for deuterium by 

p d = 4.028 x 10~ 3 {43.291 - 3.4176 — - — + 0.5783( — - — ) 2 } (4.10) 
1 1 14.6959 V 14.6959 ; 1 V ; 

These result in densities which vary from data set to data set by no more than 
~ 2%, which is negligible compared to the 6.5% uncertainty in the beam normal- 
ization. From the measured pressures, the densities of the hydrogen and deuterium 
targets were determined to be 

p h = 0.0706 g/cm 3 
pa = 0.1630 g/cm 3 . 
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As reported in the particle data book [Gro 00] , the hadronic interaction lengths 

2 2 

for hydrogen and deuterium are 50.8 g/cm and 54.7 g/cm respectively. To ex- 
press these lengths in terms of the physical length, as opposed to the "thickness" 
of the target, we need to divide by the density of the targets. For the cryogenics 
used in the first fill, and the liquid hydrogen from the second fill, this is all we 
needed to do. For the liquid deuterium from the second fill, however, we needed 
to correct for the hydrogen contamination by taking the weighted average. The 
resulting hadronic absorption lengths are 

Mh2 = 719.5 cm 

Kdf = 335.6 cm 
Ag 2 nd) = 340.6 cm 

where the superscripts indicate either the first or second fill. The total integrated 
luminosities for each data set and the entire data sample are shown in table 4.7. 

4.3 Trigger Efficiencies 

In order to properly measure the efficiencies of the hodoscopes, it was necessary 
to use a sample of events which were triggered independently of the hodoscopes 
being measured. Because we only use the information from the Y hodoscopes in 
the experiment, events which triggered on the X hodoscopes could provide this 
independent measurement. Unfortunately, the rate of events which satisfied the 
X134L-X134R trigger (a muon on either side of the spectrometer) was prohibitively 
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Table 4.7: Integrated luminosities in the E866 data sets. 



Data Set 


Integrated H 2 


Integrated D 2 




Luminosity (nucleon/nb) 


Luminosity (nucleon/nb) 


5 


4.74 x 10 6 


1.04 x 10 7 


7 


1.97 x 10 7 


4.34 x 10 7 


8 


4.04 x 10 7 


8.95 x 10 7 


9 


2.95 x 10 6 


6.65 x 10 6 


10 


4.72 x 10 5 


1.03 x 10 6 


11 


1.03 x 10 7 


2.28 x 10 7 


Total 


7.86 x 10 7 


1.74 x 10 s 
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high to take along with the physics data. Therefore, several specialized runs were 
taken periodically to determine the efficiencies of the Y hodoscopes. A specialized 
trigger required that five hodoscope planes on one side of the spectrometer were 
hit. Reconstructed single muon tracks which fired at least six hodoscopes (one 
in each plane) were then studied to determine the efficiency of the hodoscope in 
the seventh plane. The typical efficiencies measured were ~ 97%. However, these 
specialized runs did not measure the efficiencies near the edge of the acceptance 
with sufficient precision. To increase the statistics in these regions, we used events 
from the full data set. Although these events introduced a bias from the trigger 
to the efficiency measurement, they did probe the edge of the acceptance with far 
more statistics. The efficiencies measured here were checked against those of the 
specialized runs and found to be largely consistent. 

4.4 Chamber Efficiencies 

In order to study the efficiency of each detector plane, we began with a sample 
of tracks which would have been reconstructed regardless of whether or not the 
particular detector plane being studied had fired or not. This was possible due to 
the redundant tracking information at each station. Each plane was divided into 
quadrants, and the ratio of the number of tracks which fired in a given quadrant 
to the total number of tracks passing through the quadrant was used to estimate 
the efficiency of the drift chamber. 
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Table 4.8: Rate-dependent corrections to the hydrogen and deuterium yields. 
Each event is weighted by 1 — SEM x where % is the target. 



Mass Setting 


Rd 


F 


Rh = Rd/F 


Low 


-0.2035 x 10~ 4 


1.51 


-.1348 x 10~ 4 


Intermediate 


-0.4976 x 10~ 5 


1.18 


-.4152 x 10~ 5 


High 


-0.2696 x 10~ 5 


1.78 


-.1515 x lO" 5 



The use of single-hit TDC's in the DAQ electronics introduced an additional 
source of inefficiency A second charged particle which passed through a drift cell 
following the first particle in an event would not have its position recorded to 
tape, resulting in the event being more difficult to reconstruct. This effect would 
have a definite dependence on the event rate, and therefore the beam intensity. 

To correct for the rate dependence, the number of reconstructed events per 
unit beam intensity (measured in terms of SEM counts) was plotted versus the 
beam intensity and fit to a straight line. The fits for hydrogen and deuterium 
were constrained by the relative number of extra drift chamber hits in an average 
event F = j^kx™ for each target [Tow 99]. This meant that only one parameter, 

h 

the slope Rd of the efnciency-versus-intensity plot for deuterium, needed to be 
determined for each mass setting. Figure 4.5 shows the fit for the high mass data, 
and table 4.8 shows the corrections obtained by reference [Tow 99]. 
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350 - 



300 - 



100 - 



Data Set 



X 2 / ndf 


76.89 / 72 


R d 


-.2696E-05 


D 7 


309.0 


D 8 


313.5 


D„ 


293.0 


E 


.4467 




11 



250 500 750 1000 1250 1500 1750 2000 



Units of Beam Intensity (s) 



x 10 



Figure 4.5: Rate dependent inefficiency. Dimuon yield (N m ) per unit beam 
intensity (s) versus beam intensity for the high mass data, as given in reference 
[Tow 99]. 
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4.5 Monte Carlo and Acceptance 

The acceptance of the spectrometer was a complicated function involving the 
detector geometry, the magnitudes of the magnetic fields in each of the three 
spectrometer magnets, the efficiencies of the detectors and trigger, and the various 
physical processes which muons traversing the spectrometer were subject to. This 
level of complexity required the use of a Monte Carlo simulation to determine the 
acceptance of the spectrometer. 

The E866 Monte Carlo consisted of custom routines designed to simulate the 
production of dimuon events, the passage of the resulting muons through the 
spectrometer, and the detector response to their passage. The output of the Monte 
Carlo was a file whose format was identical to that produced by the DAQ during 
the experiment. This enabled us to analyze the Monte Carlo events with the 
exact same analysis chain as the real data, minimizing any systematic differences 
between the data and Monte Carlo. 

To further minimize these differences, the Monte Carlo was configured to sim- 
ulate as closely as possible the state of the spectrometer while it was taking data, 
the physical processes the muons were subject to in passing through the spec- 
trometer, and the physical processes giving rise to the dimuon events. In chapter 
3 we discussed how we obtain information about the state of the spectrometer 
(magnetic fields, detector alignments and efficiencies, etc.). We will discuss how 
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events were generated in the Monte Carlo, the simulation of the physical processes 
the resulting muons were subject to, and the final set of calculations needed to 
compute the acceptance of the spectrometer. 



4.5.1 Event Generation 

In order to generate a continuum dimuon event, we had to sample the six kine- 
matic variables M^^-, xp, px, (fi p , (fid, an d 9 d from distributions which closely 
approximate the actual differential cross section d 6 a / ' dMdxFdpTd6dd<fi p d<fid- A 
priori we know the cross section differential in the angular variables goes like, 

— ^-—ocl + cos 2 (^) (4.11) 
d6 d d(fidd(fi p 

since there is no preferred direction in the experiment to give a (fi dependence, 
and the virtual photon in a qq annihilation is produced transversely polarized 
as predicted by QED and confirmed by experiment [Con 89, McG 99, Cha 99, 
Bro 01]. 

Next-to-leading order calculations of the doubly-differential cross section d 2 a/ dMdxp 
agree well with existing measurements of the deuterium and nuclear cross sections, 
but problems arise in QCD in calculating the full triple differential cross section 
d 3 a / ' dMdxpdpT- Therefore, an empirical form [Kap 78] for the pr dependence of 
the cross section was used 

^empirical 1 10p T /p T C^NLO 

dMdxpdpT K p T (1 + (pt/Pt) 2 ) 6 X dMdx F 1 ' ' 
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where p? = 2.8 GeV was fit to both the hydrogen and deuterium data, and the 
MRST 98 [MRST 98] partons were used in the NLO calculations of d 2 a/dMdx F . 

4.5.2 Physics Simulation and Detector Response 

Once an event had been sampled from the above distributions, the kinematics of 
each muon were calculated and boosted into the lab frame, and the event vertex 
sampled from 

= jr^m- (4 ' 13) 

Inside the target, the very small amount of multiple scattering the muons were 
subject to was simulated. This multiple scattering was found to be negligible. 

After the effects on the muons in the target had been simulated, the muon's 
positions were extrapolated to the point at which the SM0/SM12 field began. 
The same magnetic field maps which were used in the analysis as well as the 
overall field strengths as computed from the data were used here. The muons 
were propagated through the magnetic field one at a time in 5 cm steps. At each 
step in the field, the magnetic field kick was calculated, using 

dp x = -qdzB y (4.14) 
dp y — q dz B x (4-15) 

which neglected only the small z component of the force. 2 



2 To be precise, we computed the changes in the x and y components of the muon's momentum, 
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In addition to the simulation of the magnetic fields, the effects of the beam 
dump and absorbing wall had to be included. Multiple scattering was simulated in 
each 5 cm step of the dump and absorbing wall, using the gmols routine from the 
standard GEANT [GEANT] simulation package. This was the only case where 
a GEANT routine was used in the Monte Carlo. Multiple scattering was also 
simulated in the materials found in the helium bags which filled SMO and SM12 
downstream of the absorber, the air and helium the muons had to pass through, 
the hodoscope and drift chamber planes, and the bulk material at station 4. 

The energy lost by muons traversing the spectrometer was sampled from a set 
of fast lookup tables generated using Von Ginniken's TRAMU program [Gin 86]. 
Twelve tables for muons with incident momenta between 10 and 600 GeV were 
interpolated to obtain the energy loss for muons of arbitraty momenta. This num- 
ber was then scaled from the length of material assumed in the table to the step 
size of 5 cm, and subtracted from the muon's momenta in each step. A new energy 
loss was sampled once the muon had traversed the length of material specified by 
the energy loss routines. This greatly simplified the energy-loss simulation in the 
dump, where the total length of material the muon would pass through was not 
known beforehand. The only other places where energy loss was simulated was in 

the teeth and walls of the SMI 2 magnet, and the hadronic and EM calorimeters 

then applied momentum conservation to calculate the final z component of the momentum. It 
can be shown that this causes a systematic error in the momentum at each step which is negligible 
over the entire magnetic field. 
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at station 4. A precise simulation of the energy loss was unnecessary here, so a 
constant energy loss was subtracted in these materials. 

Equation 4.6 indicates separate corrections for the acceptance of the spec- 
trometer and detector efficiencies. In reality the Monte Carlo simulated detector 
inefficiencies based on the measurements described above. Thus, the acceptance 
calculations contained the efficiency corrections. Only the rate-dependent ineffi- 
ciencies were corrected for separately. 

4.5.3 Acceptance Calculation 

Once a sufficient number of Monte Carlo events had been recorded, they were 
subjected to the same analysis chain and cuts described in chapters 3 and 4. The 
only difference was that spill quality cuts were not applied to the Monte Carlo. 
Once we had filtered the Monte Carlo event sample, the acceptance in any given 
bin could be calculated as the ratio of the number of surviving Monte Carlo events 
in that bin divided by the total number of events generated in that bin 

N^ c c (M,x F ,p T ) 

a(M,x F ,p T ) = gen r. (4.16) 

^mc( m . %f,Pt) 

Calculating the number of accepted and generated Monte Carlo events in each 
bin would have been a straightforward counting exercise, save for two reasons. 
First, comparisons of the data to the Monte Carlo revealed differences on the order 

of 10 to 15% which we wished to eliminate. 3 Second, in order to eliminate the 

3 These differences were consistent with the observed differences between the measured cross 
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systematic error between the hydrogen and deuterium cross sections introduced 
by finite Monte Carlo statistics in the acceptance calculation, the same sample 
of Monte Carlo events were used to calculate the acceptance for hydrogen and 
deuterium. Since the event sample had been generated using either a hydrogen or 
deuterium target, the events had to be reweighted by the cross section appropriate 
to the target being calculated. 

Differences between the data and Monte Carlo were reduced by reweighting 
the Monte Carlo events as a function oi M, x F and p T . These functions were 
determined by polynomial fits to the ratios of the data divided by the Monte 
Carlo. Furthermore, the p F cross section was seen to vary as a function of M 
and x F . This correlation was dealt with by fitting in each bin, resulting in the 
empirical form 

p° T (M, x F ) = 1.21 + 0.350M - 0.0182M 2 + 1.37x F + 2.60x 2 F (4.17) 

which was the same for both targets, within the uncertainties in the fit parameters. 
The resulting weighting functions for the Monte Carlo were consistent with the dif- 
ferences observed between the measured cross sections and the double-differential 
cross section tables used to generate the Monte Carlo. 

Monte Carlo events which had been generated using hydrogen as the target 

were reweighted with the deuterium cross section and vertex distribution when 
sections, and the NLO calculations upon which the Monte Carlo generator was based. 
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calculating the deuterium cross section. Thus, the deuterium yield calculated 
from a hydrogen Monte Carlo event sample was given by 



Hydrogen yields were calculated from deuterium Monte Carlo samples in a similar 
manner: 



This procedure was carried out for both the accepted and generated events, main- 
taining the relative normalization of the thrown and accepted Monte Carlo, and 
equation 4.16 was used to compute the acceptance. 

Figures 4.6 and 4.7 compare the data and Monte Carlo for several recon- 
structed pair- and detector-quantities, summed over both targets and averaged 
over all data sets. The agreement between the data and Monte Carlo is quite 
good for most physics and detector distributions. The discrepancies in the single 
track momentum distributions, which occur at larger momenta where the statis- 
tics are poor, are thought to be due to beam-alignment problems. The other 
major discrepancy is in the ZUNIN distributions. This is partly due to a slightly 
wider ZUNIN distribution in the Monte Carlo than is observed in the data, but 
mainly due to the centroids of the data and Monte Carlo distributions not being 

precisely aligned. 4 

4 This is a major component of the estimated systematic errors discussed below. 




(4.18) 




(4.19) 
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Figure 4.6: Comparison between data and Monte Carlo yields. The ratios of 
data to Monte Carlo yields versus several "pair" quantities (mass, xp, Pt, <f>d, 
(ftp, cos(9d) and ZUNIN) and the positive and negative track momenta. Ratios 
averaged over all targets and data sets. 
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Figure 4.7: Comparison between data and Monte Carlo yields. The ratios of 
data to Monte Carlo yields versus the positions of the positive (negative) track 
projected onto the station 1 hodoscope-, SM3 bend- and station 4 hodoscope- 
planes. Ratios averaged over all targets and data sets. 
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4.6 Combining the Data 

Once all of the data were binned and subjected to all of the corrections de- 
scribed above, we had a measurement of the triple-differential cross sections 
d 3 a(pp) /dMdx F dp T and d z a(pd) / dMdx F dp T for each of the six data sets used 
in the analysis. For the purposes of the following discussion, we denote the mea- 
sured cross section in a given M, x F , pt bin as (Xj±<5crj ±A<7j, where % represents the 
data set, 5a,i is the statistical uncertainty and Ao"j is the systematic uncertainty. 

The results from the six data sets were first combined into three results corre- 
sponding to each mass setting. We used the weighted average 



where the weights were given by Ui = * , 2 . After the average cross section for 
each mass setting was calculated, the three mass settings were combined into a 
single result, again using equation 4.20. This procedure neglected the systematic 
uncertainties between data sets and mass settings, which were sufficiently small so 
that the average was not affected much by the omission. Also, it would not have 
been appropriate to use the systematic uncertainties in computing the averages, as 
most of the systematic uncertainty in the cross sections were (by design) common 
to the hydrogen and deuterium measurements. 




(4.20) 
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4.6.1 Systematic Uncertainties 

Computing the systematic uncertainty in the combined cross section was a fairly 
straightforward exercise in error propagation. The standard error propagation 
formula for an expression f(xi,x 2 , ■■■,xn) is given by 

( A /) 2 = E(^) 2 ( A ^) 2 - ( 4 - 21 ) 

From this formula, it can be shown that the systematic uncertainty in equation 
4.20 is given by 



^ 2 Ei^(A ^) 2 



(Aaf = 7 ( 4 -22) 



where the uii are the statistical weights used in the weighted average. One need 
only compute the systematic uncertainty in each data set (mass setting) A<7j. 

Recalling equation 4.6, the systematic uncertainty on a given cross section 
measurement will, by equation 4.21, be given by 

' ' + — + V + hr^ ( 4 -23) 



We must be careful here, however. The uncertainty in the luminosity AC, for 
instance, is dominated by a ±6.5% uncertainty in the beam normalization. This 
uncertainty is common to all data sets and mass settings. It would not be appro- 
priate to treat it as an independent 6.5% uncertainty in each measurement. We 
also assume that the uncertainties in detector efficiencies Ae are common to each 
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data set in a given mass setting, since events with similar kinematics will occupy 
similar regions of the spectrometer. We discuss how the systematic uncertainties 
were treated below. 



4.6.2 Uncorrelated Systematic Uncertainties 

The systematic uncertainties, not including the ±6.5% normalization uncertainty, 
fell into three categories: uncorrelated uncertainties, uncorrelated uncertainties 
which are common to a given mass setting, and correlated uncertainties which 
are common to a given mass setting. In general, we treated these systematic 
uncertainties separately. In the end they were combined in quadrature to estimate 
the total systematic uncertainty on the averaged cross sections. 

The only uncorrelated systematic uncertainty which was not shared between 
the data sets in a given mass setting was the uncertainty in the acceptance due 
to the statistical uncertainty in the Monte Carlo. The uncertainty in each mass 
setting due to this statistical uncertainty was therefore estimated to be, according 
to the prescription of equation 4.22 

2 I A a,stat\ 2 

(A^„) 2 = E -^ 2 } (4.24) 

where Aaf ,stat = —^—^i is the uncertainty in the i th measurement due to the 
statistical uncertainty in the acceptance Actj. 

The next step combined all of the uncorrelated uncertainties which were com- 
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mon to a given mass setting. These included the uncertainties in the acceptance 
in each mass setting, calculated above, and the uncertainties in the detector effi- 
ciencies. Again we use equation 4.22 

, A x2 _ Ei^ 2 {(^af + (Ae hodo ) 2 + (Ae chambcr ) 2 + (Ae trig ) 2 } 

l^°uncorrelatcdJ — 77=^ "2 l 4 -^ 

where the sum is now over all mass settings instead of all data sets in a given 
mass setting. 

The remaining source of systematic uncertainty was the uncertainty in the ac- 
ceptance due to uncertainties in the magnetic-field tweeks as measured in chapter 
3. Unlike the previous sources of uncertainties we have discussed, we expect this 
uncertainty to be at least partially correlated from mass setting to mass setting - 
the SM3 tweek was wholly determined by measurements from the high mass data. 

An estimation of the size of this uncertainty in any given mass setting was 
performed using a Monte Carlo based on the code described above, modified to 
propagate a given event through the spectrometer multiple times with different 
settings of the magnetic fields. In this modified Monte Carlo, all random processes, 
such as energy loss, multiple scattering and detector resolutions and efficiencies, 
were turned off. The result was a highly correlated sample of Monte Carlo events 
where any differences in reconstructed yields were due solely to acceptance effects. 
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4.6.3 Correlated Systematic Uncertainties 

The systematic uncertainty in the acceptance due to the uncertainties in the 
magnetic field in any given mass setting could then be estimated by 



Here, iV nominal represented the number of Monte Carlo events subjected to the 
nominal magnetic field tweeks in traversing the spectrometer which survived the 
standard analysis chain. Similarly, iV nomina i +< 5 were the number of Monte Carlo 
events subjected to a set of augmented magnetic fields in traversing the spectrom- 
eter, which were again subjected to the standard analysis chain. The field tweeks 
in this case were all augmented by the estimates of their uncertainties, 5 which 
should result in the maximum possible uncertainty. 

This results in a systematic uncertainty in the cross section for each mass 
setting given by 



To estimate the uncertainty in the total cross section, we depart from the pre- 
scription given by equation 4.22. We do this, because equation 4.22 is based on 

the assumption that the systematic uncertainties obey Gaussian statistics. This 

5 In the case of the low mass data, the SM12 and SM3 field tweeks were augmented by 
their systematic uncertainties while the SMO field tweek was reduced by its uncertainty. This 
resulted in the maximum possible systematic uncertainty in the acceptance for that spectrometer 
configuration. 




(4.26) 




(4.27) 
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is a good assumption for many of the sources of uncorrelated systematic uncer- 
tainties. Here it is not. Therefore we treat the systematic uncertainty Acrf ' tweek as 
a potentially "missing" or "excess" cross section for each mass setting, for which 
we wish to obtain an average estimate for the total cross section. We take 

EA a.tweek 

^"correlated — l 4 -^°J 

where once again the u>i are the statistical weights of the given mass setting. Fi- 
nally, we add the correlated and uncorrelated systematic uncertainties in quadra- 
ture to obtain an estimate for the total systematic uncertainty 



Act = sj (Acr uncorrelatcd ) 2 + (Acr corrclatcd ) 2 . (4.29) 
4.6.4 Integrating the Data 

After the data were combined using equation 4.20, we integrated the triple- 
differential cross sections over to obtain the double-differential cross sections 
d 2 a(pp)/dMdx F and d 2 a(pd) / dMdx F . The integrated cross sections were given 
by 

= f _ A PT „ + f(M, XF ; a) (4.30) 

where the sum is over all N p T bins whose widths are denoted Ap Tj j and /(M, xF; a) 
represents a correction explained below. The statistical and systematic uncertain- 
ties 6 on the integrated cross sections were calculated using the usual error prop- 
correlated and uncorrelated systematic uncertainties were first treated separately, then 
combined in quadrature to obtain the total systematic uncertainty. 
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agation in equation 4.23 - in other words, the uncertainties (multiplied by Ap T ) 
were added in quadrature. 

The correction in the integral f{M,xp;cx) accounted for the fact that the 
acceptance of the spectrometer could not reach the full kinematic limit in p T . For 
each bin in mass and xf, the triple-differential cross sections were integrated only 
over the range in p T where the acceptance was nonzero. Where the acceptance was 
zero fir all data sets, we made an estimate of the missing cross section f(M, xp] cx). 
The estimate was based on the empirical px distribution given in equation 4.12 
and the parameterization of given in equation 4.17, extrapolated from the 
largest nonzero Pt bin to infinity. 

The resulting corrections were well below 1% over much of the kinematic range 
covered. Only at larger masses did the correction to the cross sections approach 
5%. While the empirical form for the p T dependence of the cross section does 
a fairly good job of reproducing the data, it tends to underestimate the tail of 
the distribution. Given this, and the uncertainties in the parameterization of p^,, 
we estimate that the correction f(M,x F ;a) is uncertain by about ±15%. This 
uncertainty was added in quadrature to the systematic uncertainties in the final 
result. 



101 



5 RESULTS 

Fermilab Experiment 866 has measured continuum dimuon production cross sec- 
tions in 800-GeV pp and pd interactions. These data, which are presented below, 
represent the most extenive study of the differential pp cross to date, while the 
pd data offer better precision over a broader range of xp than previous mea- 
surements have achieved. In this chapter we tabulate and plot the double- and 
triple-differential cross sections, while deferring a discussion of the results until 
the next chapter. 

5.1 FNAL E866/NuSea Results - M 3 d 2 a/dMdx F 

The scaling form M 3 d 2 a/dMdxF for the hydrogen and deuterium cross sections 
are tabulated below in tables 5.1 and 5.2. The data are tabulated in mass and 
xf = ^jjj- bins. The mass and xf ranges are shown in the tables, as are the 
mean mass, xp and pp in each bin. The data in the tables are plotted in figures 
5.1 through 5.16, and compared to NLO calculations of the cross sections using 
different sets of parton distributions. 
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Table 5.1: Scaling form M 3 d 2 a/dMdxF ( in units of nb GeV 2 / nucleon ) for the hydrogen 
cross section. Statistical and point-to-point systematic uncertainties are shown separately. 
The normalization is subject to an uncertainty of ±6.5%. 



Mass [GeV] < M > < x F > < p T > M 3 d 2 a/dMdx F Stat. Error Syst. Error 



-0.05 <x F < 0.05 



4.20 - 


4.70 


4.64 


0.046 


0.65 


7.92E+01 


5.60E+01 


1.06E4 


-01 


4.70 - 


5.20 


5.01 


0.044 


0.36 


1.23E+01 


6.75E+00 


l.HE-f 


-00 


5.20 - 


5.70 


5.45 


0.039 


0.48 


3.01E+01 


1.81E+01 


3.84E4 


-00 


5.70 - 


6.20 


6.07 


0.036 


1.10 


4.02E+01 


1.47E+01 


3.06E+00 


6.20 - 


6.70 


6.46 


0.029 


0.93 


1.08E+01 


3.52E+00 


5.38E- 


-01 


6.70 - 


7.20 


6.89 


0.019 


0.68 


7.46E+00 


2.07E+00 


2.78E- 


-01 


7.20 - 


7.70 


7.54 


0.024 


0.84 


6.17E+00 


2.56E+00 


4.24E- 


-01 


7.70 - 


8.20 


8.04 


0.028 


0.94 


5.50E+00 


9.03E-01 


1.91E- 


-01 


8.20 - 


8.70 


8.46 


0.016 


0.85 


3.70E+00 


4.69E-01 


1.07E- 


-01 


10.85 - 


11.85 


10.86 


0.047 


1.43 


8.02E+00 


8.02E+00 


8.85E- 


-01 



0.05 <x F < 0.10 



4.20 - 4.70 4.46 0.083 0.77 2.67E+01 2.76E+00 7.73E-01 

4.70 - 5.20 4.90 0.080 0.76 1.54E+01 2.28E+00 7.09E-01 



Table 5.1: (continued) 



Mass [GeV] < M > < x F > < p T > M 3 d 2 a/dMdx F Stat. Error Syst. Error 



5.20 - 


5.70 


5.43 


0.083 


0.86 


1.06E+01 


2.83E+00 


1.03E4 


-00 


5.70 - 


6.20 


5.95 


0.083 


0.94 


1.24E+01 


2.59E+00 


7.74E- 


-01 


6.20 - 


6.70 


6.48 


0.078 


0.98 


6.37E+00 


1.82E+00 


4.82E- 


-01 


6.70 - 


7.20 


6.93 


0.076 


0.86 


6.13E+00 


2.08E+00 


2.29E- 


-01 


7.20 - 


7.70 


7.50 


0.080 


0.98 


6.60E+00 


1.95E+00 


4.68E- 


-01 


7.70 - 


8.20 


7.95 


0.078 


1.04 


4.17E+00 


5.43E-01 


2.23E- 


-01 


8.20 - 


8.70 


8.44 


0.077 


1.03 


3.71E+00 


4.83E-01 


1.75E- 


-01 


10.85 - 


11.85 


11.17 


0.081 


1.15 


3.46E+00 


1.51E+00 


7.66E- 


-02 



0.10 <x F < 0.15 



4.20 - 


4.70 


4.43 


0.129 


0.91 


2.34E+01 


1.38E+00 


5.31E- 


-01 


4.70 - 


5.20 


4.92 


0.127 


0.90 


1.59E+01 


1.26E+00 


4.18E- 


-01 


5.20 - 


5.70 


5.41 


0.130 


0.92 


1.54E+01 


1.58E+00 


5.82E- 


-01 


5.70 - 


6.20 


5.95 


0.128 


0.91 


1.26E+01 


1.45E+00 


3.60E- 


-01 


6.20 - 


6.70 


6.44 


0.130 


1.04 


9.98E+00 


1.22E+00 


2.94E- 


-01 


6.70 - 


7.20 


6.97 


0.130 


0.99 


6.86E+00 


1.05E+00 


3.05E- 


-01 


7.20 - 


7.70 


7.46 


0.128 


1.10 


5.06E+00 


4.80E-01 


1.84E- 


-01 


7.70 - 


8.20 


7.95 


0.127 


1.10 


4.56E+00 


3.67E-01 


1.47E- 


-01 


8.20 - 


8.70 


8.44 


0.127 


1.12 


3.11E+00 


2.61E-01 


9.49E- 


-02 



Table 5.1: (continued) 



Mass [GeV] < M > < x F > < p T > M 3 d 2 a/dMdx F Stat. Error Syst. Error 



10.85 - 11.85 11.22 0.129 1.10 1.36E+00 3.42E-01 5.81E-02 













0.15 < x F 


< 0.20 










4.20 - 


4.70 


4.43 





177 


1.00 


1.92E+01 


9.66E- 


-01 


1.21E4 


-00 


4.70 - 


5.20 


4.94 





177 


0.96 


1.45E+01 


7.20E- 


-01 


2.55E- 


-01 


5.20 - 


5.70 


5.43 





177 


1.02 


1.20E+01 


9.89E- 


-01 


2.71E- 


-01 


5.70 - 


6.20 


5.91 





178 


1.02 


8.96E+00 


7.88E- 


-01 


1.99E- 


-01 


6.20 - 


6.70 


6.45 





179 


1.02 


8.04E+00 


8.68E- 


-01 


3.04E- 


-01 


6.70 - 


7.20 


6.96 





179 


1.06 


6.22E+00 


4.59E- 


-01 


1.98E- 


-01 


7.20 - 


7.70 


7.46 





177 


1.12 


4.82E+00 


2.98E- 


-01 


1.18E- 


-01 


7.70 - 


8.20 


7.94 





176 


1.12 


3.61E+00 


2.45E- 


-01 


9.41E- 


-02 


8.20 - 


8.70 


8.44 





176 


1.19 


3.00E+00 


1.98E- 


-01 


6.28E- 


-02 


10.85 - 


11.85 


11.18 





178 


1.09 


1.08E+00 


2.41E- 


-01 


7.37E- 


-02 


11.85 - 


12.85 


12.23 





177 


1.37 


1.25E+00 


5.62E- 


-01 


1.54E- 


-01 


12.85 - 


14.85 


13.31 





185 


1.24 


1.17E+01 


1.03E+ 


-01 


3.13E- 


-01 



Table 5.1: (continued) 



Mass [GeV] < M > < x F > < p T > M 3 d 2 a/dMdx F Stat. Error Syst. Error 



0.20 <x F < 0.25 



4 


20 - 


4.70 


4.43 


0.226 


1 


02 


1.68E+01 


6 


71E-01 


2 


97E- 


-01 


4 


70 - 


5.20 


4.93 


0.225 


1 


02 


1.26E+01 


6 


54E-01 


2 


35E- 


-01 


5 


20 - 


5.70 


5.43 


0.225 


1 


00 


9.38E+00 


5 


79E-01 


1 


77E- 


-01 


5 


70 - 


6.20 


5.95 


0.228 


1 


05 


8.45E+00 


5 


19E-01 


2 


48E- 


-01 


6 


20 - 


6.70 


6.46 


0.228 





99 


7.22E+00 


3 


73E-01 


1 


88E- 


-01 


6 


70 - 


7.20 


6.95 


0.227 


1 


09 


5.18E+00 


2 


72E-01 


1 


09E- 


-01 


7 


20 - 


7.70 


7.45 


0.227 


1 


12 


3.91E+00 


2 


28E-01 


1 


13E- 


-01 


7 


70 - 


8.20 


7.94 


0.225 


1 


16 


3.73E+00 


1 


97E-01 


8 


98E- 


-02 


8 


20 - 


8.70 


8.44 


0.224 


1 


19 


2.76E+00 


1 


89E-01 


6 


10E- 


-02 


10 


85 - 


11.85 


11.22 


0.231 


1 


45 


7.48E-01 


1 


34E-01 


4 


23E- 


-02 


11 


85 - 


12.85 


12.33 


0.226 


1 


33 


2.02E+00 


7 


89E-01 


6 


25E- 


-01 


12 


85 - 


14.85 


13.22 


0.226 


1 


04 


2.36E+00 


1 


82E+00 


3 


90E- 


-01 



0.25 < x F < 0.30 

4.20 - 4.70 4.44 0.274 1.04 1.23E+01 5.17E-01 2.79E-01 

4.70 - 5.20 4.95 0.276 1.04 1.01E+01 4.91E-01 1.74E-01 



Table 5.1: (continued) 



Mass [GeV] < M > < x F > < p T > M 3 d 2 a/dMdx F Stat. Error Syst. Error 



5.20 - 


5.70 


5.44 


0.276 


1 


04 


7.25E+00 


4.24E- 


-01 


2.03E- 


-01 


5.70 - 


6.20 


5.94 


0.277 


1 


01 


7.53E+00 


3.04E- 


-01 


1.93E- 


-01 


6.20 - 


6.70 


6.45 


0.276 


1 


04 


5.64E+00 


2.62E- 


-01 


1.49E- 


-01 


6.70 - 


7.20 


6.94 


0.276 


1 


09 


4.46E+00 


2.09E- 


-01 


8.33E- 


-02 


7.20 - 


7.70 


7.43 


0.275 


1 


12 


3.31E+00 


1.77E- 


-01 


7.80E- 


-02 


7.70 - 


8.20 


7.94 


0.276 


1 


17 


2.93E+00 


1.67E- 


-01 


5.92E- 


-02 


8.20 - 


8.70 


8.44 


0.275 


1 


26 


2.69E+00 


1.53E- 


-01 


5.86E- 


-02 


10.85 - 


11.85 


11.23 


0.278 


1 


19 


6.93E-01 


1.09E- 


-01 


5.55E- 


-02 


11.85 - 


12.85 


12.29 


0.282 


1 


31 


4.83E-01 


1.71E- 


-01 


5.21E- 


-02 


12.85 - 


14.85 


13.31 


0.278 


1 


57 


4.68E-01 


3.33E- 


-01 


5.77E- 


-02 



0.30 < x F < 0.35 



4.20 - 


4.70 


4.45 


0.324 


1.03 


1.00E+01 


4.63E- 


-01 


2.09E- 


-01 


4.70 - 


5.20 


4.95 


0.325 


1.01 


8.85E+00 


4.15E- 


-01 


2.02E- 


-01 


5.20 - 


5.70 


5.45 


0.326 


1.02 


6.91E+00 


3.03E- 


-01 


1.89E- 


-01 


5.70 - 


6.20 


5.95 


0.325 


1.02 


5.86E+00 


2.31E- 


-01 


1.21E- 


-01 


6.20 - 


6.70 


6.44 


0.324 


1.03 


4.72E+00 


1.98E- 


-01 


8.29E- 


-02 


6.70 - 


7.20 


6.94 


0.325 


1.11 


4.03E+00 


1.78E- 


-01 


1.01E- 


-01 


7.20 - 


7.70 


7.45 


0.325 


1.15 


3.37E+00 


1.50E- 


-01 


7.64E- 


-02 



Table 5.1: (continued) 



Mass [GeV] < M > < x F > < p T > M 3 d 2 a/dMdx F Stat. Error Syst. Error 



7.70 - 


8.20 


7.93 


0.325 


1.17 


2.51E+00 


1.41E- 


-01 


4.85E- 


-02 


8.20 - 


8.70 


8.43 


0.325 


1.19 


2.22E+00 


1.33E- 


-01 


4.33E- 


-02 


10.85 - 


11.85 


11.32 


0.326 


1.27 


6.56E-01 


8.22E- 


-02 


5.27E- 


-02 


11.85 - 


12.85 


12.26 


0.328 


1.12 


9.65E-01 


4.94E- 


-01 


4.78E- 


-02 


12.85 - 


14.85 


13.30 


0.329 


1.31 


4.92E-01 


1.78E- 


-01 


4.26E- 


-02 



0.35 <x F < 0.40 



4.20 - 


4.70 


4.46 





374 


1 


02 


7.82E+00 


4.14E- 


-01 


2.43E- 


-01 


4.70 - 


5.20 


4.95 





376 


1 


04 


6.54E+00 


3.20E- 


-01 


2.26E- 


-01 


5.20 - 


5.70 


5.45 





375 


1 


02 


5.87E+00 


3.87E- 


-01 


4.67E- 


-01 


5.70 - 


6.20 


5.96 





375 





99 


4.78E+00 


2.64E- 


-01 


1.63E- 


-01 


6.20 - 


6.70 


6.45 





375 


1 


05 


3.86E+00 


1.62E- 


-01 


6.75E- 


-02 


6.70 - 


7.20 


6.94 





375 


1 


10 


3.30E+00 


1.48E- 


-01 


9.22E- 


-02 


7.20 - 


7.70 


7.44 





376 


1 


17 


2.67E+00 


1.21E- 


-01 


8.27E- 


-02 


7.70 - 


8.20 


7.95 





375 


1 


13 


2.21E+00 


1.24E- 


-01 


6.03E- 


-02 


8.20 - 


8.70 


8.44 





374 


1 


16 


1.90E+00 


1.09E- 


-01 


4.06E- 


-02 


10.85 - 


11.85 


11.32 





375 


1 


15 


4.27E-01 


7.23E- 


-02 


5.12E- 


-03 


11.85 - 


12.85 


12.26 





374 


1 


29 


3.45E-01 


6.23E- 


-02 


1.31E- 


-02 



Table 5.1: (continued) 



Mass [GeV] 


< M > 


< Xp > 


< Pt > 


M 3 d 2 a/dMdx F 


Stat. Error 


Syst. Error 


12.85 - 14.85 


13.43 


0.377 


1.31 


3.77E-01 


1.91E-01 


1.30E-02 





0.40 <x F < 0.45 



4 


20 - 


4.70 


4.46 


0.424 





99 


6.45E+00 


3 


68E- 


-01 


2 


25E- 


-01 


4 


70 - 


5.20 


4.97 


0.423 


1 


01 


5.32E+00 


2 


50E- 


-01 


1 


52E- 


-01 


5 


20 - 


5.70 


5.46 


0.424 





97 


4.13E+00 


1 


93E- 


-01 


1 


28E- 


-01 


5 


70 - 


6.20 


5.95 


0.425 


1 


01 


3.76E+00 


1 


53E- 


-01 


8 


14E- 


-02 


6 


20 - 


6.70 


6.45 


0.424 


1 


02 


3.20E+00 


1 


24E- 


-01 


8 


05E- 


-02 


6 


70 - 


7.20 


6.95 


0.424 


1 


12 


2.50E+00 


1 


16E- 


-01 


6 


28E- 


-02 


7 


20 - 


7.70 


7.45 


0.424 


1 


12 


2.04E+00 


1 


02E- 


-01 


5 


64E- 


-02 


7 


70 - 


8.20 


7.94 


0.424 


1 


15 


1.64E+00 


9 


63E- 


-02 


5 


20E- 


-02 


8 


20 - 


8.70 


8.44 


0.424 


1 


20 


1.50E+00 


9 


25E- 


-02 


5 


29E- 


-02 


10 


85 - 


11.85 


11.28 


0.423 


1 


36 


3.63E-01 


4 


95E- 


-02 


8 


00E- 


-03 


11 


85 - 


12.85 


12.24 


0.424 


1 


29 


2.46E-01 


5 


32E- 


-02 


2 


24E- 


-02 


12 


85 - 


14.85 


13.46 


0.429 


1 


41 


6.71E-01 


5 


45E- 


-01 


2 


18E- 


-02 



Table 5.1: (continued) 



Mass [GeV] < M > < x F > < p T > M 3 d 2 a/dMdx F Stat. Error Syst. Error 



0.45 <x F < 0.50 



4 


20 - 


4.70 


4.46 


0.473 





95 


4 


20E+00 


2 


80E- 


-01 


1 


27E- 


-01 


4 


70 - 


5.20 


4.95 


0.475 





99 


3 


55E+00 


2 


12E- 


-01 


9 


38E- 


-02 


5 


20 - 


5.70 


5.46 


0.475 





99 


3 


22E+00 


1 


54E- 


-01 


6 


59E- 


-02 


5 


70 - 


6.20 


5.95 


0.474 





97 


2 


54E+00 


1 


20E- 


-01 


5 


07E- 


-02 


6 


20 - 


6.70 


6.44 


0.474 


1 


03 


1 


89E+00 


9 


63E- 


-02 


3 


76E- 


-02 


6 


70 - 


7.20 


6.94 


0.474 


1 


09 


1 


88E+00 


9 


42E- 


-02 


5 


12E- 


-02 


7 


20 - 


7.70 


7.43 


0.474 


1 


06 


1 


50E+00 


8 


66E- 


-02 


5 


23E- 


-02 


7 


70 - 


8.20 


7.94 


0.473 


1 


08 


1 


21E+00 


8 


03E- 


-02 


2 


65E- 


-02 


8 


20 - 


8.70 


8.44 


0.474 


1 


24 


8 


83E-01 


7 


12E- 


-02 


2 


05E- 


-02 


10 


85 - 


11.85 


11.29 


0.472 


1 


24 


2 


50E-01 


4 


24E- 


-02 


8 


17E- 


-03 


11 


85 - 


12.85 


12.33 


0.471 


1 


18 


1 


74E-01 


4 


25E- 


-02 


6 


59E- 


-03 


12 


85 - 


14.85 


13.49 


0.474 





90 


9 


36E-02 


3 


83E- 


-02 


8 


91E- 


-03 



0.50 <x F < 0.55 

4.20 - 4.70 4.48 0.524 1.04 2.81E+00 2.42E-01 1.13E-01 

4.70 - 5.20 4.97 0.524 1.01 2.56E+00 1.76E-01 1.14E-01 



Table 5.1: (continued) 



Mass [GeV] < M > < x F > < p T > M 3 d 2 a/dMdx F Stat. Error Syst. Error 



5.20 - 


5.70 


5.46 





523 


1 


01 


2 


28E+00 


1.23E-01 


6 


86E- 


-02 


5.70 - 


6.20 


5.94 





524 


1 


00 


1 


83E+00 


9.65E-02 


6 


27E- 


-02 


6.20 - 


6.70 


6.45 





525 


1 


03 


1 


51E+00 


7.93E-02 


3 


54E- 


-02 


6.70 - 


7.20 


6.93 





523 


1 


07 


1 


27E+00 


7.58E-02 


3 


37E- 


-02 


7.20 - 


7.70 


7.44 





523 


1 


07 


9 


84E-01 


7.11E-02 


2 


33E- 


-02 


7.70 - 


8.20 


7.94 





523 


1 


14 


9 


85E-01 


7.06E-02 


2 


27E- 


-02 


8.20 - 


8.70 


8.45 





523 


1 


18 


8 


52E-01 


6.21E-02 


1 


93E- 


-02 


10.85 - 


11.85 


11.23 





527 


1 


14 


1 


84E-01 


3.27E-02 


1 


28E- 


-03 


11.85 - 


12.85 


12.16 





525 


1 


24 


7 


93E-02 


2.51E-02 


2 


80E- 


-03 


12.85 - 


14.85 


13.46 





524 


1 


12 


1 


69E-01 


7.29E-02 


1 


30E- 


-02 


14.85 - 


16.85 


14.87 





531 





77 


1 


80E+00 


1.80E+00 


6 


50E- 


-01 



0.55 <x F < 0.60 



4.20 - 


4.70 


4.46 


0.571 


0.97 


1.45E+00 


1.61E- 


-01 


6.48E- 


-02 


4.70 - 


5.20 


4.96 


0.575 


0.97 


1.57E+00 


1.20E- 


-01 


6.31E- 


-02 


5.20 - 


5.70 


5.45 


0.573 


0.91 


1.44E+00 


9.82E- 


-02 


3.99E- 


-02 


5.70 - 


6.20 


5.96 


0.573 


0.89 


1.16E+00 


7.25E- 


-02 


3.03E- 


-02 


6.20 - 


6.70 


6.44 


0.573 


1.00 


1.15E+00 


7.83E- 


-02 


2.50E- 


-02 


6.70 - 


7.20 


6.94 


0.572 


1.06 


8.69E-01 


6.05E- 


-02 


1.99E- 


-02 



Table 5.1: (continued) 



Mass [GeV] < M > < x F > < p T > M 3 d 2 a/dMdx F Stat. Error Syst. Error 



7.20 - 


7.70 


7.42 


0.573 


1.09 


7.86E- 


-01 


5.48E- 


-02 


1.84E- 


-02 


7.70 - 


8.20 


7.94 


0.572 


1.16 


6.01E- 


-01 


5.13E- 


-02 


2.25E- 


-02 


8.20 - 


8.70 


8.44 


0.573 


1.03 


5.52E- 


-01 


5.05E- 


-02 


1.82E- 


-02 


10.85 - 


11.85 


11.26 


0.570 


1.16 


1.16E- 


-01 


2.27E- 


-02 


4.33E- 


-03 


11.85 - 


12.85 


12.21 


0.576 


1.03 


9.16E- 


-02 


2.83E- 


-02 


6.40E- 


-03 


12.85 - 


14.85 


13.28 


0.570 


1.14 


2.50E- 


-01 


2.02E- 


-01 


5.16E- 


-03 



0.60 <x F < 0.65 



4.20 - 


4.70 


4.49 


0.622 





95 


1.41E4 


-00 


1.74E- 


-01 


8 


92E- 


-02 


4.70 - 


5.20 


4.96 


0.623 





94 


9.66E- 


-01 


1.14E- 


-01 


5 


80E- 


-02 


5.20 - 


5.70 


5.46 


0.622 





93 


8.05E- 


-01 


8.03E- 


-02 


3 


26E- 


-02 


5.70 - 


6.20 


5.96 


0.623 





95 


7.15E- 


-01 


5.35E- 


-02 


2 


29E- 


-02 


6.20 - 


6.70 


6.43 


0.623 


1 


02 


6.39E- 


-01 


5.76E- 


-02 


2 


27E- 


-02 


6.70 - 


7.20 


6.93 


0.624 


1 


02 


6.79E- 


-01 


5.35E- 


-02 


2 


52E- 


-02 


7.20 - 


7.70 


7.44 


0.623 


1 


02 


5.44E- 


-01 


6.51E- 


-02 


3 


05E- 


-02 


7.70 - 


8.20 


7.93 


0.624 


1 


03 


3.85E- 


-01 


4.04E- 


-02 


1 


19E- 


-02 


8.20 - 


8.70 


8.45 


0.624 


1 


08 


3.08E- 


-01 


3.93E- 


-02 


1 


04E- 


-02 


10.85 - 


11.85 


11.24 


0.625 


1 


13 


1.09E- 


-01 


3.33E- 


-02 


6 


55E- 


-03 


11.85 - 


12.85 


12.18 


0.620 


1 


34 


6.14E- 


-02 


2.93E- 


-02 


5 


11E- 


-03 



Table 5.1: (continued) 



Mass [GeV] < M > < x F > < p T > M 3 d 2 a/dMdx F Stat. Error Syst. Error 



12.85 - 14.85 13.16 0.609 0.67 2.48E-02 1.13E-02 2.18E-03 











0.65 < x F 


< 0.70 












4.20 - 


4.70 


4.47 


0.673 


0.97 


6.90E- 


-01 


1.18E- 


-01 


6.38E- 


-02 


4.70 - 


5.20 


4.97 


0.670 


0.90 


5.91E- 


-01 


7.45E- 


-02 


3.13E- 


-02 


5.20 - 


5.70 


5.45 


0.670 


0.88 


6.39E- 


-01 


7.28E- 


-02 


3.90E- 


-02 


5.70 - 


6.20 


5.95 


0.673 


0.89 


4.48E- 


-01 


6.71E- 


-02 


1.74E- 


-02 


6.20 - 


6.70 


6.46 


0.671 


0.92 


3.80E- 


-01 


5.04E- 


-02 


1.68E- 


-02 


6.70 - 


7.20 


6.94 


0.669 


1.03 


2.48E- 


-01 


3.91E- 


-02 


1.40E- 


-02 


7.20 - 


7.70 


7.43 


0.673 


1.06 


1.97E- 


-01 


4.09E- 


-02 


1.65E- 


-02 


7.70 - 


8.20 


7.95 


0.670 


1.08 


2.21E- 


-01 


3.71E- 


-02 


1.27E- 


-02 


8.20 - 


8.70 


8.43 


0.669 


1.07 


1.04E- 


-01 


2.42E- 


-02 


5.03E- 


-03 


10.85 - 


11.85 


11.30 


0.674 


1.22 


5.25E- 


-02 


1.52E- 


-02 


2.79E- 


-03 


11.85 - 


12.85 


12.42 


0.679 


1.27 


3.87E- 


-02 


2.09E- 


-02 


3.70E- 


-03 


12.85 - 


14.85 


13.77 


0.677 


0.83 


9.03E- 


-02 


4.41E- 


-02 


3.04E- 


-03 



Table 5.1: (continued) 



Mass [GeV] < M > < x F > < p T > M 3 d 2 a/dMdx F Stat. Error Syst. Error 











0.70 < x F 


< 0.75 












4.20 - 


4.70 


4.53 


0.720 


0.80 


5.31E- 


-01 


1.26E- 


-01 


4.75E- 


-02 


4.70 - 


5.20 


5.01 


0.720 


0.90 


4.21E- 


-01 


8.97E- 


-02 


2.60E- 


-02 


5.20 - 


5.70 


5.45 


0.719 


0.94 


2.44E- 


-01 


3.87E- 


-02 


1.16E- 


-02 


5.70 - 


6.20 


5.94 


0.723 


1.04 


2.68E- 


-01 


4.36E- 


-02 


2.31E- 


-02 


6.20 - 


6.70 


6.45 


0.720 


1.00 


2.04E- 


-01 


4.72E- 


-02 


1.62E- 


-02 


6.70 - 


7.20 


6.93 


0.720 


0.86 


1.62E- 


-01 


2.65E- 


-02 


1.24E- 


-02 


7.20 - 


7.70 


7.41 


0.726 


1.00 


1.78E- 


-01 


3.30E- 


-02 


1.06E- 


-02 


7.70 - 


8.20 


7.95 


0.719 


1.07 


3.54E- 


-01 


1.61E- 


-01 


1.21E- 


-02 


8.20 - 


8.70 


8.42 


0.717 


1.13 


1.05E- 


-01 


2.43E- 


-02 


1.05E- 


-02 


10.85 - 


11.85 


11.28 


0.728 


1.30 


3.61E- 


-02 


2.34E- 


-02 


4.30E- 


-04 


12.85 - 


14.85 


13.14 


0.723 


1.33 


3.83E- 


-01 


2.72E- 


-01 


4.21E- 


-02 



0.75 <x F < 0.80 

4.20 - 4.70 4.48 0.770 1.01 3.20E-01 1.11E-01 5.33E-02 

4.70 - 5.20 4.96 0.767 0.84 4.98E-01 2.48E-01 3.98E-02 

5.20 - 5.70 5.45 0.768 0.65 1.28E-01 2.67E-02 9.63E-03 



Table 5.1: (continued) 



Mass [GeV] < M > < x F > < p T > M 3 d 2 a/dMdx F Stat. Error Syst. Error 



5.70 - 


6.20 


5.95 


0.771 


0.96 


1.76E- 


-01 


4.18E- 


-02 


1.37E- 


-02 


6.20 - 


6.70 


6.43 


0.769 


0.82 


9.89E- 


-02 


2.22E- 


-02 


5.33E- 


-03 


6.70 - 


7.20 


6.94 


0.767 


0.83 


8.01E- 


-02 


2.43E- 


-02 


5.63E- 


-03 


7.20 - 


7.70 


7.45 


0.773 


0.97 


1.09E- 


-01 


4.14E- 


-02 


1.91E- 


-02 


7.70 - 


8.20 


7.93 


0.772 


1.42 


1.04E- 


-01 


3.55E- 


-02 


9.37E- 


-03 


8.20 - 


8.70 


8.47 


0.776 


0.80 


7.31E- 


-02 


2.44E- 


-02 


6.20E- 


-03 


10.85 - 


11.85 


11.39 


0.765 


1.52 


6.65E- 


-02 


4.97E- 


-02 


5.74E- 


-03 



Table 5.2: Scaling form M 3 d 2 a/dMdxF ( in units of nb GeV 2 / nucleon ) for the deuterium 
cross section. Statistical and point-to-point systematic uncertainties are shown separately. 
The normalization is subject to an uncertainty of ±6.5%. 



Mass [GeV] < M > < x F > < p T > M 3 d 2 a/dMdx F Stat. Error Syst. Error 











-0.05 < x F 


< 0.05 






4.70 - 


5.20 


4.94 


0.044 


0.96 


4.50E+01 


2.77E+01 


5.70E-01 


5.20 - 


5.70 


5.40 


0.039 


0.55 


1.81E+01 


8.30E+00 


2.00E+00 


5.70 - 


6.20 


6.10 


0.037 


1.04 


1.21E+01 


4.70E+00 


1.46E+00 


6.20 - 


6.70 


6.48 


0.022 


1.03 


1.25E+01 


2.75E+00 


9.14E-01 


6.70 - 


7.20 


6.95 


0.025 


0.90 


9.76E+00 


1.73E+00 


4.41E-01 


7.20 - 


7.70 


7.49 


0.027 


0.97 


7.15E+00 


2.47E+00 


1.12E+00 


7.70 - 


8.20 


8.01 


0.024 


0.99 


5.62E+00 


6.72E-01 


2.47E-01 


8.20 - 


8.70 


8.45 


0.017 


0.97 


4.32E+00 


4.32E-01 


2.15E-01 


10.85 - 


11.85 


11.12 


0.012 


0.44 


8.71E-01 


3.58E-01 


1.75E-02 



0.05 <x F < 0.10 

4.20 - 4.70 4.48 0.084 0.78 2.78E+01 2.08E+00 7.52E-01 

4.70 - 5.20 4.92 0.081 0.83 2.18E+01 1.73E+00 7.37E-01 

5.20 - 5.70 5.43 0.084 0.77 1.73E+01 2.55E+00 1.11E+00 



Table 5.2: (continued) 



Mass [GeV] < M > < x F > < p T > M 3 d 2 a/dMdx F Stat. Error Syst. Error 



5.70 - 


6.20 


5.97 


0.080 


0.98 


1.30E+01 


1.81E+00 


6.84E- 


-01 


6.20 - 


6.70 


6.46 


0.079 


0.94 


7.82E+00 


1.18E+00 


4.34E- 


-01 


6.70 - 


7.20 


6.90 


0.080 


0.90 


6.23E+00 


1.32E+00 


4.38E- 


-01 


7.20 - 


7.70 


7.50 


0.082 


0.91 


6.18E+00 


1.21E+00 


5.15E- 


-01 


7.70 - 


8.20 


7.97 


0.078 


1.03 


4.18E+00 


3.05E-01 


2.02E- 


-01 


8.20 - 


8.70 


8.45 


0.076 


1.03 


3.58E+00 


2.51E-01 


1.58E- 


-01 


10.85 - 


11.85 


11.21 


0.085 


0.98 


1.84E+00 


6.56E-01 


7.84E- 


-03 



0.10 <x F < 0.15 



4.20 - 


4.70 


4.44 





130 


0.91 


2.64E+01 


9.63E- 


-01 


4.98E- 


-01 


4.70 - 


5.20 


4.93 





128 


0.86 


1.72E+01 


7.96E- 


-01 


3.55E- 


-01 


5.20 - 


5.70 


5.42 





129 


0.92 


1.45E+01 


1.01E+00 


4.16E- 


-01 


5.70 - 


6.20 


5.93 





130 


0.98 


1.02E+01 


7.86E- 


-01 


2.91E- 


-01 


6.20 - 


6.70 


6.45 





129 


0.93 


8.98E+00 


8.55E- 


-01 


2.84E- 


-01 


6.70 - 


7.20 


6.97 





130 


1.04 


7.18E+00 


6.49E- 


-01 


2.73E- 


-01 


7.20 - 


7.70 


7.49 





128 


1.07 


5.91E+00 


4.46E- 


-01 


1.97E- 


-01 


7.70 - 


8.20 


7.95 





129 


1.13 


4.92E+00 


2.49E- 


-01 


1.82E- 


-01 


8.20 - 


8.70 


8.45 





127 


1.13 


3.86E+00 


1.87E- 


-01 


1.02E- 


-01 


10.85 - 


11.85 


11.20 





126 


1.17 


1.44E+00 


3.19E- 


-01 


1.55E- 


-01 



Table 5.2: (continued) 



Mass [GeV] < M > < x F > < p T > M 3 d 2 a/dMdx F Stat. Error Syst. Error 



11.85 - 12.85 12.03 0.135 1.72 2.12E+00 1.04E+00 9.62E-02 













0.15 < x F 


< 0.20 










4.20 - 


4.70 


4.43 





177 


1.03 


2.17E+01 


6.01E- 


-01 


4.30E- 


-01 


4.70 - 


5.20 


4.93 





177 


1.00 


1.68E+01 


5.27E- 


-01 


3.45E- 


-01 


5.20 - 


5.70 


5.43 





177 


0.97 


1.17E+01 


6.14E- 


-01 


2.86E- 


-01 


5.70 - 


6.20 


5.92 





178 


1.00 


1.05E+01 


5.32E- 


-01 


2.30E- 


-01 


6.20 - 


6.70 


6.46 





179 


1.08 


7.81E+00 


4.77E- 


-01 


2.93E- 


-01 


6.70 - 


7.20 


6.98 





179 


1.08 


6.90E+00 


3.01E- 


-01 


2.73E- 


-01 


7.20 - 


7.70 


7.44 





176 


1.11 


5.35E+00 


2.25E- 


-01 


1.52E- 


-01 


7.70 - 


8.20 


7.95 





176 


1.15 


4.46E+00 


1.73E- 


-01 


1.29E- 


-01 


8.20 - 


8.70 


8.43 





176 


1.16 


3.63E+00 


1.41E- 


-01 


7.58E- 


-02 


10.85 - 


11.85 


11.24 





178 


1.11 


1.17E+00 


2.06E- 


-01 


5.42E- 


-02 


11.85 - 


12.85 


12.06 





166 


0.91 


5.16E-01 


1.91E- 


-01 


5.58E- 


-02 


12.85 - 


14.85 


13.13 





195 


1.25 


1.82E+00 


1.79E4 


-00 


5.40E- 


-02 



Table 5.2: (continued) 



Mass [GeV] < M > < x F > < p T > M 3 d 2 a/dMdx F Stat. Error Syst. Error 



0.20 <x F < 0.25 



4 


20 - 


4.70 


4.44 


0.225 


1 


04 


1.82E+01 


4 


23E- 


-01 


3 


32E- 


-01 


4 


70 - 


5.20 


4.93 


0.225 


1 


01 


1.35E+01 


4 


00E- 


-01 


2 


44E- 


-01 


5 


20 - 


5.70 


5.43 


0.226 


1 


00 


1.03E+01 


3 


93E- 


-01 


2 


06E- 


-01 


5 


70 - 


6.20 


5.94 


0.227 


1 


03 


8.80E+00 


3 


40E- 


-01 


2 


30E- 


-01 


6 


20 - 


6.70 


6.46 


0.229 


1 


05 


7.77E+00 


2 


64E- 


-01 


1 


91E- 


-01 


6 


70 - 


7.20 


6.95 


0.227 


1 


09 


5.87E+00 


1 


87E- 


-01 


1 


15E- 


-01 


7 


20 - 


7.70 


7.43 


0.226 


1 


11 


4.73E+00 


1 


55E- 


-01 


1 


16E- 


-01 


7 


70 - 


8.20 


7.94 


0.226 


1 


18 


4.06E+00 


1 


35E- 


-01 


1 


02E- 


-01 


8 


20 - 


8.70 


8.45 


0.225 


1 


17 


3.35E+00 


1 


24E- 


-01 


6 


75E- 


-02 


10 


85 - 


11.85 


11.17 


0.227 


1 


25 


7.07E-01 


7 


43E- 


-02 


1 


97E- 


-02 


11 


85 - 


12.85 


12.31 


0.226 


1 


23 


5.63E-01 


1 


38E- 


-01 


4 


95E- 


-02 


12 


85 - 


14.85 


13.19 


0.229 





85 


3.61E-01 


1 


48E- 


-01 


2 


88E- 


-02 



0.25 < x F < 0.30 

4.20 - 4.70 4.44 0.275 1.01 1.46E+01 3.45E-01 2.39E-01 

4.70 - 5.20 4.94 0.275 1.03 1.15E+01 3.20E-01 1.97E-01 



Table 5.2: (continued) 



Mass [GeV] < M > < x F > < p T > M 3 d 2 a/dMdx F Stat. Error Syst. Error 



5.20 - 


5.70 


5.45 


0.276 


1 


02 


9.32E+00 


2.96E- 


-01 


2.30E- 


-01 


5.70 - 


6.20 


5.95 


0.277 


1 


01 


7.72E+00 


1.98E- 


-01 


2.19E- 


-01 


6.20 - 


6.70 


6.45 


0.276 


1 


02 


6.47E+00 


1.72E- 


-01 


1.68E- 


-01 


6.70 - 


7.20 


6.95 


0.276 


1 


08 


5.15E+00 


1.39E- 


-01 


8.95E- 


-02 


7.20 - 


7.70 


7.44 


0.275 


1 


15 


4.11E+00 


1.21E- 


-01 


9.33E- 


-02 


7.70 - 


8.20 


7.94 


0.276 


1 


19 


3.29E+00 


1.08E- 


-01 


6.60E- 


-02 


8.20 - 


8.70 


8.45 


0.275 


1 


20 


2.85E+00 


1.01E- 


-01 


5.41E- 


-02 


10.85 - 


11.85 


11.23 


0.275 


1 


19 


8.09E-01 


6.51E- 


-02 


2.10E- 


-02 


11.85 - 


12.85 


12.23 


0.278 


1 


04 


3.41E-01 


8.92E- 


-02 


2.38E- 


-02 


12.85 - 


14.85 


13.42 


0.270 


1 


29 


4.11E-01 


2.70E- 


-01 


2.52E- 


-02 



0.30 < x F < 0.35 



4.20 - 


4.70 


4.44 


0.323 


1.04 


1.13E+01 


3.26E- 


-01 


2.56E- 


-01 


4.70 - 


5.20 


4.95 


0.325 


1.05 


8.45E+00 


2.57E- 


-01 


1.86E- 


-01 


5.20 - 


5.70 


5.45 


0.325 


1.00 


8.01E+00 


2.07E- 


-01 


1.77E- 


-01 


5.70 - 


6.20 


5.95 


0.325 


1.02 


6.21E+00 


1.54E- 


-01 


1.53E- 


-01 


6.20 - 


6.70 


6.45 


0.325 


1.07 


5.23E+00 


1.30E- 


-01 


8.99E- 


-02 


6.70 - 


7.20 


6.95 


0.325 


1.10 


4.45E+00 


1.16E- 


-01 


8.46E- 


-02 


7.20 - 


7.70 


7.44 


0.325 


1.15 


3.57E+00 


9.93E- 


-02 


7.96E- 


-02 



Table 5.2: (continued) 



Mass [GeV] < M > < x F > < p T > M 3 d 2 a/dMdx F Stat. Error Syst. Error 



7.70 - 


8.20 


7.93 


0.325 


1.22 


2.98E+00 


9.51E- 


-02 


5.55E- 


-02 


8.20 - 


8.70 


8.44 


0.325 


1.19 


2.48E+00 


8.81E- 


-02 


4.77E- 


-02 


10.85 - 


11.85 


11.29 


0.326 


1.11 


5.01E-01 


4.47E- 


-02 


2.97E- 


-02 


11.85 - 


12.85 


12.17 


0.326 


1.18 


5.05E-01 


1.10E- 


-01 


3.74E- 


-02 


12.85 - 


14.85 


13.41 


0.328 


1.16 


3.18E-01 


1.06E- 


-01 


1.92E- 


-02 



0.35 <x F < 0.40 



4.20 - 


4.70 


4.45 


0.375 


1 


00 


7.97E+00 


2.55E- 


-01 


1 


80E- 


-01 


4.70 - 


5.20 


4.96 


0.374 


1 


01 


6.47E+00 


1.93E- 


-01 


1 


81E- 


-01 


5.20 - 


5.70 


5.46 


0.375 





99 


5.90E+00 


1.58E- 


-01 


9 


73E- 


-02 


5.70 - 


6.20 


5.95 


0.374 


1 


02 


4.91E+00 


1.20E- 


-01 


8 


34E- 


-02 


6.20 - 


6.70 


6.44 


0.374 


1 


06 


4.08E+00 


1.03E- 


-01 


6 


99E- 


-02 


6.70 - 


7.20 


6.94 


0.374 


1 


10 


3.50E+00 


9.44E- 


-02 


9 


94E- 


-02 


7.20 - 


7.70 


7.44 


0.374 


1 


16 


2.92E+00 


8.24E- 


-02 


1 


01E- 


-01 


7.70 - 


8.20 


7.94 


0.375 


1 


15 


2.51E+00 


8.00E- 


-02 


1 


10E- 


-01 


8.20 - 


8.70 


8.44 


0.374 


1 


18 


2.08E+00 


7.96E- 


-02 


8 


28E- 


-02 


10.85 - 


11.85 


11.27 


0.375 


1 


16 


5.13E-01 


4.04E- 


-02 


4 


33E- 


-03 


11.85 - 


12.85 


12.29 


0.376 


1 


17 


2.05E-01 


3.65E- 


-02 


1 


12E- 


-02 


12.85 - 


14.85 


13.31 


0.378 


1 


08 


1.94E-01 


4.76E- 


-02 


1 


20E- 


-02 



Table 5.2: (continued) 



Mass [GeV] < M > < x F > < p T > M 3 d 2 a/dMdx F Stat. Error Syst. Error 



14.85 - 16.85 15.48 0.390 1.64 2.40E-01 1.77E-01 1.76E-01 













0.40 < x F 


< 0.45 












4 


20 - 


4.70 


4.45 


0.423 


1.03 


6.09E+00 


2.22E- 


-01 


1 


68E- 


-01 


4 


70 - 


5.20 


4.96 


0.424 


0.99 


5.09E+00 


1.58E- 


-01 


1 


59E- 


-01 


5 


20 - 


5.70 


5.46 


0.424 


0.99 


4.28E+00 


1.22E- 


-01 


1 


29E- 


-01 


5 


70 - 


6.20 


5.95 


0.424 


1.00 


3.78E+00 


9.71E- 


-02 


1 


10E- 


-01 


6 


20 - 


6.70 


6.45 


0.424 


1.02 


3.12E+00 


8.00E- 


-02 


7 


37E- 


-02 


6 


70 - 


7.20 


6.93 


0.424 


1.08 


2.60E+00 


7.43E- 


-02 


6 


73E- 


-02 


7 


20 - 


7.70 


7.44 


0.424 


1.12 


2.10E+00 


6.57E- 


-02 


4 


68E- 


-02 


7 


70 - 


8.20 


7.94 


0.425 


1.13 


1.86E+00 


6.46E- 


-02 


6 


32E- 


-02 


8 


20 - 


8.70 


8.44 


0.424 


1.18 


1.63E+00 


6.21E- 


-02 


4 


68E- 


-02 


10 


85 - 


11.85 


11.26 


0.425 


1.17 


4.33E-01 


3.40E- 


-02 


1 


10E- 


-02 


11 


85 - 


12.85 


12.23 


0.427 


1.27 


2.70E-01 


3.32E- 


-02 


1 


10E- 


-02 


12 


85 - 


14.85 


13.56 


0.426 


1.15 


1.39E-01 


2.81E- 


-02 


6 


47E- 


-03 


14 


85 - 


16.85 


15.59 


0.434 


1.29 


4.24E-01 


2.41E- 


-01 


1 


20E- 


-01 



Table 5.2: (continued) 



Mass [GeV] < M > < x F > < p T > M 3 d 2 a/dMdx F Stat. Error Syst. Error 



0.45 <x F < 0.50 



4 


20 - 


4.70 


4.46 


0.473 


1 


00 


4.49E+00 


2 


03E- 


-01 


3 


45E- 


-01 


4 


70 - 


5.20 


4.96 


0.474 





97 


3.68E+00 


1 


31E- 


-01 


1 


16E- 


-01 


5 


20 - 


5.70 


5.46 


0.474 


1 


01 


3.07E+00 


9 


86E- 


-02 


6 


13E- 


-02 


5 


70 - 


6.20 


5.95 


0.474 





99 


2.75E+00 


7 


76E- 


-02 


5 


63E- 


-02 


6 


20 - 


6.70 


6.44 


0.475 


1 


02 


2.29E+00 


6 


88E- 


-02 


5 


10E- 


-02 


6 


70 - 


7.20 


6.94 


0.474 


1 


07 


1.91E+00 


6 


02E- 


-02 


4 


62E- 


-02 


7 


20 - 


7.70 


7.44 


0.473 


1 


16 


1.57E+00 


5 


62E- 


-02 


3 


89E- 


-02 


7 


70 - 


8.20 


7.93 


0.474 


1 


13 


1.33E+00 


5 


24E- 


-02 


2 


79E- 


-02 


8 


20 - 


8.70 


8.44 


0.474 


1 


08 


1.21E+00 


7 


74E- 


-02 


1 


68E- 


-02 


10 


85 - 


11.85 


11.31 


0.476 


1 


19 


3.62E-01 


3 


11E- 


-02 


1 


06E- 


-02 


11 


85 - 


12.85 


12.26 


0.475 


1 


09 


2.01E-01 


2 


96E- 


-02 


7 


79E- 


-03 


12 


85 - 


14.85 


13.43 


0.471 


1 


08 


1.01E-01 


2 


71E- 


-02 


4 


90E- 


-03 


14 


85 - 


16.85 


16.62 


0.466 





96 


7.46E-02 


7 


46E- 


-02 


2 


04E- 


-02 



0.50 <x F < 0.55 

4.20 - 4.70 4.46 0.523 0.97 3.08E+00 1.60E-01 1.55E-01 



Table 5.2: (continued) 



Mass [GeV] < M > < x F > < p T > M 3 d 2 a/dMdx F Stat. Error Syst. Error 



4 


70 - 


5.20 


4.97 





525 





99 


2 


57E+00 


1.11E- 


-01 


1.79E- 


-01 


5 


20 - 


5.70 


5.46 





524 





95 


2 


08E+00 


7.41E- 


-02 


7.49E- 


-02 


5 


70 - 


6.20 


5.95 





524 





97 


1 


95E+00 


6.28E- 


-02 


5.96E- 


-02 


6 


20 - 


6.70 


6.45 





524 





97 


1 


66E+00 


5.34E- 


-02 


3.90E- 


-02 


6 


70 - 


7.20 


6.95 





524 


1 


07 


1 


29E+00 


4.83E- 


-02 


3.00E- 


-02 


7 


20 - 


7.70 


7.44 





523 


1 


05 


1 


16E+00 


4.65E- 


-02 


2.67E- 


-02 


7 


70 - 


8.20 


7.93 





523 


1 


11 


9 


72E-01 


4.40E- 


-02 


3.70E- 


-02 


8 


20 - 


8.70 


8.44 





524 


1 


16 


8 


48E-01 


4.03E- 


-02 


2.75E- 


-02 


10 


85 - 


11.85 


11.30 





522 


1 


14 


2 


00E-01 


2.14E- 


-02 


8.15E- 


-03 


11 


85 - 


12.85 


12.27 





527 


1 


03 


1 


23E-01 


2.46E- 


-02 


8.66E- 


-03 


12 


85 - 


14.85 


13.55 





524 


1 


19 


8 


12E-02 


1.94E- 


-02 


3.91E- 


-03 



0.55 <x F < 0.60 



4.20 - 


4.70 


4.47 


0.572 


1.04 


2.17E+00 


1.37E- 


-01 


9.50E- 


-02 


4.70 - 


5.20 


4.97 


0.574 


0.96 


1.75E+00 


8.45E- 


-02 


6.55E- 


-02 


5.20 - 


5.70 


5.45 


0.574 


1.00 


1.51E+00 


6.37E- 


-02 


4.53E- 


-02 


5.70 - 


6.20 


5.95 


0.574 


0.96 


1.29E+00 


4.89E- 


-02 


2.80E- 


-02 


6.20 - 


6.70 


6.44 


0.573 


1.02 


1.07E+00 


4.29E- 


-02 


2.27E- 


-02 


6.70 - 


7.20 


6.94 


0.572 


1.02 


9.56E-01 


3.99E- 


-02 


2.22E- 


-02 



Table 5.2: (continued) 



Mass [GeV] < M > < x F > < p T > M 3 d 2 a/dMdx F Stat. Error Syst. Error 



7.20 - 


7.70 


7.44 


0.573 


1.08 


7.75E- 


-01 


3.52E- 


-02 


1.78E- 


-02 


7.70 - 


8.20 


7.93 


0.573 


1.04 


6.78E- 


-01 


3.54E- 


-02 


2.12E- 


-02 


8.20 - 


8.70 


8.44 


0.573 


1.10 


6.17E- 


-01 


3.44E- 


-02 


2.01E- 


-02 


10.85 - 


11.85 


11.26 


0.574 


1.12 


1.48E- 


-01 


1.71E- 


-02 


2.14E- 


-03 


11.85 - 


12.85 


12.35 


0.574 


1.23 


1.11E- 


-01 


1.96E- 


-02 


5.64E- 


-03 


12.85 - 


14.85 


13.51 


0.572 


1.11 


6.24E- 


-02 


1.77E- 


-02 


8.75E- 


-03 


14.85 - 


16.85 


15.19 


0.597 


0.56 


6.89E- 


-03 


6.89E- 


-03 


2.44E- 


-03 



0.60 < x F < 0.65 



4.20 - 


4.70 


4.48 


0.624 


1.01 


1.13E-f 


-00 


9.93E- 


-02 


5.28E- 


-02 


4.70 - 


5.20 


4.96 


0.623 


0.92 


1.07E+00 


6.80E- 


-02 


4.01E- 


-02 


5.20 - 


5.70 


5.47 


0.624 


0.90 


8.64E- 


-01 


4.54E- 


-02 


3.58E- 


-02 


5.70 - 


6.20 


5.96 


0.622 


0.94 


7.77E- 


-01 


4.28E- 


-02 


2.64E- 


-02 


6.20 - 


6.70 


6.45 


0.623 


0.92 


6.47E- 


-01 


3.39E- 


-02 


2.70E- 


-02 


6.70 - 


7.20 


6.94 


0.622 


0.99 


5.81E- 


-01 


3.25E- 


-02 


1.73E- 


-02 


7.20 - 


7.70 


7.44 


0.624 


1.08 


4.47E- 


-01 


3.15E- 


-02 


1.69E- 


-02 


7.70 - 


8.20 


7.96 


0.622 


1.03 


4.30E- 


-01 


3.25E- 


-02 


1.30E- 


-02 


8.20 - 


8.70 


8.45 


0.625 


1.09 


3.49E- 


-01 


2.76E- 


-02 


9.92E- 


-03 


10.85 - 


11.85 


11.28 


0.621 


1.19 


9.51E- 


-02 


1.42E- 


-02 


3.28E- 


-03 



Table 5.2: (continued) 



Mass [GeV] < M > < x F > < p T > M 3 d 2 a/dMdx F Stat. Error Syst. Error 



11.85 - 12.85 12.33 0.621 1.21 8.06E-02 1.61E-02 6.38E-03 

12.85 - 14.85 13.38 0.622 1.25 2.90E-02 7.26E-03 1.73E-03 

14.85 - 16.85 15.89 0.623 3.48 4.67E-02 4.67E-02 1.78E-01 



0.65 <x F < 0.70 



4.20 - 


4.70 


4.48 


0.669 


0.95 


5.85E- 


-01 


7.05E- 


-02 


5.31E- 


-02 


4.70 - 


5.20 


4.96 


0.670 


0.93 


6.19E- 


-01 


4.74E- 


-02 


3.13E- 


-02 


5.20 - 


5.70 


5.47 


0.673 


0.90 


5.69E- 


-01 


3.99E- 


-02 


2.79E- 


-02 


5.70 - 


6.20 


5.96 


0.673 


0.88 


4.49E- 


-01 


2.94E- 


-02 


1.59E- 


-02 


6.20 - 


6.70 


6.45 


0.673 


1.02 


4.01E- 


-01 


2.91E- 


-02 


1.33E- 


-02 


6.70 - 


7.20 


6.95 


0.673 


0.99 


3.65E- 


-01 


2.89E- 


-02 


1.63E- 


-02 


7.20 - 


7.70 


7.44 


0.671 


1.11 


2.81E- 


-01 


2.54E- 


-02 


1.25E- 


-02 


7.70 - 


8.20 


7.93 


0.671 


1.10 


2.44E- 


-01 


3.18E- 


-02 


1.53E- 


-02 


8.20 - 


8.70 


8.43 


0.670 


1.06 


1.91E- 


-01 


2.09E- 


-02 


9.33E- 


-03 


10.85 - 


11.85 


11.33 


0.674 


0.97 


5.01E- 


-02 


1.46E- 


-02 


3.46E- 


-03 


11.85 - 


12.85 


12.26 


0.675 


0.95 


4.39E- 


-02 


1.34E- 


-02 


4.85E- 


-03 


12.85 - 


14.85 


13.52 


0.679 


0.95 


2.01E- 


-02 


8.56E- 


-03 


5.89E- 


-03 



Table 5.2: (continued) 



Mass [GeV] < M > < x F > < p T > M 3 d 2 a/dMdx F Stat. Error Syst. Error 











0.70 < x F 


< 0.75 












4.20 - 


4.70 


4.50 


0.721 


1.03 


3.58E- 


-01 


6.46E- 


-02 


3.68E- 


-02 


4.70 - 


5.20 


4.96 


0.722 


0.91 


3.17E- 


-01 


3.93E- 


-02 


2.03E- 


-02 


5.20 - 


5.70 


5.47 


0.722 


0.87 


2.08E- 


-01 


2.19E- 


-02 


1.06E- 


-02 


5.70 - 


6.20 


5.94 


0.723 


0.82 


2.31E- 


-01 


3.15E- 


-02 


1.44E- 


-02 


6.20 - 


6.70 


6.43 


0.723 


0.90 


1.88E- 


-01 


2.83E- 


-02 


1.72E- 


-02 


6.70 - 


7.20 


6.95 


0.722 


0.98 


1.69E- 


-01 


1.82E- 


-02 


8.32E- 


-03 


7.20 - 


7.70 


7.45 


0.722 


1.00 


1.35E- 


-01 


1.79E- 


-02 


7.68E- 


-03 


7.70 - 


8.20 


7.93 


0.721 


0.98 


1.00E- 


-01 


1.42E- 


-02 


5.09E- 


-03 


8.20 - 


8.70 


8.42 


0.722 


0.98 


8.28E- 


-02 


1.37E- 


-02 


4.85E- 


-03 


10.85 - 


11.85 


11.28 


0.719 


1.16 


3.11E- 


-02 


7.74E- 


-03 


3.80E- 


-03 


11.85 - 


12.85 


12.25 


0.733 


0.74 


3.29E- 


-02 


1.42E- 


-02 


3.08E- 


-03 


12.85 - 


14.85 


13.68 


0.717 


1.73 


4.94E- 


-03 


4.94E- 


-03 


1.78E- 


-03 



0.75 <x F < 0.80 

4.20 - 4.70 4.47 0.777 1.20 1.89E-01 6.62E-02 2.48E-02 

4.70 - 5.20 4.95 0.773 0.90 1.52E-01 2.82E-02 1.37E-02 



Table 5.2: (continued) 



Mass [GeV] < M > < x F > < p T > M 3 d 2 a/dMdx F Stat. Error Syst. Error 



5.20 - 


5.70 


5.44 


0.769 


0.75 


1.02E- 


-01 


1.62E- 


-02 


6.22E- 


-03 


5.70 - 


6.20 


5.98 


0.770 


0.83 


1.19E- 


-01 


2.45E- 


-02 


8.09E- 


-03 


6.20 - 


6.70 


6.45 


0.772 


0.90 


1.27E- 


-01 


5.28E- 


-02 


8.69E- 


-03 


6.70 - 


7.20 


6.92 


0.769 


0.93 


8.04E- 


-02 


1.20E- 


-02 


6.24E- 


-03 


7.20 - 


7.70 


7.47 


0.769 


0.83 


4.95E- 


-02 


1.09E- 


-02 


5.26E- 


-03 


7.70 - 


8.20 


7.88 


0.769 


1.06 


6.40E- 


-02 


1.42E- 


-02 


6.45E- 


-03 


8.20 - 


8.70 


8.44 


0.775 


0.93 


3.33E- 


-02 


1.06E- 


-02 


1.48E- 


-03 


10.85 - 


11.85 


11.30 


0.770 


1.01 


1.36E- 


-02 


6.22E- 


-03 


2.71E- 


-03 


11.85 - 


12.85 


12.30 


0.765 


1.35 


2.08E- 


-02 


1.20E- 


-02 


3.21E- 


-03 


12.85 - 


14.85 


13.59 


0.773 


0.61 


1.09E- 


-02 


1.09E- 


-02 


4.06E- 


-03 
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Figure 5.1: Scaling form M 3 d 2 cr / 'dMdxp of the hydrogen (left) and deuterium (right) cross sections 
for —0.05 < x F < 0.05. The data are compared to NLO calculations based on several different sets of 
parton distributions. The inset shows the ratio of the experimentally measured cross sections to NLO 
calculations based on the CTEQ 5 partons. 
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Figure 5.2: Scaling form M 3 d 2 'a j 'dM dx p of the hydrogen (left) and deuterium (right) cross sections for 
0.05 < x F < 0.1. The data are compared to NLO calculations based on several different sets of parton 
distributions. The inset shows the ratio of the experimentally measured cross sections to NLO calculations 
based on the CTEQ 5 partons. 



10 



a 

J3 




E866 
CTEQ5 
CTEQ6 
MRST 98 
MRST 01 
GRV 



10 



a 

J3 



x = 0.125 

F 




E866 
CTEQ5 
CTEQ6 
MRST 98 
MRST 01 
GRV 



12 



16 




12 



16 



M 



M 



Figure 5.3: Scaling form M 3 d 2 'a j 'dM dx p of the hydrogen (left) and deuterium (right) cross sections for 
0.1 < x F < 0.15. The data are compared to NLO calculations based on several different sets of parton 
distributions. The inset shows the ratio of the experimentally measured cross sections to NLO calculations 
based on the CTEQ 5 partons. 
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Figure 5.4: Scaling form M 3 d 2 a/dMdxp of the hydrogen (left) and deuterium (right) cross sections for 
0.15 < x F < 0.2. The data are compared to NLO calculations based on several different sets of parton 
distributions. The inset shows the ratio of the experimentally measured cross sections to NLO calculations 
based on the CTEQ 5 partons. 
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Figure 5.5: Scaling form M 3 d 2 'a / 'dM dx p of the hydrogen (left) and deuterium (right) cross sections for 
0.2 < x F < 0.25. The data are compared to NLO calculations based on several different sets of parton 
distributions. The inset shows the ratio of the experimentally measured cross sections to NLO calculations 
based on the CTEQ 5 partons. 
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Figure 5.6: Scaling form M 3 d 2 a/dMdxp of the hydrogen (left) and deuterium (right) cross sections for 
0.25 < x F < 0.3. The data are compared to NLO calculations based on several different sets of parton 
distributions. The inset shows the ratio of the experimentally measured cross sections to NLO calculations 
based on the CTEQ 5 partons. 
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Figure 5.7: Scaling form M 3 d 2 'a / ' dM dx p of the hydrogen (left) and deuterium (right) cross sections for 
0.3 < x F < 0.35. The data are compared to NLO calculations based on several different sets of parton 
distributions. The inset shows the ratio of the experimentally measured cross sections to NLO calculations 
based on the CTEQ 5 partons. 




Figure 5.8: Scaling form M 3 d 2 a/dMdxp of the hydrogen (left) and deuterium (right) cross sections for 
0.35 < x F < 0.4. The data are compared to NLO calculations based on several different sets of parton 
distributions. The inset shows the ratio of the experimentally measured cross sections to NLO calculations 
based on the CTEQ 5 partons. 
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Figure 5.9: Scaling form M 3 d 2 a/dMdxF of the hydrogen (left) and deuterium (right) cross sections for 
0.4 < xf < 0.45. The data are compared to NLO calculations based on several different sets of parton 
distributions. The inset shows the ratio of the experimentally measured cross sections to NLO calculations 
based on the CTEQ 5 partons. 
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Figure 5.10: Scaling form M 3 d 2 a / 'dMdxp of the hydrogen (left) and deuterium (right) cross sections 
for 0.45 < x F < 0.5. The data are compared to NLO calculations based on several different sets of 
parton distributions. The inset shows the ratio of the experimentally measured cross sections to NLO 
calculations based on the CTEQ 5 partons. 




Figure 5.11: Scaling form M 3 d 2 a / 'dMdxp of the hydrogen (left) and deuterium (right) cross sections 
for 0.5 < x F < 0.55. The data are compared to NLO calculations based on several different sets of 
parton distributions. The inset shows the ratio of the experimentally measured cross sections to NLO 
calculations based on the CTEQ 5 partons. 
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Figure 5.12: Scaling form M 3 d 2 a/dMdxp of the hydrogen (left) and deuterium (right) cross sections 
for 0.55 < x F < 0.6. The data are compared to NLO calculations based on several different sets of 
parton distributions. The inset shows the ratio of the experimentally measured cross sections to NLO 
calculations based on the CTEQ 5 partons. 
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Figure 5.13: Scaling form M 3 d 2 a/dMdxp of the hydrogen (left) and deuterium (right) cross sections 
for 0.6 < x F < 0.65. The data are compared to NLO calculations based on several different sets of 
parton distributions. The inset shows the ratio of the experimentally measured cross sections to NLO 
calculations based on the CTEQ 5 partons. 
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Figure 5.14: Scaling form M 3 d 2 a/dMdxp of the hydrogen (left) and deuterium (right) cross sections 
for 0.65 < x F < 0.7. The data are compared to NLO calculations based on several different sets of 
parton distributions. The inset shows the ratio of the experimentally measured cross sections to NLO 
calculations based on the CTEQ 5 partons. 
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Figure 5.15: Scaling form M 3 d 2 a / 'dMdxp of the hydrogen (left) and deuterium (right) cross sections 
for 0.7 < x F < 0.75. The data are compared to NLO calculations based on several different sets of 
parton distributions. The inset shows the ratio of the experimentally measured cross sections to NLO 
calculations based on the CTEQ 5 partons. 




Figure 5.16: Scaling form M 3 d 2 a/dMdxp of the hydrogen (left) and deuterium (right) cross sections 
for 0.75 < x F < 0.8. The data are compared to NLO calculations based on several different sets of 
parton distributions. The inset shows the ratio of the experimentally measured cross sections to NLO 
calculations based on the CTEQ 5 partons. 



5.2 FNAL E866/NuSea Results - Ed 3 a /dp 3 

The pr-distributions are presented in terms of the invariant cross section 

E fa = 2E d 2 a 
dp 3 TCy/sdxFdp^ 

where the dependence has been averaged over. The cross sections are tabu- 
lated in tables 5.3 through 5.10, in xp bins, integrated over the specified range in 
invariant mass. The mean mass, xf, Pt and E in each bin are also tabulated. 
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Table 5.3: Invariant form Ed 3 a/dp 3 (in units of pb / GeV 2 / nucleon) for the hydrogen 
cross section over the range -0.05 < xf < 0.15. Statistical and point-to-point systematic 
uncertainties are shown separately. An additional ±6.5% uncertainty in the normalization 
is common to all data points. 



p T [GeV] < p T > < M > < x F > < E > Ed 3 a /dp 3 Stat. Error. Syst. Error 



4.20 < M„+„- < 5.20 






00- 





25 





16 


4 


64 





11 


6.277 


1.33E+01 


1.56E4 


-00 


8.80E- 


-01 





25- 





50 





39 


4 


66 





12 


6.544 


1.21E+01 


1.20E4 


-00 


3.62E- 


-01 





50- 





75 





63 


4 


65 





12 


6.533 


9.72E+00 


8.45E- 


-01 


2.93E- 


-01 





75- 


1 


00 





87 


4 


64 





12 


6.641 


8.85E+00 


6.40E- 


-01 


2.80E- 


-01 


1 


00- 


1 


25 


1 


12 


4 


61 





12 


6.727 


5.70E+00 


5.60E- 


-01 


1.84E- 


-01 


1 


25- 


1 


50 


1 


36 


4 


64 





12 


6.741 


4.86E+00 


4.72E- 


-01 


1.57E- 


-01 


1 


50- 


1 


75 


1 


64 


4 


61 





13 


6.923 


2.01E+00 


3.70E- 


-01 


8.46E- 


-02 


1 


75- 


2 


00 


1 


85 


4 


58 





12 


6.782 


1.67E+00 


3.43E- 


-01 


7.15E- 


-02 


2 


00- 


2 


50 


2 


19 


4 


54 





12 


6.866 


7.72E-01 


1.45E- 


-01 


3.40E- 


-02 


2 


50- 


3 


00 


2 


69 


4 


77 





13 


7.434 


3.06E-01 


1.00E- 


-01 


1.91E- 


-02 


3 


00- 


■3 


50 


3 


13 


4 


75 





13 


7.657 


1.22E-01 


7.95E- 


-02 


1.34E- 


-02 



Table 5.3: (continued) 



p T [GeV] < pr> < M > < x F > < E > Ed 3 a /dp 3 Stat. Error. Syst. Error 









5.20 


< M„+„- 


< 6.20 










0.00-0.25 


0.19 


5.61 


0.13 


7.449 


1.97E+00 


8.11E- 


-01 


2.00E- 


-01 


0.25-0.50 


0.39 


5.67 


0.12 


7.305 


4.75E+00 


6.04E- 


-01 


1.80E- 


-01 


0.50-0.75 


0.62 


5.56 


0.12 


7.221 


4.72E+00 


4.68E- 


-01 


1.74E- 


-01 


0.75-1.00 


0.88 


5.65 


0.12 


7.475 


2.88E+00 


3.85E- 


-01 


1.11E- 


-01 


1.00-1.25 


1.12 


5.66 


0.12 


7.496 


2.84E+00 


3.24E- 


-01 


1.10E- 


-01 


1.25-1.50 


1.35 


5.61 


0.13 


7.566 


1.01E+00 


2.40E- 


-01 


7.08E- 


-02 


1.50-1.75 


1.60 


5.70 


0.12 


7.498 


1.05E+00 


2.01E- 


-01 


8.97E- 


-02 


1.75-2.00 


1.86 


5.66 


0.12 


7.577 


8.38E-01 


1.92E- 


-01 


7.77E- 


-02 


2.00-2.50 


2.16 


5.67 


0.11 


7.518 


1.73E-01 


6.68E- 


-02 


2.25E- 


-02 


2.50-3.00 


2.76 


5.68 


0.11 


7.656 


2.45E-01 


8.65E- 


-02 


2.72E- 


-02 


3.00-3.50 


3.28 


5.79 


0.12 


8.119 


5.91E-02 


4.21E- 


-02 


1.32E- 


-02 



6.20 < M M+ ^- < 7.20 

0.00-0.25 0.16 6.66 0.10 7.877 1.58E+00 3.93E-01 1.59E-01 

0.25-0.50 0.38 6.67 0.09 7.804 1.52E+00 3.17E-01 8.64E-02 

0.50-0.75 0.62 6.69 0.10 7.865 1.42E+00 2.43E-01 9.54E-02 



Table 5.3: (continued) 



p T [GeV] < pr> < M > < x F > < E > Ed 3 a /dp 3 Stat. Error. Syst. Error 



0.75-1.00 


0.88 


6.70 


0.10 


7.982 


1.11E+00 


1.60E- 


-01 


7.37E- 


-02 


1.00-1.25 


1.10 


6.75 


0.11 


8.178 


8.44E-01 


1.63E- 


-01 


9.10E- 


-02 


1.25-1.50 


1.36 


6.66 


0.12 


8.212 


3.80E-01 


1.38E- 


-01 


6.89E- 


-02 


1.50-1.75 


1.62 


6.75 


0.11 


8.266 


2.92E-01 


1.16E- 


-01 


5.42E- 


-02 


1.75-2.00 


1.84 


6.80 


0.10 


8.086 


1.32E-01 


5.25E- 


-02 


1.00E- 


-02 


2.00-2.50 


2.16 


6.56 


0.12 


8.390 


1.22E-01 


6.79E- 


-02 


2.00E- 


-02 


2.50-3.00 


2.56 


6.47 


0.12 


8.333 


5.16E-02 


3.09E- 


-02 


1.01E- 


-02 



7.20 < M„+„- < 8.70 



0.00- 


■0.25 


0.16 


8.14 


0.09 


8.988 


8.26E- 


-01 


1.35E- 


-01 


5.02E- 


-02 


0.25- 


-0.50 


0.39 


8.09 


0.09 


9.029 


1.08E+00 


9.70E- 


-02 


3.41E- 


-02 


0.50- 


-0.75 


0.63 


8.05 


0.09 


8.982 


7.51E- 


-01 


6.72E- 


-02 


3.68E- 


-02 


0.75- 


-1.00 


0.87 


8.06 


0.10 


9.072 


7.33E- 


-01 


5.25E- 


-02 


2.39E- 


-02 


1.00- 


■1.25 


1.12 


8.02 


0.10 


9.102 


5.26E- 


-01 


4.69E- 


-02 


2.08E- 


-02 


1.25- 


■1.50 


1.36 


8.05 


0.10 


9.179 


3.09E- 


-01 


3.35E- 


-02 


1.27E- 


-02 


1.50- 


■1.75 


1.61 


8.10 


0.10 


9.265 


2.39E- 


-01 


2.96E- 


-02 


1.03E- 


-02 


1.75- 


■2.00 


1.86 


8.05 


0.11 


9.386 


1.50E- 


-01 


2.55E- 


-02 


7.80E- 


-03 


2.00- 


■2.50 


2.17 


8.05 


0.11 


9.397 


5.58E- 


-02 


1.15E- 


-02 


2.08E- 


-03 


2.50- 


■3.00 


2.66 


8.10 


0.11 


9.569 


3.27E- 


-02 


8.53E- 


-03 


2.38E- 


-03 



Table 5.3: (continued) 



Pt [GeV] 


< Pt > 


< M > 


< Xf > 


< E > 


Ed 3 a /dp 3 


Stat. Error. 


Syst. Error 


3.00-3.50 
3.50-4.00 


3.26 
3.78 


7.98 
7.89 


0.10 
0.11 


9.524 
9.762 


1.22E-02 
1.88E-02 


4.51E-03 
9.54E-03 


1.88E-03 
3.18E-03 





10.85 < M„+„- < 12.85 



0.25-0.50 


0.40 


11.25 


0.08 


11.627 


6.61E- 


-02 


6.61E- 


-02 


1.55E- 


-02 


0.50-0.75 


0.64 


11.35 


0.13 


12.525 


7.12E- 


-02 


3.07E- 


-02 


2.13E- 


-02 


0.75-1.00 


0.89 


11.28 


0.10 


11.989 


6.28E- 


-02 


3.63E- 


-02 


1.21E- 


-02 


1.00-1.25 


1.16 


11.04 


0.11 


11.870 


4.23E- 


-02 


1.97E- 


-02 


4.45E- 


-03 


1.25-1.50 


1.37 


10.98 


0.11 


11.866 


2.99E- 


-02 


1.50E- 


-02 


4.27E- 


-03 


1.50-1.75 


1.59 


11.23 


0.13 


12.355 


3.57E- 


-02 


2.52E- 


-02 


6.28E- 


-03 


1.75-2.00 


1.89 


11.25 


0.10 


12.048 


9.40E- 


-02 


4.70E- 


-02 


1.96E- 


-02 



Table 5.4: Invariant form Ed 3 a/dp 3 (in units of pb / GeV 2 / nucleon) for the hydrogen 
cross section over the range 0.15 < xp < 0.35. Statistical and point-to-point systematic 
uncertainties are shown. An additional ±6.5% uncertainty in the normalization is common 
to all data points. 



p T [GeV] < p T > < M > < x F > < E > Ed 3 a /dp 3 Stat. Error. Syst. Error 



4.20 < M„+„- < 5.20 






00- 





25 





17 


4 


67 





25 


10.767 


1 


15E+01 


7 


46E- 


-01 


4 


71E- 


-01 





25- 





50 





39 


4 


66 





25 


10.840 


1 


31E+01 


5 


19E- 


-01 


3 


94E- 


-01 





50- 





75 





63 


4 


66 





25 


10.848 


1 


16E+01 


3 


96E- 


-01 


3 


20E- 


-01 





75- 


1 


00 





87 


4 


66 





25 


10.880 


8 


05E+00 


2 


91E- 


-01 


2 


22E- 


-01 


1 


00- 


1 


25 


1 


12 


4 


66 





25 


10.807 


5 


99E+00 


2 


21E- 


-01 


1 


70E- 


-01 


1 


25- 


1 


50 


1 


37 


4 


67 





25 


10.902 


3 


87E+00 


1 


82E- 


-01 


1 


15E- 


-01 


1 


50- 


1 


75 


1 


61 


4 


65 





25 


10.877 


2 


54E+00 


1 


39E- 


-01 


6 


89E- 


-02 


1 


75- 


2 


00 


1 


87 


4 


68 





25 


11.050 


1 


49E+00 


1 


06E- 


-01 


5 


25E- 


-02 


2 


00- 


2 


50 


2 


21 


4 


63 





26 


11.252 


6 


55E-01 


5 


78E- 


-02 


2 


25E- 


-02 


2 


50- 


3 


00 


2 


72 


4 


62 





26 


11.532 


3 


01E-01 


3 


44E- 


-02 


2 


07E- 


-02 


3 


00- 


3 


50 


3 


21 


4 


69 





26 


11.599 


1 


07E-01 


3 


09E- 


-02 


1 


33E- 


-02 


3 


50- 


4 


00 


3 


72 


4 


65 





27 


12.160 


5 


40E-02 


2 


34E- 


-02 


1 


28E- 


-02 


4 


00- 


5 


00 


4 


32 


4 


71 





28 


12.512 


3 


92E-02 


1 


39E- 


-02 


5 


32E- 


-03 



Table 5.4: (continued) 



p T [GeV] < pr> < M > < x F > < E > Ed 3 a /dp 3 Stat. Error. Syst. Error 



5.20 < M.+U- < 6.20 






00- 





25 





16 


5 


69 





27 


12.024 


4 


26E+00 


3 


56E- 


-01 


1 


85E- 


-01 





25- 





50 





39 


5 


68 





27 


12.053 


5 


55E+00 


2 


39E- 


-01 


1 


86E- 


-01 





50- 





75 





63 


5 


69 





27 


12.023 


4 


32E+00 


1 


68E- 


-01 


1 


37E- 


-01 





75- 


1 


00 





87 


5 


69 





27 


11.996 


3 


25E+00 


1 


32E- 


-01 


9 


49E- 


-02 


1 


00- 


1 


25 


1 


12 


5 


69 





27 


12.090 


2 


34E+00 


1 


03E- 


-01 


6 


63E- 


-02 


1 


25- 


1 


50 


1 


37 


5 


70 





27 


12.065 


1 


59E+00 


8 


55E- 


-02 


4 


59E- 


-02 


1 


50- 


1 


75 


1 


62 


5 


70 





27 


12.115 


1 


11E+00 


6 


37E- 


-02 


3 


29E- 


-02 


1 


75- 


2 


00 


1 


87 


5 


71 





28 


12.359 


6 


55E-01 


5 


27E- 


-02 


2 


15E- 


-02 


2 


00- 


2 


50 


2 


21 


5 


63 





27 


12.071 


3 


15E-01 


2 


56E- 


-02 


1 


09E- 


-02 


2 


50- 


3 


00 


2 


69 


5 


70 





27 


12.127 


9 


87E-02 


1 


65E- 


-02 


4 


90E- 


-03 


3 


00- 


-3 


50 


3 


24 


5 


75 





28 


12.880 


3 


94E-02 


9 


16E- 


-03 


6 


33E- 


-03 


3 


50- 


4 


00 


3 


71 


5 


66 





28 


12.947 


2 


74E-02 


8 


02E- 


-03 


6 


96E- 


-03 


4 


00- 


5 


00 


4 


35 


5 


84 





27 


12.853 


2 


06E-02 


6 


10E- 


-03 


4 


55E- 


-03 



6.20 < M M+ ^- < 7.20 
0.00-0.25 0.16 6.69 0.27 12.511 1.81E+00 1.44E-01 1.22E-01 



Table 5.4: (continued) 



p T [GeV] < pr> < M > < x F > < E > Ed 3 a /dp 3 Stat. Error. Syst. Error 



0.25-0.50 


0.38 


6.67 


0.27 


12.631 


2.33E+00 


1.05E- 


-01 


9.07E- 


-02 


0.50-0.75 


0.63 


6.69 


0.27 


12.488 


1.82E+00 


7.45E- 


-02 


7.60E- 


-02 


0.75-1.00 


0.88 


6.68 


0.27 


12.638 


1.41E+00 


5.89E- 


-02 


6.05E- 


-02 


1.00-1.25 


1.12 


6.70 


0.27 


12.618 


1.02E+00 


4.52E- 


-02 


3.56E- 


-02 


1.25-1.50 


1.36 


6.70 


0.27 


12.667 


6.21E-01 


3.53E- 


-02 


1.81E- 


-02 


1.50-1.75 


1.61 


6.73 


0.27 


12.751 


4.55E-01 


2.84E- 


-02 


1.88E- 


-02 


1.75-2.00 


1.86 


6.71 


0.27 


12.678 


2.56E-01 


2.12E- 


-02 


1.34E- 


-02 


2.00-2.50 


2.21 


6.72 


0.28 


12.922 


1.26E-01 


1.02E- 


-02 


4.00E- 


-03 


2.50-3.00 


2.69 


6.78 


0.27 


12.943 


4.26E-02 


5.60E- 


-03 


2.05E- 


-03 


3.00-3.50 


3.25 


6.71 


0.29 


13.499 


2.09E-02 


3.86E- 


-03 


1.38E- 


-03 


3.50-4.00 


3.71 


6.74 


0.25 


12.325 


9.14E-03 


2.80E- 


-03 


8.29E- 


-04 


4.00-5.00 


4.40 


6.91 


0.26 


13.002 


5.71E-03 


2.03E- 


-03 


1.02E- 


-03 


5.00-6.00 


5.74 


6.82 


0.32 


15.200 


2.79E-03 


2.79E- 


-03 


1.59E- 


-03 



7.20 < M„+„- < 8.70 



0.00-0.25 0.17 7.87 

0.25-0.50 0.39 7.87 

0.50-0.75 0.63 7.86 

0.75-1.00 0.87 7.87 



0.26 12.844 8.59E-01 

0.26 12.854 1.22E+00 

0.26 12.931 9.48E-01 

0.26 12.878 7.37E-01 



7.78E-02 3.89E-02 

5.56E-02 3.62E-02 

4.04E-02 2.73E-02 

3.02E-02 3.04E-02 



Table 5.4: (continued) 



p T [GeV] < pr> < M > < x F > < E > Ed 3 a /dp 3 Stat. Error. Syst. Error 



1.00-1.25 


1.12 


7.88 


0.26 


12.846 


5.22E- 


-01 


2.29E- 


-02 


1.94E- 


-02 


1.25-1.50 


1.37 


7.88 


0.26 


12.835 


3.86E- 


-01 


1.79E- 


-02 


1.07E- 


-02 


1.50-1.75 


1.62 


7.91 


0.26 


12.933 


2.35E- 


-01 


1.34E- 


-02 


7.07E- 


-03 


1.75-2.00 


1.86 


7.92 


0.26 


13.026 


1.61E- 


-01 


1.05E- 


-02 


4.82E- 


-03 


2.00-2.50 


2.21 


7.93 


0.26 


13.090 


8.23E- 


-02 


5.21E- 


-03 


3.49E- 


-03 


2.50-3.00 


2.71 


7.95 


0.26 


13.225 


3.39E- 


-02 


3.25E- 


-03 


2.56E- 


-03 


3.00-3.50 


3.23 


8.04 


0.27 


13.744 


1.33E- 


-02 


1.98E- 


-03 


1.08E- 


-03 


3.50-4.00 


3.69 


7.95 


0.28 


14.005 


5.60E- 


-03 


1.57E- 


-03 


5.36E- 


-04 


4.00-5.00 


4.43 


7.79 


0.28 


14.172 


5.35E- 


-04 


3.54E- 


-04 


1.30E- 


-04 


5.00-6.00 


5.27 


8.18 


0.26 


14.262 


1.22E- 


-03 


7.15E- 


-04 


2.92E- 


-04 









10.85 


< AW- 


< 12.85 












0.00-0.25 


0.17 


11.35 


0.28 


15.856 


5.19E- 


-02 


1.79E- 


-02 


6.69E- 


-03 


0.25-0.50 


0.36 


11.55 


0.28 


15.850 


8.50E- 


-02 


1.90E- 


-02 


3.76E- 


-03 


0.50-0.75 


0.65 


11.47 


0.28 


15.745 


6.75E- 


-02 


1.37E- 


-02 


5.43E- 


-03 


0.75-1.00 


0.86 


11.52 


0.27 


15.750 


4.67E- 


-02 


9.80E- 


-03 


3.43E- 


-03 


1.00-1.25 


1.11 


11.43 


0.28 


15.749 


5.49E- 


-02 


8.02E- 


-03 


3.49E- 


-03 


1.25-1.50 


1.36 


11.53 


0.26 


15.492 


3.29E- 


-02 


5.92E- 


-03 


3.02E- 


-03 


1.50-1.75 


1.62 


11.48 


0.28 


15.898 


2.97E- 


-02 


6.12E- 


-03 


4.83E- 


-03 



Table 5.4: (continued) 



p T [GeV] < pr> < M > < x F > < E > Ed 3 a /dp 3 Stat. Error. Syst. Error 



1.75-2.00 


1.86 


11.45 


0.28 


16.073 


9.18E- 


-03 


3.09E- 


-03 


2.44E- 


-03 


2.00-2.50 


2.17 


11.39 


0.26 


15.376 


6.91E- 


-03 


1.44E- 


-03 


1.12E- 


-03 


2.50-3.00 


2.75 


11.35 


0.29 


16.260 


4.82E- 


-03 


1.21E- 


-03 


1.22E- 


-03 


3.00-3.50 


3.24 


11.60 


0.29 


16.438 


2.84E- 


-03 


1.08E- 


-03 


5.82E- 


-04 









12.85 


< M M+M - 


< 16.85 












0.00-0.25 


0.22 


13.08 


0.22 


15.641 


4.16E- 


-02 


3.71E- 


-02 


1.57E- 


-02 


0.25-0.50 


0.41 


13.47 


0.26 


16.864 


4.77E- 


-02 


3.37E- 


-02 


8.62E- 


-03 


0.50-0.75 


0.67 


13.40 


0.27 


17.111 


4.03E- 


-02 


4.03E- 


-02 


8.31E- 


-03 


0.75-1.00 


0.93 


13.09 


0.29 


17.194 


1.67E- 


-02 


1.10E- 


-02 


1.52E- 


-03 


1.00-1.25 


1.09 


13.18 


0.24 


16.172 


2.02E- 


-02 


1.23E- 


-02 


2.20E- 


-03 


1.25-1.50 


1.35 


14.11 


0.31 


18.589 


2.47E- 


-02 


1.44E- 


-02 


3.57E- 


-03 


1.50-1.75 


1.61 


13.38 


0.32 


18.430 


2.77E- 


-02 


1.96E- 


-02 


4.81E- 


-03 


2.00-2.50 


2.29 


13.53 


0.25 


16.762 


5.19E- 


-03 


3.67E- 


-03 


8.84E- 


-04 


2.50-3.00 


2.79 


13.05 


0.25 


16.712 


5.77E- 


-03 


4.09E- 


-03 


1.33E- 


-03 



Table 5.5: Invariant form Ed 3 a/dp 3 (in units of pb / GeV 2 / nucleon) for the hydrogen 
cross section over the range 0.35 < xp < 0.55. Statistical and point-to-point systematic 
uncertainties are shown. An additional ±6.5% uncertainty in the normalization is common 
to all data points. 



p T [GeV] < p T > < M > < x F > < E > Ed 3 a /dp 3 Stat. Error. Syst. Error 



4.20 < M„+„- < 5.20 






00- 





25 





17 


4 


72 





43 


17.356 


7 


12E+00 


6 


41E- 


-01 


5 


75E- 


-01 





25- 





50 





39 


4 


73 





43 


17.272 


9 


00E+00 


4 


58E- 


-01 


3 


12E- 


-01 





50- 





75 





63 


4 


71 





43 


17.439 


7 


93E+00 


3 


27E- 


-01 


2 


30E- 


-01 





75- 


1 


00 





87 


4 


74 





43 


17.357 


5 


28E+00 


2 


47E- 


-01 


1 


94E- 


-01 


1 


00- 


1 


25 


1 


12 


4 


74 





43 


17.360 


3 


76E+00 


1 


82E- 


-01 


1 


70E- 


-01 


1 


25- 


1 


50 


1 


37 


4 


74 





43 


17.470 


2 


34E+00 


1 


42E- 


-01 


1 


16E- 


-01 


1 


50- 


1 


75 


1 


61 


4 


74 





43 


17.542 


1 


25E+00 


1 


08E- 


-01 


4 


09E- 


-02 


1 


75- 


2 


00 


1 


86 


4 


74 





43 


17.555 


7 


64E-01 


7 


55E- 


-02 


2 


69E- 


-02 


2 


00- 


2 


50 


2 


20 


4 


75 





43 


17.313 


3 


03E-01 


3 


40E- 


-02 


1 


14E- 


-02 


2 


50- 


3 


00 


2 


71 


4 


70 





42 


17.237 


9 


48E-02 


2 


02E- 


-02 


4 


83E- 


-03 


3 


00- 


3 


50 


3 


23 


4 


75 





42 


17.204 


4 


71E-02 


1 


14E- 


-02 


3 


79E- 


-03 


3 


50- 


4 


00 


3 


73 


4 


73 





41 


17.056 


3 


28E-02 


1 


01E- 


-02 


7 


09E- 


-03 


4 


00- 


5 


00 


4 


32 


4 


81 





44 


18.352 


1 


12E-02 


4 


90E- 


-03 


1 


85E- 


-02 


5 


00- 


6 


00 


5 


11 


4 


74 





37 


15.808 


4 


89E-03 


4 


89E- 


-03 


9 


98E- 


-03 



Table 5.5: (continued) 



p T [GeV] < pr> < M > < x F > < E > Ed 3 a /dp 3 Stat. Error. Syst. Error 









^ 9fl 


< M _L 












n nn n 9^ 


n ifi 

U.1U 


^ 71 


fl 44 

U .tit: 


1 s rm 




9 17E- 


-01 
yj J. 


1 97F,_ 


-D1 


0.25-0.50 


0.39 


5.72 


0.44 


17.976 


3.56E+00 


1.36E- 


-01 


1.22E- 


-01 


0.50-0.75 


0.63 


5.72 


0.43 


17.825 


3.06E+00 


1.07E- 


-01 


9.12E- 


-02 


0.75-1.00 


0.87 


5.71 


0.44 


17.900 


2.02E+00 


7.66E- 


-02 


6.00E- 


-02 


1.00-1.25 


1.12 


5.71 


0.44 


18.059 


1.43E+00 


5.79E- 


-02 


3.89E- 


-02 


1.25-1.50 


1.37 


5.73 


0.44 


17.975 


9.52E-01 


4.57E- 


-02 


2.72E- 


-02 


1.50-1.75 


1.62 


5.72 


0.44 


17.940 


5.82E-01 


3.53E- 


-02 


1.74E- 


-02 


1.75-2.00 


1.87 


5.71 


0.44 


18.033 


3.54E-01 


2.64E- 


-02 


1.13E- 


-02 


2.00-2.50 


2.22 


5.71 


0.44 


17.961 


1.42E-01 


1.29E- 


-02 


5.06E- 


-03 


2.50-3.00 


2.68 


5.65 


0.43 


17.837 


5.01E-02 


7.23E- 


-03 


2.59E- 


-03 


3.00-3.50 


3.26 


5.71 


0.43 


18.108 


1.65E-02 


4.88E- 


-03 


1.48E- 


-03 


3.50-4.00 


3.73 


5.64 


0.43 


17.960 


7.44E-03 


2.65E- 


-03 


2.12E- 


-03 


4.00-5.00 


4.41 


5.74 


0.48 


19.844 


2.03E-03 


9.28E- 


-04 


7.55E- 


-04 


5.00-6.00 


5.44 


5.59 


0.38 


16.530 


3.53E-03 


2.50E- 


-03 


2.75E- 


-03 


6.00-7.00 


6.82 


6.08 


0.53 


22.453 


4.82E-04 


4.82E- 


-04 


9.83E- 


-04 



Table 5.5: (continued) 



p T [GeV] < pr> < M > < x F > < E > Ed 3 a /dp 3 Stat. Error. Syst. Error 









u.zu 




< 7.20 












0.00-0.25 


0.17 


6.63 


0.44 


18.205 


1.25E+00 


8.81E- 


-02 


4.28E- 


-02 


0.25-0.50 


0.38 


6.67 


0.44 


18.246 


1.41E+00 


6.54E- 


-02 


4.28E- 


-02 


0.50-0.75 


0.63 


6.67 


0.43 


18.145 


1.19E-f 


-00 


4.62E- 


-02 


3.33E- 


-02 


0.75-1.00 


0.87 


6.66 


0.44 


18.319 


8.88E- 


-01 


3.66E- 


-02 


2.41E- 


-02 


1.00-1.25 


1.12 


6.68 


0.44 


18.231 


6.87E- 


-01 


2.71E- 


-02 


1.87E- 


-02 


1.25-1.50 


1.37 


6.69 


0.44 


18.247 


4.30E- 


-01 


2.10E- 


-02 


1.20E- 


-02 


1.50-1.75 


1.62 


6.70 


0.43 


18.140 


2.45E- 


-01 


1.56E- 


-02 


7.90E- 


-03 


1.75-2.00 


1.86 


6.70 


0.44 


18.376 


1.59E- 


-01 


1.21E- 


-02 


5.80E- 


-03 


2.00-2.50 


2.20 


6.71 


0.43 


18.284 


7.74E- 


-02 


5.60E- 


-03 


2.67E- 


-03 


2.50-3.00 


2.70 


6.71 


0.43 


18.307 


2.80E- 


-02 


3.09E- 


-03 


1.54E- 


-03 


3.00-3.50 


3.21 


6.72 


0.45 


18.919 


7.96E- 


-03 


1.93E- 


-03 


5.68E- 


-04 


3.50-4.00 


3.67 


6.64 


0.43 


18.472 


4.06E- 


-03 


1.16E- 


-03 


5.21E- 


-04 


4.00-5.00 


4.49 


6.74 


0.44 


18.966 


2.23E- 


-03 


6.52E- 


-04 


4.63E- 


-04 


5.00-6.00 


5.23 


6.71 


0.46 


19.640 


4.46E- 


-04 


4.46E- 


-04 


2.20E- 


-04 



Table 5.5: (continued) 



p T [GeV] < pr> < M > < x F > < E > Ed 3 a /dp 3 Stat. Error. Syst. Error 









7 9D 
i .zu 




- - O. I u 












0.00-0.25 


0.16 


7.80 


0.44 


18.755 


6.65E- 


-01 


5.96E- 


-02 


4.80E- 


-02 


0.25-0.50 


0.39 


7.84 


0.44 


18.642 


8.30E- 


-01 


4.19E- 


-02 


2.80E- 


-02 


0.50-0.75 


0.63 


7.82 


0.44 


18.725 


6.88E- 


-01 


2.99E- 


-02 


2.27E- 


-02 


0.75-1.00 


0.87 


7.86 


0.44 


18.712 


5.48E- 


-01 


2.25E- 


-02 


1.50E- 


-02 


1.00-1.25 


1.12 


7.87 


0.44 


18.827 


3.60E- 


-01 


1.73E- 


-02 


9.79E- 


-03 


1.25-1.50 


1.37 


7.88 


0.44 


18.856 


2.59E- 


-01 


1.30E- 


-02 


7.09E- 


-03 


1.50-1.75 


1.61 


7.85 


0.43 


18.694 


1.59E- 


-01 


9.36E- 


-03 


4.49E- 


-03 


1.75-2.00 


1.87 


7.94 


0.43 


18.630 


1.05E- 


-01 


7.17E- 


-03 


3.19E- 


-03 


2.00-2.50 


2.21 


7.89 


0.44 


18.858 


5.49E- 


-02 


3.57E- 


-03 


1.65E- 


-03 


2.50-3.00 


2.70 


7.95 


0.43 


18.752 


1.66E- 


-02 


1.83E- 


-03 


5.96E- 


-04 


3.00-3.50 


3.21 


7.79 


0.42 


18.363 


5.26E- 


-03 


1.10E- 


-03 


3.61E- 


-04 


3.50-4.00 


3.70 


7.71 


0.44 


19.029 


2.20E- 


-03 


5.98E- 


-04 


2.31E- 


-04 


4.00-5.00 


4.26 


7.76 


0.42 


18.531 


7.84E- 


-04 


3.14E- 


-04 


1.22E- 


-04 


5.00-6.00 


5.51 


7.29 


0.41 


18.390 


1.45E- 


-04 


1.45E- 


-04 


9.02E- 


-05 



Table 5.5: (continued) 



p T [GeV] < pr> < M > < x F > < E > Ed 3 a /dp 3 Stat. Error. Syst. Error 









1 n 
















0.00-0.25 


0.18 


11.71 


0.42 


20.021 


2.17E- 


-02 


1.73E- 


-02 


2.43E- 


-03 


0.25-0.50 


0.38 


11.56 


0.44 


20.773 


5.98E- 


-02 


1.08E- 


-02 


4.66E- 


-03 


0.50-0.75 


0.63 


11.56 


0.43 


20.459 


4.32E- 


-02 


7.20E- 


-03 


3.31E- 


-03 


0.75-1.00 


0.87 


11.52 


0.44 


20.784 


4.28E- 


-02 


6.62E- 


-03 


3.99E- 


-03 


1.00-1.25 


1.12 


11.43 


0.43 


20.407 


2.58E- 


-02 


4.00E- 


-03 


1.22E- 


-03 


1.25-1.50 


1.37 


11.44 


0.44 


20.586 


1.94E- 


-02 


3.39E- 


-03 


9.90E- 


-04 


1.50-1.75 


1.62 


11.49 


0.45 


20.970 


6.01E- 


-03 


2.93E- 


-03 


4.91E- 


-04 


1.75-2.00 


1.87 


11.58 


0.44 


20.665 


1.14E- 


-02 


2.54E- 


-03 


9.49E- 


-04 


2.00-2.50 


2.21 


11.57 


0.43 


20.299 


4.87E- 


-03 


8.33E- 


-04 


3.53E- 


-04 


2.50-3.00 


2.68 


11.54 


0.45 


21.151 


1.90E- 


-03 


4.60E- 


-04 


1.10E- 


-04 


3.00-3.50 


3.18 


11.33 


0.47 


21.718 


5.87E- 


-04 


3.70E- 


-04 


5.84E- 


-05 


3.50-4.00 


3.65 


11.49 


0.41 


19.832 


4.52E- 


-04 


3.39E- 


-04 


6.75E- 


-05 


4.00-5.00 


4.67 


11.81 


0.41 


20.319 


1.10E- 


-04 


1.10E- 


-04 


3.26E- 


-05 


5.00-6.00 


5.64 


11.40 


0.41 


20.415 


5.57E- 


-05 


5.57E- 


-05 


1.13E- 


-04 



Table 5.5: (continued) 



p T [GeV] < pr> < M > < x F > < E > Ed 3 a /dp 3 Stat. Error. Syst. Error 















12. 


85 


< M„+„- 


< 16.85 















00- 





25 


0.06 


13.22 





49 


23.133 


5.19E- 


-02 


5.19E- 


-02 


9.21E- 


-03 





25- 





50 


0.41 


13.25 





51 


23.884 


7.28E- 


-02 


5.15E- 


-02 


8.29E- 


-03 





50- 





75 


0.67 


13.37 





42 


21.241 


1.61E- 


-02 


5.47E- 


-03 


9.73E- 


-04 





75- 


1 


00 


0.84 


13.68 





48 


23.158 


1.09E- 


-02 


4.43E- 


-03 


6.29E- 


-04 


1 


00- 


1 


25 


1.08 


13.51 





48 


23.098 


9.68E- 


-03 


4.26E- 


-03 


6.40E- 


-04 


1 


25- 


1 


50 


1.40 


14.02 





45 


22.405 


4.22E- 


-03 


2.15E- 


-03 


2.87E- 


-04 


1 


50- 


1 


75 


1.64 


13.39 





43 


21.321 


5.02E- 


-03 


2.31E- 


-03 


3.69E- 


-04 


1 


75- 


2 


00 


1.95 


13.36 





48 


23.111 


1.61E- 


-02 


1.14E- 


-02 


1.96E- 


-03 


2 


00- 


2 


50 


2.09 


13.24 





40 


20.624 


6.51E- 


-03 


4.61E- 


-03 


9.44E- 


-04 


2 


50- 


3 


00 


2.77 


13.56 





52 


24.424 


1.53E- 


-03 


1.53E- 


-03 


3.07E- 


-04 


3 


00- 


3 


50 


3.17 


13.43 





42 


21.367 


5.40E- 


-03 


5.40E- 


-03 


3.18E- 


-03 



Table 5.6: Invariant form Ed 3 a/dp 3 (in units of pb / GeV 2 / nucleon) for the hydrogen 
cross section over the range 0.55 < xp < 0.8. Statistical and point-to-point systematic 
uncertainties are shown. An additional ±6.5% uncertainty in the normalization is common 
to all data points. 



p T [GeV] < p T > < M > < x F > < E > Ed 3 a /dp 3 Stat. Error. Syst. Error 



4.20 < M„+„- < 5.20 






00- 





25 





17 


4 


87 





62 


24.586 


1 


86E+00 


3 


00E- 


-01 


2 


19E- 


-01 





25- 





50 





39 


4 


78 





63 


24.948 


2 


61E+00 


2 


48E- 


-01 


1 


54E- 


-01 





50- 





75 





63 


4 


80 





62 


24.641 


1 


71E+00 


1 


59E- 


-01 


8 


19E- 


-02 





75- 


1 


00 





87 


4 


79 





63 


24.800 


1 


35E+00 


1 


18E- 


-01 


5 


93E- 


-02 


1 


00- 


1 


25 


1 


12 


4 


78 





61 


24.314 


9 


03E-01 


8 


65E- 


-02 


3 


87E- 


-02 


1 


25- 


1 


50 


1 


36 


4 


77 





62 


24.673 


4 


90E-01 


6 


05E- 


-02 


2 


06E- 


-02 


1 


50- 


1 


75 


1 


61 


4 


81 





62 


24.731 


2 


00E-01 


3 


54E- 


-02 


9 


80E- 


-03 


1 


75- 


2 


00 


1 


86 


4 


72 





60 


23.822 


1 


03E-01 


3 


21E- 


-02 


7 


25E- 


-03 


2 


00- 


2 


50 


2 


20 


4 


77 





62 


24.448 


6 


73E-02 


1 


40E- 


-02 


3 


81E- 


-03 


2 


50- 


3 


00 


2 


66 


4 


70 





62 


24.646 


1 


60E-02 


9 


68E- 


-03 


1 


57E- 


-03 


3 


00- 


3 


50 


3 


30 


4 


93 





67 


26.580 


1 


57E-02 


9 


58E- 


-03 


2 


52E- 


-03 


3 


50- 


4 


00 


3 


67 


4 


88 





63 


25.254 


1 


64E-02 


1 


16E- 


-02 


1 


17E- 


-02 



Table 5.6: (continued) 



p T [GeV] < pr> < M > < x F > < E > Ed 3 a /dp 3 Stat. Error. Syst. Error 













<r f\ 9D 
- - u.z,u 












00-0 25 


0.17 


5.67 


0.63 


25.027 


8.73E- 


-01 


9.90E- 


-02 


5.51E- 


-02 


0.25-0.50 


0.38 


5.69 


0.62 


24.668 


9.47E- 


-01 


6.90E- 


-02 


4.34E- 


-02 


0.50-0.75 


0.63 


5.71 


0.62 


24.597 


7.51E- 


-01 


4.52E- 


-02 


3.21E- 


-02 


0.75-1.00 


0.87 


5.70 


0.62 


24.795 


5.30E- 


-01 


3.87E- 


-02 


1.95E- 


-02 


1.00-1.25 


1.12 


5.73 


0.62 


24.646 


3.41E- 


-01 


2.75E- 


-02 


1.15E- 


-02 


1.25-1.50 


1.36 


5.71 


0.62 


24.700 


2.30E- 


-01 


2.03E- 


-02 


8.16E- 


-03 


1.50-1.75 


1.61 


5.75 


0.63 


25.074 


1.08E- 


-01 


1.43E- 


-02 


5.43E- 


-03 


1.75-2.00 


1.88 


5.69 


0.61 


24.352 


5.90E- 


-02 


8.88E- 


-03 


3.68E- 


-03 


2.00-2.50 


2.17 


5.72 


0.61 


24.255 


2.55E- 


-02 


3.88E- 


-03 


1.12E- 


-03 


2.50-3.00 


2.73 


5.69 


0.61 


24.371 


5.26E- 


-03 


1.94E- 


-03 


7.89E- 


-04 


3.00-3.50 


3.33 


5.85 


0.66 


26.333 


4.22E- 


-03 


1.57E- 


-03 


7.21E- 


-04 


3.50-4.00 


3.74 


5.40 


0.64 


25.676 


3.07E- 


-03 


2.18E- 


-03 


8.02E- 


-04 


4.00-5.00 


4.20 


5.83 


0.71 


28.484 


1.50E- 


-03 


1.50E- 


-03 


5.70E- 


-04 










6.20 


< AW 


< 7.20 












0.00-0.25 


0.16 


6.65 


0.62 


24.963 


3.38E- 


-01 


4.39E- 


-02 


2.69E- 


-02 



Table 5.6: (continued) 



p T [GeV] < pr> < M > < x F > < E > Ed 3 a /dp 3 Stat. Error. Syst. Error 



0.25-0.50 


0.39 


6.64 


0.63 


25.195 


3.69E- 


-01 


3.60E- 


-02 


1.48E- 


-02 


0.50-0.75 


0.63 


6.67 


0.62 


24.953 


3.24E- 


-01 


2.41E- 


-02 


1.09E- 


-02 


0.75-1.00 


0.87 


6.65 


0.62 


24.933 


2.79E- 


-01 


1.82E- 


-02 


9.13E- 


-03 


1.00-1.25 


1.12 


6.68 


0.62 


25.010 


1.65E- 


-01 


1.31E- 


-02 


6.72E- 


-03 


1.25-1.50 


1.37 


6.67 


0.62 


24.847 


8.91E- 


-02 


1.01E- 


-02 


4.80E- 


-03 


1.50-1.75 


1.61 


6.70 


0.61 


24.564 


5.63E- 


-02 


6.72E- 


-03 


2.32E- 


-03 


1.75-2.00 


1.86 


6.64 


0.62 


24.857 


3.50E- 


-02 


4.69E- 


-03 


1.45E- 


-03 


2.00-2.50 


2.20 


6.70 


0.61 


24.819 


8.22E- 


-03 


1.77E- 


-03 


4.07E- 


-04 


2.50-3.00 


2.70 


6.64 


0.61 


24.863 


4.86E- 


-03 


1.53E- 


-03 


2.53E- 


-04 


3.00-3.50 


3.23 


6.73 


0.61 


24.780 


1.87E- 


-03 


7.67E- 


-04 


1.44E- 


-04 


3.50-4.00 


3.62 


6.76 


0.60 


24.599 


9.29E- 


-04 


6.17E- 


-04 


3.42E- 


-04 


4.00-5.00 


4.29 


6.85 


0.59 


24.152 


3.88E- 


-04 


3.07E- 


-04 


1.08E- 


-04 



7.20 < M„+„- < 8.70 



0.00-0.25 


0.18 


7.87 


0.25-0.50 


0.38 


7.79 


0.50-0.75 


0.64 


7.83 


0.75-1.00 


0.88 


7.85 


1.00-1.25 


1.12 


7.82 



0.63 25.571 2.20E-01 

0.62 25.230 2.34E-01 

0.62 25.198 1.74E-01 

0.61 25.090 1.46E-01 

0.62 25.364 9.27E-02 



3.02E-02 1.18E-02 

2.10E-02 8.44E-03 

1.50E-02 5.38E-03 

1.11E-02 4.86E-03 

8.34E-03 3.10E-03 



Table 5.6: (continued) 



p T [GeV] < pr> < M > < x F > < E > Ed 3 a /dp 3 Stat. Error. Syst. Error 



1.25-1.50 


1.37 


7.91 


0.61 


25.161 


5.32E- 


-02 


5.74E- 


-03 


1.83E- 


-03 


1.50-1.75 


1.61 


7.81 


0.62 


25.266 


3.45E- 


-02 


4.14E- 


-03 


1.40E- 


-03 


1.75-2.00 


1.86 


7.82 


0.61 


25.098 


2.41E- 


-02 


3.17E- 


-03 


9.33E- 


-04 


2.00-2.50 


2.22 


7.82 


0.62 


25.486 


4.76E- 


-03 


1.01E- 


-03 


2.28E- 


-04 


2.50-3.00 


2.72 


7.78 


0.62 


25.371 


2.11E- 


-03 


7.78E- 


-04 


8.39E- 


-05 


3.00-3.50 


3.28 


8.01 


0.62 


25.553 


1.64E- 


-03 


5.08E- 


-04 


9.57E- 


-05 


3.50-4.00 


3.72 


7.83 


0.61 


25.046 


9.17E- 


-04 


4.58E- 


-04 


1.03E- 


-04 


4.00-5.00 


4.31 


7.92 


0.60 


24.851 


2.62E- 


-04 


4.27E- 


-04 


1.50E- 


-04 


5.00-6.00 


5.51 


7.44 


0.58 


24.190 


2.08E- 


-04 


2.08E- 


-04 


9.57E- 


-05 



10.85 < M M+M - < 12.85 



0.00-0.25 


0.21 


11.20 


0.62 


26.684 


8.62E- 


-03 


4.34E- 


-03 


7.43E- 


-04 


0.25-0.50 


0.38 


11.33 


0.59 


25.404 


8.13E- 


-03 


2.76E- 


-03 


3.99E- 


-04 


0.50-0.75 


0.64 


11.50 


0.60 


26.056 


9.56E- 


-03 


2.71E- 


-03 


4.70E- 


-04 


0.75-1.00 


0.87 


11.43 


0.63 


26.839 


9.35E- 


-03 


2.61E- 


-03 


1.67E- 


-03 


1.00-1.25 


1.10 


11.38 


0.63 


26.864 


3.26E- 


-03 


1.41E- 


-03 


1.54E- 


-04 


1.25-1.50 


1.35 


11.76 


0.63 


27.083 


3.75E- 


-03 


1.19E- 


-03 


2.38E- 


-04 


1.50-1.75 


1.57 


11.49 


0.61 


26.375 


3.71E- 


-03 


1.12E- 


-03 


2.21E- 


-04 


1.75-2.00 


1.86 


11.30 


0.63 


26.843 


1.58E- 


-03 


5.06E- 


-04 


1.09E- 


-04 



Table 5.6: (continued) 



p T [GeV] < pr> < M > < x F > < E > Ed 3 a /dp 3 Stat. Error. Syst. Error 



2.00-2.50 


2.31 


11.49 


0.64 


27.573 


4.19E- 


-04 


2.21E- 


-04 


2.37E- 


■04 


2.50-3.00 


2.59 


11.02 


0.60 


25.991 


3.00E- 


-04 


3.00E- 


-04 


3.60E- 


-05 


3.00-3.50 


3.14 


11.24 


0.61 


26.422 


6.62E- 


-04 


6.62E- 


-04 


1.07E- 


-04 


3.50-4.00 


3.60 


10.95 


0.55 


24.276 


3.14E- 


-05 


3.14E- 


-05 


3.82E- 


-06 









12.85 


< M M+M - 


< 16.85 












0.00-0.25 


0.20 


13.03 


0.60 


26.813 


1.97E- 


-02 


1.97E- 


-02 


1.27E- 


-02 


0.25-0.50 


0.38 


13.73 


0.67 


29.495 


7.02E- 


■03 


4.11E- 


-03 


8.60E- 


-04 


0.50-0.75 


0.66 


13.23 


0.61 


27.224 


6.90E- 


-03 


3.45E- 


-03 


1.06E- 


-03 


0.75-1.00 


0.88 


13.89 


0.63 


28.152 


1.69E- 


-03 


1.16E- 


-03 


1.52E- 


-04 


1.00-1.25 


1.12 


13.55 


0.58 


26.315 


1.37E- 


-03 


1.09E- 


-03 


1.49E- 


-04 


1.25-1.50 


1.33 


13.01 


0.64 


28.196 


9.64E- 


-04 


1.13E- 


-03 


1.12E- 


-04 


1.50-1.75 


1.61 


13.02 


0.56 


25.240 


5.00E- 


-04 


5.00E- 


-04 


6.03E- 


-05 


2.00-2.50 


2.24 


13.28 


0.71 


30.505 


3.34E- 


-04 


3.34E- 


-04 


4.89E- 


-05 



Table 5.7: Invariant form Ed 3 a /dp 3 (in units of pb / GeV 2 / nucleon) for the deuterium 
cross section over the range -0.05 < xf < 0.15. Statistical and point-to-point systematic 
uncertainties are shown. An additional ±6.5% uncertainty in the normalization is common 
to all data points. 



p T [GeV] < p T > < M > < x F > < E > Ed 3 a /dp 3 Stat. Error. Syst. Error 



4.20 < M„+„- < 5.20 






00- 





25 





16 


4 


66 





11 


6.456 


1.35E+01 


1.03E+00 


8.93E- 


-01 





25- 





50 





38 


4 


68 





12 


6.569 


1.44E+01 


7.84E- 


-01 


4.32E- 


-01 





50- 





75 





62 


4 


66 





12 


6.588 


1.25E+01 


5.79E- 


-01 


3.76E- 


-01 





75- 


1 


00 





87 


4 


65 





12 


6.644 


9.01E+00 


4.29E- 


-01 


2.85E- 


-01 


1 


00- 


1 


25 


1 


12 


4 


64 





12 


6.694 


6.82E+00 


3.81E- 


-01 


2.15E- 


-01 


1 


25- 


1 


50 


1 


37 


4 


63 





12 


6.844 


5.05E+00 


3.24E- 


-01 


1.64E- 


-01 


1 


50- 


1 


75 


1 


63 


4 


66 





12 


6.771 


2.16E+00 


2.33E- 


-01 


9.04E- 


-02 


1 


75- 


2 


00 


1 


87 


4 


61 





12 


6.826 


1.74E+00 


2.22E- 


-01 


7.72E- 


-02 


2 


00- 


2 


50 


2 


22 


4 


59 





12 


6.956 


9.53E-01 


1.10E- 


-01 


4.21E- 


-02 


2 


50- 


-3 


00 


2 


65 


4 


67 





12 


7.068 


1.91E-01 


5.35E- 


-02 


1.36E- 


-02 


3 


50- 


4 


00 


3 


70 


4 


60 





12 


7.442 


2.20E-01 


1.10E- 


-01 


6.12E- 


-02 



Table 5.7: (continued) 



p T [GeV] < pr> < M > < x F > < E > Ed 3 a /dp 3 Stat. Error. Syst. Error 









5.20 


< M„+„- 


< 6.20 










0.00-0.25 


0.16 


5.55 


0.11 


7.116 


4.36E+00 


6.50E- 


-01 


2.97E- 


-01 


0.25-0.50 


0.39 


5.61 


0.11 


7.188 


4.28E+00 


3.69E- 


-01 


1.75E- 


-01 


0.50-0.75 


0.63 


5.60 


0.12 


7.249 


4.11E+00 


2.90E- 


-01 


1.54E- 


-01 


0.75-1.00 


0.87 


5.62 


0.12 


7.395 


3.11E+00 


2.47E- 


-01 


1.15E- 


-01 


1.00-1.25 


1.12 


5.65 


0.12 


7.384 


1.63E+00 


1.66E- 


-01 


6.78E- 


-02 


1.25-1.50 


1.37 


5.59 


0.12 


7.537 


1.23E+00 


1.63E- 


-01 


6.56E- 


-02 


1.50-1.75 


1.62 


5.65 


0.12 


7.527 


9.74E-01 


1.31E- 


-01 


8.31E- 


-02 


1.75-2.00 


1.86 


5.61 


0.12 


7.547 


7.12E-01 


1.17E- 


-01 


8.39E- 


-02 


2.00-2.50 


2.25 


5.63 


0.12 


7.644 


4.42E-01 


6.96E- 


-02 


5.71E- 


-02 


2.50-3.00 


2.61 


5.78 


0.13 


8.023 


1.42E-01 


4.76E- 


-02 


3.26E- 


-02 


3.00-3.50 


3.15 


5.54 


0.13 


8.066 


3.99E-02 


2.36E- 


-02 


1.30E- 


-02 



6.20 < M M+ ^- < 7.20 

0.00-0.25 0.18 6.69 0.10 7.819 1.31E+00 2.28E-01 1.08E-01 

0.25-0.50 0.38 6.69 0.11 7.981 1.79E+00 1.91E-01 1.19E-01 

0.50-0.75 0.63 6.74 0.10 7.960 1.34E+00 1.49E-01 7.59E-02 



Table 5.7: (continued) 



p T [GeV] < pr> < M > < x F > < E > Ed 3 a /dp 3 Stat. Error. Syst. Error 



0.75-1.00 


0.87 


6.70 


0.10 


7.999 


1.31E-f 


-00 


1.18E- 


-01 


7.05E- 


-02 


1.00-1.25 


1.11 


6.71 


0.10 


8.011 


6.17E- 


-01 


8.45E- 


-02 


6.31E- 


-02 


1.25-1.50 


1.37 


6.68 


0.10 


8.089 


5.56E- 


-01 


8.99E- 


-02 


3.65E- 


-02 


1.50-1.75 


1.62 


6.82 


0.12 


8.464 


3.18E- 


-01 


6.52E- 


-02 


4.46E- 


-02 


1.75-2.00 


1.88 


6.70 


0.10 


8.069 


1.89E- 


-01 


4.09E- 


-02 


1.50E- 


-02 


2.00-2.50 


2.13 


6.66 


0.09 


8.017 


9.92E- 


-02 


3.31E- 


-02 


9.85E- 


-03 


2.50-3.00 


2.69 


6.91 


0.12 


8.873 


4.07E- 


-02 


1.59E- 


-02 


7.55E- 


-03 


3.00-3.50 


3.21 


6.91 


0.12 


8.902 


4.33E- 


-02 


2.20E- 


-02 


1.74E- 


-02 



7.20 < M„+„- < 8.70 



0.00-0.25 


0.16 


8.10 


0.09 


8.987 


8.63E- 


-01 


8.43E- 


-02 


5.27E- 


-02 


0.25-0.50 


0.39 


8.06 


0.09 


9.003 


9.77E- 


-01 


5.96E- 


-02 


4.12E- 


-02 


0.50-0.75 


0.63 


8.05 


0.09 


9.015 


9.64E- 


-01 


4.71E- 


-02 


2.86E- 


-02 


0.75-1.00 


0.87 


8.08 


0.10 


9.074 


7.73E- 


-01 


3.58E- 


-02 


2.51E- 


-02 


1.00-1.25 


1.12 


8.08 


0.10 


9.162 


5.83E- 


-01 


3.06E- 


-02 


2.30E- 


-02 


1.25-1.50 


1.37 


8.11 


0.10 


9.275 


3.60E- 


-01 


2.23E- 


-02 


1.46E- 


-02 


1.50-1.75 


1.61 


8.06 


0.10 


9.283 


2.54E- 


-01 


1.90E- 


-02 


1.09E- 


-02 


1.75-2.00 


1.86 


8.07 


0.11 


9.369 


1.43E- 


-01 


1.51E- 


-02 


7.08E- 


-03 


2.00-2.50 


2.23 


8.13 


0.11 


9.558 


8.54E- 


-02 


8.31E- 


-03 


3.24E- 


-03 



Table 5.7: (continued) 



p T [GeV] < pr> < M > < x F > < E > Ed 3 a /dp 3 Stat. Error. Syst. Error 



2.50-3.00 


2.69 


7.96 


0.10 


9.393 


3.09E- 


-02 


5.21E- 


-03 


2.22E- 


-03 


3.00-3.50 


3.21 


8.10 


0.11 


9.708 


1.24E- 


-02 


3.19E- 


-03 


1.91E- 


-03 


3.50-4.00 


3.74 


7.99 


0.12 


10.020 


6.94E- 


-03 


3.25E- 


-03 


1.30E- 


-03 


4.00-5.00 


4.25 


8.57 


0.10 


10.455 


2.11E- 


-03 


1.22E- 


-03 


5.46E- 


-04 










10.85 


< M^- 


< 12.85 












0.00-0.25 


0.19 


11.13 


0.04 


11.595 


4.74E- 


-02 


2.94E- 


-02 


9.09E- 


-03 


0.25-0.50 


0.41 


11.32 


0.10 


12.003 


9.01E- 


-02 


2.85E- 


-02 


1.92E- 


-02 


0.50-0.75 


0.63 


11.13 


0.10 


11.857 


1.01E- 


-01 


2.45E- 


-02 


3.02E- 


-02 


0.75-1.00 


0.87 


11.31 


0.11 


12.136 


8.56E- 


-02 


2.53E- 


-02 


1.57E- 


-02 


1.00-1.25 


1.12 


11.15 


0.11 


12.091 


4.04E- 


-02 


1.36E- 


-02 


4.23E- 


-03 


1.25-1.50 


1.34 


11.49 


0.11 


12.402 


2.88E- 


-02 


1.13E- 


-02 


3.67E- 


-03 


1.50-1.75 


1.61 


11.23 


0.12 


12.307 


3.24E- 


-02 


1.28E- 


-02 


4.91E- 


-03 


1.75-2.00 


1.87 


11.27 


0.12 


12.370 


3.27E- 


-02 


1.37E- 


-02 


1.02E- 


-02 


2.00-2.50 


2.26 


11.65 


0.13 


12.853 


1.22E- 


-02 


6.57E- 


-03 


1.04E- 


-02 


2.50-3.00 


2.71 


11.33 


0.14 


12.882 


6.74E- 


-02 


4.77E- 


-02 


8.51E- 


-02 


3.00-3.50 


3.49 


11.13 


0.12 


12.537 


3.48E- 


-03 


3.48E- 


-03 


2.55E- 


-03 



Table 5.8: Invariant form Ed 3 a/dp 3 (in units of pb / GeV 2 / nucleon) for the deuterium 
cross section over the range 0.15 < xp < 0.35. Statistical and point-to-point systematic 
uncertainties are shown. An additional ±6.5% uncertainty in the normalization is common 
to all data points. 



p T [GeV] < p T > < M > < x F > < E > Ed 3 a /dp 3 Stat. Error. Syst. Error 



4.20 < M„+„- < 5.20 






00- 





25 





17 


4 


65 





25 


10.788 


1 


27E+01 


4 


92E- 


-01 


4 


75E- 


-01 





25- 





50 





39 


4 


66 





25 


10.855 


1 


42E+01 


3 


41E- 


-01 


4 


20E- 


-01 





50- 





75 





63 


4 


66 





25 


10.844 


1 


24E+01 


2 


58E- 


-01 


3 


41E- 


-01 





75- 


1 


00 





87 


4 


66 





25 


10.844 


8 


91E+00 


1 


90E- 


-01 


2 


41E- 


-01 


1 


00- 


1 


25 


1 


12 


4 


65 





25 


10.813 


6 


51E+00 


1 


47E- 


-01 


1 


76E- 


-01 


1 


25- 


1 


50 


1 


37 


4 


65 





25 


10.752 


4 


28E+00 


1 


17E- 


-01 


1 


20E- 


-01 


1 


50- 


1 


75 


1 


61 


4 


64 





25 


10.856 


2 


87E+00 


9 


27E- 


-02 


8 


09E- 


-02 


1 


75- 


2 


00 


1 


87 


4 


64 





25 


10.903 


1 


81E+00 


7 


28E- 


-02 


5 


54E- 


-02 


2 


00- 


2 


50 


2 


21 


4 


66 





26 


11.229 


8 


43E-01 


3 


86E- 


-02 


2 


84E- 


-02 


2 


50- 


3 


00 


2 


71 


4 


65 





26 


11.449 


3 


39E-01 


2 


42E- 


-02 


1 


96E- 


-02 


3 


00- 


3 


50 


3 


21 


4 


63 





26 


11.694 


1 


01E-01 


1 


76E- 


-02 


1 


31E- 


-02 


3 


50- 


4 


00 


3 


75 


4 


61 





28 


12.354 


4 


07E-02 


1 


27E- 


-02 


7 


20E- 


-03 


4 


00- 


5 


00 


4 


42 


4 


57 





28 


12.559 


1 


69E-02 


5 


90E- 


-03 


2 


56E- 


-03 


5 


00- 


6 


00 


5 


38 


4 


78 





27 


12.854 


5 


47E-03 


5 


47E- 


-03 


2 


55E- 


-03 



Table 5.8: (continued) 



p T [GeV] < pr> < M > < x F > < E > Ed 3 a /dp 3 Stat. Error. Syst. Error 









o.zu 














0.00-0.25 


0.16 


5.66 


0.27 


12.052 


4.85E+00 


2.24E- 


-01 


2.13E- 


-01 


0.25-0.50 


0.39 


5.67 


0.27 


12.021 


5.95E+00 


1.57E- 


-01 


2.00E- 


-01 


0.50-0.75 


0.63 


5.68 


0.27 


12.058 


4.86E+00 


1.13E- 


-01 


1.53E- 


-01 


0.75-1.00 


0.87 


5.68 


0.27 


11.968 


3.70E+00 


8.79E- 


-02 


1.07E- 


-01 


1.00-1.25 


1.12 


5.69 


0.27 


12.061 


2.81E+00 


7.05E- 


-02 


7.91E- 


-02 


1.25-1.50 


1.37 


5.70 


0.27 


12.061 


1.67E+00 


5.41E- 


-02 


4.85E- 


-02 


1.50-1.75 


1.61 


5.70 


0.27 


12.075 


1.11E+00 


4.14E- 


-02 


3.31E- 


-02 


1.75-2.00 


1.86 


5.66 


0.28 


12.291 


6.86E-01 


3.38E- 


-02 


2.22E- 


-02 


2.00-2.50 


2.21 


5.70 


0.27 


12.328 


3.26E-01 


1.65E- 


-02 


1.21E- 


-02 


2.50-3.00 


2.72 


5.67 


0.27 


12.371 


1.16E-01 


1.08E- 


-02 


5.51E- 


-03 


3.00-3.50 


3.21 


5.60 


0.27 


12.230 


3.27E-02 


5.34E- 


-03 


7.88E- 


-03 


3.50-4.00 


3.72 


5.75 


0.29 


13.320 


3.02E-02 


5.59E- 


-03 


8.07E- 


-03 


4.00-5.00 


4.37 


5.62 


0.31 


13.879 


5.33E-03 


2.33E- 


-03 


6.22E- 


-03 


5.00-6.00 


5.36 


5.61 


0.30 


13.873 


7.66E-03 


7.66E- 


-03 


3.92E- 


-03 



Table 5.8: (continued) 



p T [GeV] < pr> < M > < x F > < E > Ed 3 a /dp 3 Stat. Error. Syst. Error 









u.zu 




<r 7 9D 










0.00-0.25 


0.17 


6.68 


0.27 


12.585 


1.94E+00 


9.54E- 


-02 


1.30E- 


-01 


0.25-0.50 


0.39 


6.69 


0.27 


12.582 


2.43E+00 


6.90E- 


-02 


1.01E- 


-01 


0.50-0.75 


0.63 


6.69 


0.28 


12.674 


1.98E+00 


4.98E- 


-02 


8.25E- 


-02 


0.75-1.00 


0.87 


6.70 


0.27 


12.606 


1.61E+00 


3.93E- 


-02 


6.81E- 


-02 


1.00-1.25 


1.12 


6.70 


0.27 


12.637 


1.15E+00 


3.03E- 


-02 


4.01E- 


-02 


1.25-1.50 


1.36 


6.71 


0.27 


12.660 


7.19E-01 


2.30E- 


-02 


2.10E- 


-02 


1.50-1.75 


1.62 


6.72 


0.27 


12.675 


4.95E-01 


1.83E- 


-02 


2.03E- 


-02 


1.75-2.00 


1.87 


6.72 


0.28 


12.793 


2.94E-01 


1.42E- 


-02 


1.57E- 


-02 


2.00-2.50 


2.21 


6.72 


0.28 


12.934 


1.66E-01 


7.35E- 


-03 


5.23E- 


-03 


2.50-3.00 


2.70 


6.74 


0.27 


12.947 


6.35E-02 


4.49E- 


-03 


2.65E- 


-03 


3.00-3.50 


3.20 


6.66 


0.27 


12.911 


2.66E-02 


2.98E- 


-03 


1.69E- 


-03 


3.50-4.00 


3.73 


6.82 


0.27 


13.075 


4.84E-03 


1.37E- 


-03 


8.00E- 


-04 


4.00-5.00 


4.23 


6.73 


0.29 


13.745 


2.99E-03 


9.76E- 


-04 


4.94E- 


-04 


5.00-6.00 


5.06 


7.09 


0.35 


16.081 


1.34E-03 


1.45E- 


-03 


3.15E- 


-03 



Table 5.8: (continued) 



p T [GeV] < pr> < M > < x F > < E > Ed 3 a /dp 3 Stat. Error. Syst. Error 









7 9D 

1 .ZU 




< 8.70 












0.00-0.25 


0.16 


7.90 


0.26 


12.756 


1.07E+00 


5.25E- 


-02 


4.67E- 


-02 


0.25-0.50 


0.39 


7.89 


0.26 


12.777 


1.23E+00 


3.59E- 


-02 


4.29E- 


-02 


0.50-0.75 


0.63 


7.87 


0.26 


12.805 


1.12E-f 


-00 


2.71E- 


-02 


3.21E- 


-02 


0.75-1.00 


0.87 


7.87 


0.26 


12.845 


8.73E- 


-01 


2.03E- 


-02 


3.59E- 


-02 


1.00-1.25 


1.12 


7.87 


0.26 


12.899 


6.09E- 


-01 


1.53E- 


-02 


2.23E- 


-02 


1.25-1.50 


1.37 


7.89 


0.26 


12.933 


4.30E- 


-01 


1.19E- 


-02 


1.19E- 


-02 


1.50-1.75 


1.62 


7.90 


0.26 


13.020 


2.88E- 


-01 


9.13E- 


-03 


8.20E- 


-03 


1.75-2.00 


1.87 


7.93 


0.26 


13.036 


1.90E- 


-01 


7.04E- 


-03 


5.61E- 


-03 


2.00-2.50 


2.21 


7.89 


0.26 


13.051 


9.25E- 


-02 


3.44E- 


-03 


3.91E- 


-03 


2.50-3.00 


2.72 


7.93 


0.26 


13.186 


4.03E- 


-02 


2.28E- 


-03 


3.00E- 


-03 


3.00-3.50 


3.21 


7.94 


0.26 


13.457 


1.59E- 


-02 


1.43E- 


-03 


1.31E- 


-03 


3.50-4.00 


3.68 


8.22 


0.25 


13.366 


3.37E- 


-03 


8.14E- 


-04 


2.26E- 


-04 


4.00-5.00 


4.45 


7.87 


0.28 


14.118 


2.59E- 


-03 


5.15E- 


-04 


7.51E- 


-04 


5.00-6.00 


5.30 


7.91 


0.27 


14.123 


8.93E- 


-04 


4.65E- 


-04 


7.22E- 


-04 



Table 5.8: (continued) 



p T [GeV] < pr> < M > < x F > < E > Ed 3 a /dp 3 Stat. Error. Syst. Error 













1 n rz> 

1U.OO 
















0.00- 


-0.25 


0.15 


11 


49 


0.26 


15.472 


6.21E- 


-02 


1.40E- 


-02 


8.00E- 


-03 


0.25- 


-0.50 


0.38 


11 


42 


0.27 


15.690 


1.12E- 


-01 


1.33E- 


-02 


1.01E- 


-02 


0.50- 


-0.75 


0.64 


11 


39 


0.27 


15.487 


7.81E- 


-02 


9.09E- 


-03 


4.52E- 


-03 


0.75- 


-1.00 


0.87 


11 


44 


0.27 


15.719 


7.01E- 


-02 


6.93E- 


-03 


5.63E- 


-03 


1.00- 


-1.25 


1.12 


11 


41 


0.28 


15.759 


5.44E- 


-02 


5.24E- 


-03 


6.89E- 


-03 


1.25- 


-1.50 


1.38 


11 


41 


0.26 


15.393 


3.13E- 


-02 


3.72E- 


-03 


2.87E- 


-03 


1.50- 


-1.75 


1.62 


11 


39 


0.27 


15.729 


2.19E- 


-02 


3.34E- 


-03 


3.08E- 


-03 


1.75- 


-2.00 


1.87 


11 


47 


0.29 


16.177 


1.57E- 


-02 


2.63E- 


-03 


3.28E- 


-03 


2.00- 


-2.50 


2.22 


11 


42 


0.26 


15.592 


7.04E- 


-03 


1.01E- 


-03 


1.14E- 


-03 


2.50- 


-3.00 


2.73 


11 


25 


0.25 


15.216 


2.06E- 


-03 


5.37E- 


-04 


5.23E- 


-04 


3.00- 


-3.50 


3.20 


11 


35 


0.25 


15.369 


9.92E- 


-04 


4.02E- 


-04 


2.63E- 


-04 


3.50- 


-4.00 


3.66 


11 


37 


0.31 


17.000 


5.29E- 


-04 


3.12E- 


-04 


9.67E- 


-05 


4.00- 


-5.00 


4.21 


11 


17 


0.33 


17.418 


2.25E- 


-04 


2.25E- 


-04 


2.20E- 


-04 



12.85 < M M+M - < 16.85 
0.00-0.25 0.21 13.63 0.26 17.097 2.41E-02 2.23E-02 9.04E-03 



Table 5.8: (continued) 



p T [GeV] < pr> < M > < x F > < E > Ed 3 a /dp 3 Stat. Error. Syst. Error 



0.25- 


-0.50 


0.34 


13.16 


0.28 


17.161 


1.50E- 


-02 


1.36E- 


-02 


4.06E- 


-03 


0.50- 


-0.75 


0.64 


13.10 


0.28 


17.006 


2.33E- 


-02 


2.77E- 


-02 


3.15E- 


-03 


0.75- 


-1.00 


0.88 


13.48 


0.31 


18.132 


4.91E- 


-02 


1.49E- 


-02 


1.35E- 


-03 


1.00- 


-1.25 


1.10 


13.58 


0.24 


16.468 


1.09E- 


-02 


5.88E- 


-03 


2.04E- 


-03 


1.25- 


-1.50 


1.34 


13.49 


0.28 


17.375 


2.08E- 


-02 


8.25E- 


-03 


2.68E- 


-03 


1.50- 


-1.75 


1.58 


13.31 


0.30 


17.868 


6.07E- 


-03 


3.36E- 


-03 


1.73E- 


-03 


1.75- 


-2.00 


1.86 


13.45 


0.30 


17.818 


3.73E- 


-02 


2.64E- 


-02 


1.30E- 


-02 


2.00- 


-2.50 


2.45 


13.25 


0.35 


19.118 


2.86E- 


-03 


2.86E- 


-03 


2.60E- 


-04 


2.50- 


-3.00 


2.95 


13.08 


0.27 


17.073 


2.02E- 


-02 


2.02E- 


-02 


5.17E- 


-04 



Table 5.9: Invariant form Ed 3 a/dp 3 (in units of pb / GeV 2 / nucleon) for the deuterium 
cross section over the range 0.35 < xp < 0.55. Statistical and point-to-point systematic 
uncertainties are shown. An additional ±6.5% uncertainty in the normalization is common 
to all data points. 



p T [GeV] < p T > < M > < x F > < E > Ed 3 a /dp 3 Stat. Error. Syst. Error 

















4.20 


< m m+m - 


< 


5.2 

















00-0 


25 





16 


4 


74 


0.43 


17.465 


8 


63E+00 


4 


41E- 


-01 


6 


65E- 


-01 





25-0 


50 





38 


4 


75 


0.43 


17.338 


9 


61E+00 


3 


03E- 


-01 


3 


35E- 


-01 





50-0 


75 





63 


4 


74 


0.43 


17.430 


7 


31E+00 


2 


04E- 


-01 


2 


13E- 


-01 





75-1 


00 





87 


4 


73 


0.43 


17.416 


5 


68E+00 


1 


60E- 


-01 


1 


75E- 


-01 


1 


00-1 


25 


1 


12 


4 


74 


0.43 


17.371 


3 


67E+00 


1 


16E- 


-01 


1 


64E- 


-01 


1 


25-1 


50 


1 


37 


4 


73 


0.43 


17.455 


2 


33E+00 


8 


90E- 


-02 


1 


11E- 


-01 


1 


50-1 


75 


1 


61 


4 


72 


0.43 


17.482 


1 


28E+00 


6 


67E- 


-02 


4 


15E- 


-02 


1 


75-2 


00 


1 


87 


4 


74 


0.43 


17.396 


7 


94E-01 


4 


85E- 


-02 


2 


77E- 


-02 


2 


00-2 


50 


2 


20 


4 


72 


0.43 


17.313 


3 


12E-01 


2 


15E- 


-02 


1 


30E- 


-02 


2 


50-3 


00 


2 


70 


4 


73 


0.42 


17.216 


1 


14E-01 


1 


32E- 


-02 


5 


38E- 


-03 


3 


00-3 


50 


3 


22 


4 


73 


0.44 


17.841 


5 


37E-02 


8 


42E- 


-03 


3 


55E- 


-03 


3 


50-4 


00 


3 


69 


4 


70 


0.40 


16.494 


1 


71E-02 


5 


07E- 


-03 


2 


73E- 


-03 


4 


00-5 


00 


4 


42 


4 


64 


0.43 


17.756 


4 


00E-03 


2 


04E- 


-03 


3 


08E- 


-03 



Table 5.9: (continued) 



p T [GeV] < pr> < M > < x F > < E > Ed 3 a /dp 3 Stat. Error. Syst. Error 









^ 9f1 
o.zu 




<r f\ 9D 
- - u.z,u 










0.00-0.25 


0.16 


5.71 


0.44 


17.995 


3.23E+00 


1.33E- 


-01 


1.15E- 


-01 


0.25-0.50 


0.39 


5.71 


0.44 


17.907 


3.81E+00 


9.07E- 


-02 


1.14E- 


-01 


0.50-0.75 


0.63 


5.72 


0.44 


17.875 


3.10E+00 


6.83E- 


-02 


9.26E- 


-02 


0.75-1.00 


0.87 


5.72 


0.44 


17.881 


2.20E+00 


4.97E- 


-02 


6.17E- 


-02 


1.00-1.25 


1.12 


5.72 


0.44 


17.983 


1.59E+00 


3.86E- 


-02 


4.35E- 


-02 


1.25-1.50 


1.37 


5.72 


0.44 


17.948 


9.18E-01 


2.86E- 


-02 


2.62E- 


-02 


1.50-1.75 


1.61 


5.72 


0.43 


17.844 


5.93E-01 


2.22E- 


-02 


1.77E- 


-02 


1.75-2.00 


1.87 


5.70 


0.43 


17.877 


3.35E-01 


1.64E- 


-02 


1.17E- 


-02 


2.00-2.50 


2.22 


5.73 


0.43 


17.876 


1.62E-01 


8.15E- 


-03 


5.91E- 


-03 


2.50-3.00 


2.71 


5.76 


0.43 


18.006 


5.66E-02 


4.68E- 


-03 


2.66E- 


-03 


3.00-3.50 


3.19 


5.72 


0.44 


18.145 


1.05E-02 


2.38E- 


-03 


9.71E- 


-04 


3.50-4.00 


3.68 


5.67 


0.44 


18.217 


6.48E-03 


1.72E- 


-03 


8.23E- 


-04 


4.00-5.00 


4.35 


5.74 


0.43 


18.107 


1.84E-03 


5.86E- 


-04 


6.55E- 


-04 


6.00-7.00 


6.08 


6.10 


0.40 


17.758 


2.38E-04 


2.51E- 


-04 


4.85E- 


-04 



Table 5.9: (continued) 



p T [GeV] < pr> < M > < x F > < E > Ed 3 a /dp 3 Stat. Error. Syst. Error 













^ 7 on 










0.00-0.25 


0.17 


6.65 


0.44 


18.312 


1.36E+00 


5.98E- 


-02 


4.62E- 


-02 


0.25-0.50 


0.39 


6.64 


0.44 


18.275 


1.50E+00 


4.17E- 


-02 


4.46E- 


-02 


0.50-0.75 


0.63 


6.67 


0.44 


18.292 


1.30E+00 


3.06E- 


-02 


3.62E- 


-02 


0.75-1.00 


0.87 


6.67 


0.44 


18.305 


9.66E-01 


2.35E- 


-02 


2.63E- 


-02 


1.00-1.25 


1.12 


6.68 


0.44 


18.218 


7.17E-01 


1.77E- 


-02 


1.95E- 


-02 


1.25-1.50 


1.37 


6.68 


0.43 


18.205 


4.26E-01 


1.31E- 


-02 


1.21E- 


-02 


1.50-1.75 


1.61 


6.68 


0.43 


18.202 


2.90E-01 


1.03E- 


-02 


9.40E- 


-03 


1.75-2.00 


1.87 


6.70 


0.43 


18.157 


1.57E-01 


7.45E- 


-03 


5.68E- 


-03 


2.00-2.50 


2.21 


6.71 


0.43 


18.244 


7.50E-02 


3.47E- 


-03 


2.60E- 


-03 


2.50-3.00 


2.70 


6.72 


0.43 


18.169 


2.60E-02 


1.97E- 


-03 


1.39E- 


-03 


3.00-3.50 


3.19 


6.71 


0.43 


18.449 


1.08E-02 


1.33E- 


-03 


7.07E- 


-04 


3.50-4.00 


3.69 


6.68 


0.42 


18.088 


3.55E-03 


7.31E- 


-04 


3.86E- 


-04 


4.00-5.00 


4.26 


6.68 


0.45 


19.144 


1.45E-03 


3.46E- 


-04 


2.83E- 


-04 


5.00-6.00 


5.55 


6.58 


0.41 


17.946 


1.75E-04 


1.75E- 


-04 


8.61E- 


-05 



Table 5.9: (continued) 



p T [GeV] < pr> < M > < x F > < E > Ed 3 a /dp 3 Stat. Error. Syst. Error 









7 9D 
i .zu 




- - O. I u 












0.00-0.25 


0.16 


7.84 


0.44 


18.819 


7.58E- 


-01 


3.95E- 


-02 


5.35E- 


-02 


0.25-0.50 


0.38 


7.88 


0.44 


18.765 


9.06E- 


-01 


2.80E- 


-02 


3.11E- 


-02 


0.50-0.75 


0.63 


7.85 


0.44 


18.781 


7.89E- 


-01 


2.02E- 


-02 


2.26E- 


-02 


0.75-1.00 


0.88 


7.85 


0.44 


18.734 


5.83E- 


-01 


1.46E- 


-02 


1.59E- 


-02 


1.00-1.25 


1.12 


7.87 


0.44 


18.731 


4.24E- 


-01 


1.14E- 


-02 


1.15E- 


-02 


1.25-1.50 


1.37 


7.86 


0.43 


18.653 


2.89E- 


-01 


8.56E- 


-03 


7.86E- 


-03 


1.50-1.75 


1.62 


7.89 


0.43 


18.680 


1.65E- 


-01 


6.00E- 


-03 


4.59E- 


-03 


1.75-2.00 


1.87 


7.85 


0.44 


18.779 


1.17E- 


-01 


4.85E- 


-03 


3.51E- 


-03 


2.00-2.50 


2.21 


7.89 


0.44 


18.804 


5.34E- 


-02 


2.24E- 


-03 


1.54E- 


-03 


2.50-3.00 


2.70 


7.87 


0.43 


18.748 


1.95E- 


-02 


1.28E- 


-03 


7.02E- 


-04 


3.00-3.50 


3.19 


7.93 


0.43 


18.890 


6.28E- 


-03 


7.44E- 


-04 


2.66E- 


-04 


3.50-4.00 


3.71 


7.88 


0.43 


18.772 


2.50E- 


-03 


4.32E- 


-04 


1.85E- 


-04 


4.00-5.00 


4.36 


7.88 


0.43 


18.936 


1.17E- 


-03 


2.42E- 


-04 


1.53E- 


-04 


5.00-6.00 


5.26 


7.98 


0.43 


19.251 


2.03E- 


-04 


9.31E- 


-05 


1.21E- 


-04 



Table 5.9: (continued) 



p T [GeV] < pr> < M > < x F > < E > Ed 3 a /dp 3 Stat. Error. Syst. Error 









iu.oo 
















0.00-0.25 


0.17 


11.56 


0.45 


20.901 


4.24E- 


-02 


1.13E- 


-02 


4.63E- 


-03 


0.25-0.50 


0.39 


11.51 


0.44 


20.469 


7.26E- 


-02 


7.40E- 


-03 


3.50E- 


-03 


0.50-0.75 


0.62 


11.58 


0.45 


20.859 


6.28E- 


-02 


5.58E- 


-03 


4.80E- 


-03 


0.75-1.00 


0.88 


11.43 


0.45 


20.832 


4.29E- 


-02 


4.11E- 


-03 


4.31E- 


-03 


1.00-1.25 


1.12 


11.59 


0.45 


20.858 


3.51E- 


-02 


2.98E- 


-03 


2.61E- 


-03 


1.25-1.50 


1.37 


11.47 


0.44 


20.595 


2.36E- 


-02 


2.43E- 


-03 


1.22E- 


-03 


1.50-1.75 


1.63 


11.55 


0.44 


20.653 


1.30E- 


-02 


1.91E- 


-03 


8.16E- 


-04 


1.75-2.00 


1.85 


11.41 


0.44 


20.675 


1.01E- 


-02 


1.53E- 


-03 


1.60E- 


-03 


2.00-2.50 


2.21 


11.52 


0.44 


20.770 


4.41E- 


-03 


5.59E- 


-04 


6.29E- 


-04 


2.50-3.00 


2.71 


11.37 


0.44 


20.761 


1.40E- 


-03 


2.73E- 


-04 


2.65E- 


-04 


3.00-3.50 


3.18 


11.69 


0.45 


21.149 


4.98E- 


-04 


2.34E- 


-04 


3.85E- 


-05 


3.50-4.00 


3.57 


11.34 


0.40 


19.469 


3.51E- 


-04 


2.07E- 


-04 


5.26E- 


-05 


4.00-5.00 


4.23 


11.83 


0.41 


20.241 


6.69E- 


-05 


4.32E- 


-05 


1.62E- 


-05 



12.85 < M M+M - < 16.85 
0.00-0.25 0.15 13.83 0.45 22.325 4.21E-02 1.54E-02 5.61E-03 



Table 5.9: (continued) 



p T [GeV] < pr> < M > < x F > < E > Ed 3 a /dp 3 Stat. Error. Syst. Error 



0.25- 


-0.50 


0.40 


13.34 


0.46 


22.267 


2.18E- 


-02 


9.75E- 


-03 


6.32E- 


-04 


0.50- 


-0.75 


0.64 


13.76 


0.46 


22.449 


1.85E- 


-02 


4.48E- 


-03 


1.08E- 


-03 


0.75- 


-1.00 


0.87 


13.93 


0.45 


22.507 


9.92E- 


-03 


2.88E- 


-03 


5.70E- 


-04 


1.00- 


-1.25 


1.08 


13.33 


0.44 


21.690 


5.75E- 


-03 


2.25E- 


-03 


3.79E- 


-04 


1.25- 


-1.50 


1.37 


13.72 


0.44 


22.071 


8.60E- 


-03 


2.21E- 


-03 


5.42E- 


-04 


1.50- 
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0.45 
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-03 
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0.47 
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0.43 
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0.45 
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3.52E- 
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2.14E- 


-04 


6.36E- 


-05 


3.00- 


-3.50 


3.04 


13.99 


0.51 


24.267 


5.10E- 


-04 


3.61E- 


-04 


9.74E- 


-05 



Table 5.10: Invariant form Ed 3 a/dp 3 (in units of pb / GeV 2 / nucleon) for the deuterium 
cross section over the range 0.55 < xp < 0.8. Statistical and point-to-point systematic 
uncertainties are shown. An additional ±6.5% uncertainty in the normalization is common 
to all data points. 
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Table 5.10: (continued) 



p T [GeV] < pr> < M > < x F > < E > Ed 3 a /dp 3 Stat. Error. Syst. Error 
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Table 5.10: (continued) 
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Table 5.10: (continued) 



p T [GeV] < pr> < M > < x F > < E > Ed 3 a /dp 3 Stat. Error. Syst. Error 
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Table 5.10: (continued) 



p T [GeV] < pr> < M > < x F > < E > Ed 3 a /dp 3 Stat. Error. Syst. Error 
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Table 5.10: (continued) 
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6 DISCUSSION 

In the thirty years since the seminal experiment by Christenson et al. [Chr 70], 
the Drell-Yan process for continuum dilepton production in hadronic interactions 
remains an active area of experimental and theoretical investigation. Several 
experiments have studied both pion- and proton-induced dileptons, providing a 
wealth of data utilizing various beams, targets and center-of-mass energies. When 
these data are analyzed within the framework of perturbative QCD along with 
data from other types of hadronic interactions, they provide a powerful tool for 
extracting information about the parton distribution functions. 

In the preceding chapters we have described the procedures used in and re- 
ported the results of a measurement of continuum diumon cross sections in 800- 
GeV pp and pd interactions. In this chapter, we compare our results to exist- 
ing dimuon data from pA and pd measurements performed by the CERN NA3 
[Bad 83], FNAL E605 [Mor 91] and FNAL E772 [McG 94] experiments. Fer- 
milab E866 has already published [Haw 98, Tow 01] results on the deuterium- 
to-hydrogen cross section ratio, and we will also demonstrate that our absolute 
measurements are consistent with those previous results. We will also discuss the 
level of agreement between our results and next-to-leading order calculations of 
the cross sections based on various sets of parton distribution fits. We will high- 
light some important differences between data and theory, and discuss the likely 
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impact our results will have on future fits to the parton distributions. 
6.1 Comparison with Other Experiments 

Following the work of Christenson et al, several experiments measured Drell-Yan 
cross sections using various beams and targets. The first high-statistics mea- 
surements to study the Drell-Yan process over an extended range of kinematics 
examined dimuon production in pA interactions. CERN NA3 [Bad 83] studied 
pPt at y/s = 27.4 GeV, while Fermilab E605 [Mor 91] studied pCu at y/s = 38.8 
GeV. Less extensive measurements of the pp and pd cross sections in a narrow 
xf range centered around xf = were also carried out at CERN. CERN R-209 
[Ant 82] studied the pp cross sections at y/s = 44 and 62 GeV, while CERN NA51 
[Bal 94] studied the pp and pd cross sections 1 at y/s = 29 GeV. In the early 1990s, 
the pd cross sections were measured over a wide range of mass and xp by the Fer- 
milab E772 collaboration [McG 94]. They published results spanning the range 
0.05 < xf < 0.7 based on an 83, 000 event dimuon sample. It is important to note 
that the E605 and E772 measurements were conducted using essentially the same 
apparatus as E866. Thus, these two experiments provide important cross checks 
on the consistency of our results. 

Figures 6.1 and 6.2 compare the E866 measurements of the scaling form of 

1 CERN NA51 did not publish absolute measurements, only the asymmetry between pp and 
pd cross sections. 
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the cross section M 3 d 2 a/dMdx F in pd interactions to the NA3, E605 and E772 
measurements. The data are plotted versus the scaling variable yfr to enable a 
direct comparison between the results measured at different y/s. The E866 data, 
which are subject to an additional ±6.5% global normalization uncertainty, are 
in good agreement with these previous dimuon experiments over a broad range 
of mass and x F . The E605, E772 and E866 data typically agree to within less 
than the quoted normalization uncertainties 2 of the three experiments. The only 
major area of disagreement is at larger values of x F and smaller masses, where the 
E866 measurement differs significantly from the E772 results. In this kinematic 
region, the E772 measurement is subject to larger systematic uncertainties (~ 
20%) [McG 94]. The major differences between E772 and E866 in this region 
are thought to be primarily due to uncertainties in the E772 magnetic field maps 
[McG]. 

Because the E772 and E866 measurements used essentially the same apparatus 
it is worth exploring this difference in a little more detail. Figure 6.3 shows the 
results from the low- and high-mass data settings used in E866, and the published 
E772 results. They are shown over the range 0.4 < xp < 0.5 where the discrepancy 
is most pronounced. The two E866 mass settings have very different acceptances. 

The acceptance in the high-mass data is plunging rapidly with decreasing mass, 

2 We expect the systematic uncertainties in the normalization of these three experiments to 
be correlated at some level, since they used the same beam line monitors to calibrate the data, 
and some of the same calibration measurements are common. 
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Figure 6.1: Comparison between E866 and previous dimuon experiments. The 
scaling form M 3 d 2 a/dMdxp is plotted versus y/r for the E866 pd, E772 pd 
[McG 94], E605 pCu [Mor 91] and NA3 pPt [Bad 83] measurements. Errors are 
statistical only, except for the E866 uncertainties which are the linear sum of sta- 
tistical and systematic uncertainties. The E866 data contain an additional ±6.5% 
systematic uncertainty in the normalization. 
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Figure 6.2: Comparison between E866 and previous dimuon experiments. The 
scaling form M 3 d 2 'a / ' dM dx p is plotted versus y/r for the E866 pd, E772 pd 
[McG 94], and NA3 pPt [Bad 83] measurements. Errors are statistical only, ex- 
cept for the E866 uncertainties which are the linear sum of the statistical and 
systematic uncertainties. The E866 data contain an additional ±6.5% systematic 
uncertainty in the normalization. 
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while the low-mass acceptance is fairly large in the same region. The agreement 
between these data sets is quite good, suggesting that the E866 data is generally 
free of the systematic uncertainties present in the E772 data in this region. 

Fermilab E772 has also published measurements of the triple-differential cross 
section Ed 3 a/dp 3 , over a more restricted kinematic range (0.1 < xf < 0.3) than 
their double-differential measurements. This range is below the xp range where 
the E772 and E866 measurements disagree, and was imposed on E772 by cor- 
relations between their xp and pp acceptance [McG], which was much less of a 
concern in the E866 data. In the previous chapter we tabulated Ed 3 a /dp 3 in var- 
ious xp bins covering our entire xp range. To make a direct comparison between 
our results and those of E772, we recalculated our cross sections using the E772 
binning. These results are shown in figure 6.4 along with the E772 measurement. 
In general, there is excellent agreement between the two experiments. Only as 
Pt — > GeV do the two experiments disagree. The cross sections in this region 
are highly sensitive to the precise alignment of the beam. Uncertainties in either 
or both the beam position and angle would lead to systematic shifts in the pp of 
each event. Thus, the differences between the two measurements may be due to 
beam alignment issues in the two experiments. 
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Figure 6.3: Detailed comparison between E866 and E772 in the range 0.4 < 
xf < 0.5. Deuterium results from the low- and high-mass E866 data sets are 
shown along with the E772 measurement. The error bars represent statistical 
uncertainties only. The high-mass E866 data has been offset by 0.05 GeV. 
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Figure 6.4: Comparison of the E866 and E772 results for Ed 3 a/dp 3 . Data are 
in the range < xf < 0.3. Open (closed) circles are the E866 results for pp (pd). 
Open diamonds represent the E772 results. The E866 error bars are the linear 
sum of the statistical and systematic uncertainties. 
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6.2 Cross Section Ratios 

We turn now to a comparison of our absolute measurements with the previously 
published cross section ratios. The primary goal of Fermilab E866 was to measure 
the deuterium-to-hydrogen cross section ratio as a function of x 2 , which provides 
a clean probe of the Bjorken-x dependence of the ratio of anti-down to anti- 
up quarks in the nucleon sea. 3 Those results [Haw 98, Tow 01] have already 
been published and incorporated into the various global analyses of the parton 
distribution functions [CTEQ, MRST, GRV]. 

Although the previous analysis and this current work began with the same 
data sample, the tighter event selection criteria required for absolute normalization 
resulted in the elimination of much of the data above x 2 ~ 0.15. Direct comparison 
between the results presented in chapter 5 and those of Towell et al. [Tow 01] are 
complicated by this. Further complicating matters, the results of Towell et al. are 
differential only in a single variable (x 2 ) , 4 while our new results are differential in 
the mass and xp of the dimuon pair. 

In figure 6.5, we present the deuterium-to-hydrogen cross-section ratio versus 

mass in four xf bins. The ratios were calculated by integrating the data in tables 

3 Although not explicitly shown in chapter 1, the ratio of the leading-order pd and pp cross 

sections can be written for large Xf as — (M, ij?;^) ~ j — \ + d z , " rrr — z\ (l + =(#2) ), where 

ct pp U+4j(ii)j(j;2)J V w J 

xi is evaluated according to equation 1.3. 

4 Towcll et al. did look at the cross section ratio versus other variables, but none which are 
directly comparable to our double-differential cross sections. 
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5.3 through 5.10 over all px and taking the ratio of deuterium to hydrogen in each 
bin. Next-to-leading order calculations of the cross section ratio, based on the 
CTEQ 5 and CTEQ 6 5 parton distributions, are also shown in the figure. Since 
the E866 cross section ratios were used in the CTEQ 5 and CTEQ 6 global anal- 
yses, they provide a convenient means to check the consistency between the two 
analyses. We also superimpose the previous cross section ratio results of Towell 
et al. on figure 6.5. Although those results are differential only in x 2 , reference 
[Tow 01] does quote a mean mass and x F in each x 2 bin. We therefore plot the 
cross section ratios in each x 2 bin at the mean mass quoted in the reference, and 
only in the xf bins appropriate to their mean xf- Overall, the level of agreement 
between the data presented in chapter 5 and the cross section ratios previously 
published by E866 is quite good. 

To provide a more detailed and direct comparison with the previously pub- 
lished results, we calculated the absolute cross section d 2 a /dx\dx 2 for hydrogen 
and deuterium using the same x 2 binning as reference [Tow 01] and the same pro- 
cedures as outlined in chapter 4. For this study, the event selection criteria were 
loosened to admit the large-a^ events which were otherwise problematic in the ab- 
solute cross sections analysis. 6 Figure 6.7 shows the deuterium-to-hydrogen cross 

5 The MRST 98, MRST 2001 and GRV 98 parton distributions give similar results. 

Specifically, these events tended to occupy regions of the spectrometer with cither known 
efficiency problems or uncertainties in reconstructing the event kinematics. These effects should 
almost completely cancel in ratio. We did not include any data from intermediate mass data sets 
3 and 4, which account for approximately 1/3 of the total number of events used in reference 



197 



section ratio versus x<i as calculated in this analysis in three different x\ bins. As 
with the previous analysis, there appears to be no significant x\ dependence on 
the cross section ratio. Averaging the ratios over all x\ we obtain a cross section 
ratio versus x 2 only which is compared with the results of Towell et al. in figure 
6.7. The only differences between the two analyses are almost completely due to 
the larger event sample available to reference [Tow 01]. 

6.3 Data versus Theory 

Detailed information on how our data will impact the parton distribution functions 
must inevitably await their incorporation into future global analyses. It is still 
useful to compare our results with theoretical calculations of continuum dimuon 
cross sections. Such comparisons will demonstrate whether the existing sets of 
parton distributions are consistent with our data in kinematic regions already 
constrained by previous measurements (E605), and reveal where we expect our 
data to substantially impact upon the parton distributions in the next round of 
fits. 

As discussed back in chapter 1, continuum dimuon production is sensitive to 
the gluon distributions through QCD processes at and beyond next-to-leading- 
order. Dimuon production is dominated by gluon Compton scattering for pt > 

\M [Ber 98], and the E866 results satisfy this kinematic condition over an ex- 
[Tow 01]. 
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Figure 6.5: Ratio of pc? and (per nucleon) cross sections. The results of this 
analysis (circles) are compared with the results (triangles) of the previous analysis 
of Towell et al. [Tow 01]. Although the results of Towell et al. are differential 
only in x 2 j the mean mass and xp of the events in each £ 2 bin was tabulated in the 
reference. We plot those results at those mean mass values, and only in bins 
appropriate to the mean xp. NLO calculations of the ratio based on the CTEQ 5 
and CTEQ 6 parton distributions are also shown. 
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Figure 6.7: Ratio of pd and pp (per nucleon) cross sections plotted versus X2, at 
the center of each x<i bin. The data are compared with the previous E866 results, 
which are plotted at the mean X2 in each bin. Both analyses are subject to an 
additional ±1% systematic uncertainty in the relative normalization. 
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tended range in mass and Xf- Thus, in addition to information on the light- 
antiquark distributions, our data should provide a new and valuable source of 
information about the gluon distributions. Because this is still an area of ongo- 
ing theoretical investigation, in the following discussion we limit ourselves to a 
discussion of the quark and antiquark distributions. 

6.3.1 Scaling Form: M 3 d 2 a/dMdx F 

In the preceding chapter, the scaling form M 3 d 2 a/dMdx F was tabulated and 
plotted in mass and x F bins. Figures 5.1 through 5.16 plotted the pp and pd 
cross section in x F bins versus the invariant mass of the dimuon pair. In each xf 
bin we also compared the data to perturbative QCD calculations performed to 
next-to-leading order. These calculations were based on several of the currently 
available and most commonly used sets of fits to the parton distributions. In 
general, they reproduce the shape and normalization of the data well, save for the 
GRV 98 distributions which systematically overestimate the data by ~ 20%. To 
better illustrate this, the theoretical calculations were fit to the data to determine 
a "If -factor" , which we define as the ratio of the experimentally measured dimuon 
continuum to NLO theory. 

The if-factors for hydrogen and deuterium are shown in table 6.1. All of the 
parton distributions studied, with the exception of the GRV 98 set, reproduce the 
normalization of the measured cross sections to better than our systematic uncer- 
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Table 6.1: if'-factors, where K' = -^nu-, for 800-GeV pp and pd dimuon produc- 
tion. The errors represent the statistical and point-to-point systematic uncertain- 
ties in the data. The data are subject to an additional ±6.5% uncertainty in the 
normalization. 
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X 2 /ndf 


CTEQ5 


0.9753 ±0.0037 
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0.9633 ±0.0017 


2.56 


CTEQ6 


1.0150 ±0.0028 


1.44 


1.0003 ±0.0004 
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MRST 98 


0.9732 ±0.0022 


1.43 


0.9604 ±0.0028 


2.42 


MRST 2001 


0.9799 ±0.0034 


1.50 


0.9664 ±0.0078 


2.48 


GRV 98 


0.8107 ±0.0043 


2.11 


0.8075 ±0.0002 


4.24 



tainties. This is not unexpected. The valence distributions are well constrained 
by DIS measurements, and the E605 Drell-Yan data (with which we agree) are 
an important input into the fits, constraining the antiquark sea over much of the 
range covered by E866 [CTEQ, MRST, GRV]. 

It is interesting to note the improvement in the normalization of the CTEQ 
6 [CTEQ 02] based calculations relative to CTEQ 5 [CTEQ 99]. The main dif- 
ference between the CTEQ 6 distributions and the previous CTEQ fit is a stiffer 
large-x gluon, the result of including Tevatron jet data in the newer analysis 
[CTEQ 02]. The same data has also been incorporated into the latest MRST fit 
[MRST 01], but without as significant a change in the gluon distributions. Though 
there is improvement in the calculations based on the MRST 2001 distributions 
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over the previous MRST analysis, it is not as dramatic as that seen in the CTEQ 
distributions. 

The relatively poor quality of the fits shown in table 6.1 can be explained by 
behavior observed in the cross sections which was not anticipated by the parton 
distribution functions. The agreement between data and theory is quite good 
below xp ~ 0.3, which is not surprising given the degree of agreement with the 
E605 data previously discussed. As x? increases beyond the range constrained 
by E605, the agreement between the E866 data and the NLO calculations deteri- 
orates, especially at low mass. While this kinematic range probes smaller values 
of X2 than reached in E605, it also probes values of x\ (the momentum fraction 
of the annihilating quark in the parton model) where the estimated uncertainties 
on the valence quark distributions are about ~ 20%. Thus, we expect our data 
to provide important constraints on both the valence and sea quarks in the next 
round of fits. 

6.3.2 d 2 a / dx\dx2 

To explore the impact of our data on future fits in more detail, we now examine the 
cross sections written in the form d 2 a/dxidx 2 - While higher-order QCD processes 
confuse the interpretation of x\ and x 2 here, one should recall that almost half 
of the cross section is due to the leading-order qq annihilation for which x\ and 
X2 have their usual meaning. Writing the pp and pd cross sections in the limit of 
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large i^we have 




oc -u(xi)u(x 2 ) + -d(x 1 )d(x 2 ) 



(6.1) 



dx\dx2 



Thus, the shape of the cross sections relative to theoretical calculations should 
help us understand the effect of our data on the parton distributions. 

In the comparisons between data and theory which follow, the cross sections 
of the form d 2 o jdx\dx 2 were determined using the same procedures outlined in 
previous chapters. However, certain limitations are imposed on this analysis by 
the Monte Carlo event sample, which was generated over a limited range in mass 
and xf- Because of this, calculating the acceptance in the lowest-lying x±-x 2 bins 
is problematic - the x±-x 2 bins extend beyond the kinematic limits of the Monte 
Carlo sample. This mainly affects the data below x 2 = 0.03, which partially 
overlaps the low-mass, large-a^ data discussed in the previous section. 

To study the likely impact of our data on future fits to the antiquark dis- 
tributions, we examine the behavior of data and theory as a function of x 2 , the 
momentum fraction of the annihilating antiquark in the parton model. After the 
cross sections differential in the momentum fractions of the interacting partons 
was determined, the ratio of those cross sections to next-to-leading-order calcula- 
tions based on the CTEQ 6 [CTEQ 02] partons was calculated in each X\-x 2 bin. 
Figure 6.8 shows these ratios in two x 1 bins, each containing roughly the same 
number of events. The difference in shape (if any) between data and theory versus 
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Figure 6.8: Ratio of data to NLO calculations based on CTEQ 6 [CTEQ 02] 
versus x 2 - The left panel shows the range 0.225 < x\ < 0.375, and the right panel 
shows the range 0.375 < X\ < 0.775. 

X2 appears to be consistent over the two ranges in x\. 

In figure 6.9 we average together the ratios of data and theory over all Xi bins 
to obtain the ratio of data to theory as a function of x 2 only. The theoretical 
calculations clearly reproduce our data well, to within the quoted statistical and 
systematic errors in the data, over the x 2 range covered. This indicates that the 
CTEQ 6 fits to the light antiquarks are consistent with our data. Since these fits 
are largely constrained by the E605 Drell-Yan measurements, this result was not 
unexpected. 



The same analysis can be carried out to examine what our data imply for 
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future fits to the valence quark distributions. Figure 6.10 shows data-divided- by- 
theory plotted versus xi, over two x 2 ranges each representing about half of the 
available statistics. There is an x± dependence which is consistent in both x 2 bins. 
Combining the results over all x 2 in figure 6.11 we see that at large x\ the data 
falls about 20% below the theoretical predictions. This seems to suggest that 
the valence quark distributions are smaller at large Bjorken-x than in current 
parton distributions, though not outside the estimated uncertainties of the fits 
[CTEQ 02]. Parton distribution fits to DIS data at next-to-next-to- leading-order 
have also recently become available [MRST 02], and they also suggest a reduced 
valence distribution, relative to current NLO fits. 

6.4 Future Prospects 

Continuum dimuon production has been an active area of reseach for over thirty 
years, and the outlook for future measurements remains bright. Based on the 
cross section ratios published by E866 [Haw 98, Tow 01], Fermilab has accepted 
a proposal [E906] for a new dimuon experiment, Fermilab Experiment 906, to be 
conducted using the 120-GeV proton beam from the Main Injector. Like E866, the 
E906 measurements will examine the differences in continuum dimuon production 
in pp and pd interactions. Their primary goal is to measured out to x ~ 0.5. 
The extended range in x is made possible in large measure due to the seven-fold 
increase in cross section using 120-GeV protons relative to 800-GeV protons. 
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Flavor asymmetry is but one of the topics to be adressed with the new exper- 
iment. Absolute measurements of the hydrogen and deuterium cross sections will 
also be possible, greatly extending the kinematic range and statistical precision 
available to constrain the parton distributions. Furthermore, these measurements 
may also be able to demonstrate an important property of QCD. Logarithmic 
scaling violations have long been established in DIS measurements, but have not 
been conclusively shown using the Drell-Yan process. The reasons for this include 
the relatively high q 2 probed in dimuon experiments, and the partial cancelation 
of the effect between the large-x behavior of the valence quarks, and the low-rr 
behavior of the antiquarks [Fre 90, McG 99]. Given the statistical precision and 
overall understanding of systematic uncertainties attained with the E866 results, 
comparisons with the future measurements of E906 may well be able to demon- 
strate scaling violations in the Drell-Yan process. 

6.5 Conclusion 

Fermilab Experiment 866 has performed an absolute measurement of dimuon cross 
sections in 800-GeV pp and pd interactions. These data represent the most ex- 
tensive measurement of the Drell-Yan process to date, covering a wider range 
of kinematics with better statistical precision than has been achieved in previ- 
ous measurements. The results are in good agreement with existing pA and pd 
measurements, as well as NLO calculations based on the various sets of parton 
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distribution functions. The overall level of agreement between data and theory- 
indicates that the antiquark distributions are well understood over the kinematic 
range relevant to E866. However, differences between data and theory suggest 
that the valence distributions are overestimated at large x in current parameteri- 
zations. Given the statistical precision of the data, we can expect a major impact 
on future fits to the valence quarks. 
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